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This  work  represents  an  analysis  of  the  thermal  balance  of  a 
fusioning  plasma  from  a  control  system  perspective.   By  applying  the 
techniques  of  classical  control  theory,  the  requirements  for  stability 
and  the  response  characteristics  of  the  thermal  balance  have  been 
evaluated.   The  results  show  that  open-loop  equilibria  are 
characterized  by  restrictively  narrow  stable  operating  temperature 
regimes  and  generally  poor  system  performance.   The  results  also  show 
that  closed-loop  proportional  feedback  using  fuel  feedrate  and 
injection  energy  can  be  used  to  extend  the  stable  operating  temperature 
regime  and  significantly  improve  the  system  response.   Additional 
improvement  is  available  using  proportional-integral-derivative  (PID) 
control. 

The  stabilized  temperature  regime  is  a  function  of  the  governing 
confinement  law,  the  injection  energy  of  incoming  fuel,  and  the 
feedback  coefficients.   A  broad  range  of  feedback  fractions  is 


available  with  which  to  stabilize  the  fusion  power  balance.   The 
feedback  can  thus  be  chosen  so  as  to  improve  the  system  response.  With 
constant  confinement,  sluggish  open-loop  performance  with  no  injection 
energy  can  be  moderately  improved,  and  the  high  overshoots  present  with 
75  kev  injection  energy  can  be  significantly  reduced.   These 
improvements  are  possible  only  at  temperatures  above  20  keV  and  require 
high  feedback  fractions.   With  Bohm  confinement,  high  open-loop 
overshoots  present  with  zero  injection  energy  can  be  reduced  to 
acceptable  levels  at  temperatures  as  low  as  20  keV,  with  a  simultaneous 
decrease  in  settling  times  to  under  30  seconds.   With  75  keV  injection 
energy,  acceptable  overshoot  can  be  obtained  at  temperatures  as  low  as 
10  keV,  with  the  time-to-peak  below  20  seconds  and  settling  times  less 
than  30  seconds.   Even  with  Bohm  confinement,  it  is  still  difficult  to 
simultaneously  satisfy  overshoot  and  speed  of  response  requirements  at 
low  temperatures  and  low  feedback  fractions. 

Though  the  system  response  is  significantly  improved  with  only 
proportional  feedback,  additional  improvement  is  desirable. 
Preliminary  indications  of  the  benefits  of  PIO  control  are  that  more 
efficient  control  is  possible  while  simultaneously  improving  system 
response. 


CHAPTER  1 
INTRODUCTION 


1.1.   Motivation 


The  magnetic  confinement  of  a  plasma  is  presently  the  leading 
approach  to  achieving  the  conditions  necessary  to  sustain  thermonuclear 
fusion  reactions.   The  confined  plasma,  consisting  of  ionized  deuterium 
and  tritium  fuel,  is  heated  to  a  sufficiently  high  temperature  that 
fusion  events  occur,  releasing  energy  in  the  form  of  neutrons  and 
charged  particles.   Though  the  scientific  feasibility  of  fusion,  i.e., 
that  the  energy  released  is  greater  than  the  energy  used  in  initiating 
the  fusion  events,  has  yet  to  be  demonstrated,  there  is  strong  belief 
in  the  scientific  community  that  this  proof  will  come  within  the  next 
two  years  at  the  Tokamak  Fusion  Test  Reactor  (TFTR)  presently  in  its 
final  year  of  construction  at  Princeton,  New  Jersey. 

Beyond  scientific  proof-of-principle ,  fusion  will  become  a 
practical  energy  source  when  the  engineering  and  economic  issues  of 
producing  net  energy  from  magnetic-confinement  fusion  are  resolved. 
Foremost  among  the  issues  of  scientific  feasibility  is  the  question  of 
plasma  stability.   A  magnetically  confined  plasma  is  out  of 
thermodynamic  equilibrium,  but  attempts  to  relax  toward  equilibrium  via 
Coulomb  collisions  and  collective  plasma  instabilities.   If  Coulomb 
collisions  were  the  only  equilibration  mechanism,  viable  reactor 
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designs  would  be  possible.   However,  plasma  Instabilities  associated 
with  magnetic  confinement  represent  sufficiently  enhanced  energy  loss 
modes  as  well  as  mechanisms  for  undesirable  system  behavior  that 
control  is  essential  in  order  to  demonstrate  the  engineering 
feasibility  of  fusion. 

Plasma  instabilities  can  be  categorized  into  three  groups: 
macroinstability,  microinstability,  and  thermal  instability.   The 
raacroinstabilities  are  often  referred  to  as  MHD  instabilities  because 
they  are  derivable  from  magnetohydrodynamic  (MHD)  fluid  models  of  the 
plasma.   The  MHD  instabilities  result  in  gross  movement  of  the  plasma 
and  rapid  loss  of  confinement  on  the  microsecond  time  scale.   These 
instabilities  have  received  major  theoretical  treatment  and  fortunately 
can  be  controlled  by  proper  design  of  the  magnetic  field  topology, 
machine  configuration,  plasma  positioning  circuitry,  and  by  proper 
selection  of  the  machine  operating  parameters.   Present  theory  predicts 
the  existence  of  magnetic  configurations  of  the  Tokamak  type,  with 
operating  parameters  permitting  an  attractive  reactor,  that  are 
absolutely  stable  to  all  the  important  macroinstabilities  [1]. 

The  microinstabilities  are  not  detrimental  to  overall  plasma 
confinement  as  is  the  case  with  the  macroinstabilities.   Instead,  these 
instabilities  represent  rapid,  small-scale  plasma  transport  processes 
which  lead  to  losses  which  are  in  excess  of  classical  predictions. 
Despite  the  abundance  of  measured  data  and  sophisticated  diagnostic  and 
analysis  tools,  the  mechanism  for  many  of  these  enhanced  diffusive 
losses  remains  unexplained.   The  resultant  diffusion  is  thus  generally 
referred  to  as  "anomalous  transport."  The  net  effect  of 
microinstabilities  is  a  degradation  of  plasma  confinement,  increasing 
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the  difficulty  of  achieving  engineering  break-even.  Much  of  the 
present  experimental  effort  in  fusion  plasma  physics  is  devoted  to 
resolving  the  issues  arising  from  the  microinstabilities . 

The  subject  of  this  dissertation  is  the  control  of  the  thermal 
instability.   This  instabilitv  is  associated  with  considerations 
relating  to  the  behavior  of  the  global  power  balance  as  determined  by 
plasma  mass  and  energy  conservation  considerations.   In  particular,  the 
thermal  instability  refers  to  the  proclivity  of  the  plasma  temperature 
to  either  rise  uncontrollably  to  a  high  temperature,  or  extinguish  the 
fusion  burn,  at  the  onset  of  perturbations  from  the  desired  operating 
conditions  which  may  arise  from  system  and  component  disturbances. 

The  temperature  above  which  the  plasma  energy  balance  is  naturally 
stable  to  perturbations  in  the  parameters  governing  the  thermal 
equilibrium,  i.e.,  above  which  it  does  not  require  external  control,  is 
defined  as  the  critical  temperature.   Control  of  the  energy  balance  to 
ensure  thermal  stability  is  of  importance  to  the  field  of  thermonuclear 
fusion  because  the  operating  temperature  range  which  would  result  in 
peak,  power  density  in  the  plasma  has  been  found  to  be  unstable,  i.e., 
it  is  below  the  critical  temperature.   Operation  at  or  near  peak  power 
density  is  desirable  because  fusion  is  expected  to  be  a  highly  capital 
cost  intensive  energy  source  with  a  significant  sensitivity  of  the 
plant  cost  to  the  plasma  power  density.   In  addition  to  the  economic 
penalty,  operation  at  high  temperatures  introduces  engineering 
feasibility  issues.   High  operating  temperatures  require  materials  and 
components  operating  at  or  beyond  the  limits  of  present  day  or  even 
near  term  technology.   In  addition  to  the  economic  and  engineering 
issues,  control  over  the  thermal  balance  is  required  to  accommodate 
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load  changes,  input  and  system  perturbations,  and  component 
malfunctions . 


1.2.   Previous  Work 

In  contrast  to  the  macroinstabilities,  which  have  been  studied 
since  the  1950's,  and  the  microinstabilities ,  which  have  been  studied 
since  the  1960's,  the  thermal  instability  was  not  recognized  until  the 
1970's.   The  first  considerations  of  the  thermal  instability  and  its 
control  were  provided  by  Mills  [2,3],  and  Ohta  et^  al .  [4].   Using 
dynamic  computer  simulations  of  the  mass  and  energy  conservation 
equations,  Mills'  early  papers  demonstrated  that  if  the  equilibrium 
plasma  temperature  was  below  a  certain  value,  later  termed  the  critical 
temperature  by  Ohta  et^  al_.  [4],  perturbations  in  the  parameters 
governing  the  power  balance  would  lead  to  a  thermal  runaway  in  which 
the  plasma  temperature  would  either  decrease  and  extinguish  the 
reaction  rate,  or  increase  to  a  higher,  thermally  stable  but 
undesirable  value. 

This  behavior  can  be  explained  as  follows.   The  steady  state  of 
the  thermal  balance  is  determined  by  the  careful  balance  of  the  energy 
production  rate  via  fusion  and  the  energy  loss  mechanisms  of  optical 
radiation,  conduction,  and  diffusion  (convection).   In  the  case  of 
thermally  unstable  equilibria,  a  perturbation  resulting  in  an  increase 
in  the  plasma  temperature  results  in  an  increase  in  the  energy 
production  rate  which  is  greater  than  the  increase  in  the  energy  loss 
rates,  which  in  turn  results  in  a  further  increase  in  the  temperature, 
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etc.   As  the  temperature  continues  to  rise,  fusion  power  production 
eventually  decreases,  whereas  the  sum  of  the  loss  terms,  particularly 
radiation,  increases,  until  a  new,  stable  equilibrium  is  established  at 
a  temperature  much  higher  that  the  original.   An  analogous  scenario  for 
an  initial  decrease  in  the  plasma  temperature  can  be  constructed 
wherein  the  final  temperature  extinguishes  the  plasma  burn. 

To  obtain  thermal  stability,  Mills  suggested  several  methods  by 
which  control  might  permit  operation  below  the  critical  temperature. 
Though  not  completely  satisfactory,  these  early  analyses  were  very 
useful,  and  provided  a  platform  for  the  work  of  Ohta  et_  al .  [4]  In 
their  work,  Ohta  et  al.  [4]  supplemented  dynamic  simulations  with  a 
perturbation  analysis  of  the  nonlinear  particle  and  energy  balance 
equations  to  analytically  predict  and  verify  the  stability  of  an 
equilibrium.   They  quantified  the  critical  temperature  for  a  number  of 
confinement  laws  and  were  the  first  to  suggest  feedback  stabilization 
for  operation  in  the  unstable  temperature  regime.   Their  work  still 
continues  to  serve  as  a  standard  against  which  improved  analyses  are 
compared. 

After  the  work  of  Mills  [3]  and  Ohta  et_  al .  [4],  the  thrust  of 
research  was  directed  at  improving  the  representative  models  of  the 
plasma  processes  and  extending  the  analyses  to  additional  control 
mechanisms.   Further  details,  and  representative  models  for  the  plasma 
processes  and  terms  contributing  to  the  energy  balance  are  presented  in 
Chapter  2,  and  an  extensive  review  of  the  literature  on  thermal 
stability  is  presented  in  Appendix  A. 
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1.3.   Control  for  Fusion  Thermal  Stability 

As  the  literature  review  of  Appendix  A  shows,  a  number  of  control 
mechanisms  have  been  proposed,  implemented  and  analyzed  for  the  fusion 
thermal  balance  that  have  permitted  operation  in  otherwise  unstable 
temperature  regimes.   By  the  nature  of  the  methods  used  to  obtain  them, 
the  types  of  mechanisms  explored  have  been  hit-or-miss  in  character, 
have  left  gaps  in  the  system  analysis,  and  although  some  have  achieved 
the  goal  of  stabilizing  the  thermal  balance,  none  has  adequately 
addressed  the  issues  of  overall  system  response.   Certainly  stability 
is  a  necessary  condition  for  operation,  but  perturbations  in  component 
parts  and  systems  disturbances  will  additionally  require  that  the 
integrated  system  possess  certain  desirable  characteristics  and 
properties,  such  as  acceptable  overshoot,  adequate  speed  of  response, 
etc.   A  control  systems  approach  to  the  problem  of  thermal  stability  is 
necessary,  and  a  knowledge  of  the  system  behavior  is  required  in  order 
to  effect  suitable  control. 

The  work  described  herein  is  the  first  formal  classical  control 
theory  formulation  of  the  thermal  stability  problem  which  extends  the 
application  of  control  theory  beyond  a  simple  eigenvalue  or  roots  of 
the  characteristic  equation  analysis.   The  present  formulation  is  able 
to  duplicate  nearly  all  previous  work  regarding  the  determination  of 
the  critical  temperature  for  the  onset  of  the  thermal  instability,  with 
generally  good  agreement.   In  addition  to  determining  the  stability  of 
the  system,  the  extended  analysis  performed  on  the  plasma  model  has 
resulted  in  the  first  evaluation  of  the  response  characteristics  of  the 
thermal  balance  of  a  fusioning  plasma.   Thus,  the  expected  temporal 


-7- 
behavior  of  the  plasma  following  both  normal  system  changes  and 
external  disturbances  can  be  assessed  and  compared  to  reference  systems 
using  standard  inputs  and  accepted  weighting  criteria. 

A  model  of  the  fusioning  plasma  thermal  balance  was  developed 
which  incorporates  improvements  over  many  of  the  models  found  in  the 
literature  while  still  permitting  physical  interpretation  of  the 
results.   Included  were  the  major  processes  of  fusion  burnup  and  energy 
production,  bremsstrahlung  radiation  emission,  and  particle  and  energy 
diffusion.   Processes  which  would  be  present  in  a  fusioning  plasma  but 
not  expected  to  contribute  to  the  question  of  thermal  stability,  such 
as  finite  energy  transfer  mechanisms  among  particle  species,  impurity 
radiation  and  synchrotron  radiation  were  not  included. 

The  selected  plasma  model  was  cast  in  the  form  of  a  classical 
second  order  system  with  a  finite  zero,  and  analytical  expressions  for 
the  system  characteristics  and  specifications  were  derived  in  terms  of 
the  plasma  conditions.   Physical  associations  were  drawn  between  the 
processes  present  in  a  fusioning  plasma  and  the  system  characteristics 
and  transfer  function  coefficients.   Expressions  were  derived  for  the 
system  characteristics  of  gain,  the  finite  zero,  the  damped  and 
undamped  natural  frequencies,  the  damping  ratio,  and  the  damping 
factor.   Expressions  were  derived  for  the  impulse  and  step  response,  as 
well  as  the  system  specifications  of  time-to-peak.,  peak  overshoot, 
settling  time,  bandwidth,  resonance  frequency,  resonance  peak,  and 
phase  margin.   These  were  evaluated  for  the  cases  of  zero  and  75  keV 
injection  energy,  assuming  both  constant  and  Bohm  confinement  laws. 
The  injection  energy  case  was  examined  because  this  is  currently  the 
leading  form  of  plasma  heating  [1].   Bohm  confinement  was  evaluated 


because  Che  particle  and  energy  transport  in  the  reactor  regime  is 
expected  to  exhibit  a  Bohra-type  of  dependence  [5].   The  analytical 
expressions  so  derived  were  used  to  provide  the  initial  evaluation  of 
the  system  characteristics  and  specifications  of  the  open-loop  plasma 
thermal  balance.   Similarly,  expressions  for  the  closed-loop  transfer 
function  using  proportional  feedback  were  derived  and  evaluated.   The 
results  were  used  to  determine  the  feedback  fractions  necessary  to 
guarantee  stability  and  acceptable  system  performance.   Feedback  using 
advanced  control  techniques  was  also  investigated. 

The  format  of  this  study  on  the  control  of  a  fusioning  plasma  to 
achieve  thermal  stability  and  acceptable  system  response  proceeds  in  a 
straightforward  manner.   The  next  chapter  presents  the  assumptions  made 
in  modeling  the  plasma  thermal  balance  and  describes  the  range  of 
applicability.   It  also  develops  the  system  transfer  functions  that  are 
used  in  the  open-  and  closed-loop  analysis  of  Chapters  3  and  4. 
Open-loop  characteristics,  such  as  natural  frequency  and  damping,  and 
the  system  specifications,  such  as  peak  overshoot  and  settling  time,  of 
the  fusioning  plasma  are  presented  in  Chapter  3.   The  analogous 
closed-loop  behavior  is  presented  in  Chapter  4.   Chapter  5  uses  the 
analyses  and  results  of  the  previous  chapters  to  develop  criteria  for 
the  desired  plasma  behavior  and  implements  advanced  control  techniques 
to  further  improve  the  system  response.   The  conclusions  of  this  study 
and  the  recommendations  for  future  work  are  also  presented  in  Chapter 
5.      The  background  information  on  the  major  topics  of  this  dissertation 
is  contained  in  the  Appendices.   Appendix  A  presents  an  introduction  to 
the  aspects  of  fusion  which  are  central  to  the  issues  of  the  thermal 
instability.   It  also  presents  a  comprehensive  literature  review  of  the 
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field  of  plasma  thermal  stability.   Appendix  B  presents  a  review  of 
control  theory.   And  Appendix  C  derives  analytical  expressions  for  the 
system  specifications  of  the  transfer  function  describing  the  plasma 
thermal  balance. 


CHAPTER  2 
PLASMA  MODEL  AND  TRANSFER  FUNCTION  DEVELOPMEN" 


2.1.   Model  Assumptions 

The  governing  dynamics  of  the  thermal  balance  of  a  magnetically 
confined  fusioning  plasma  are  aptly  described  by  a  zero-dimensional 
formulation  of  the  particle  and  energy  conservation  equations.   The 
application  of  certain  simplifying  assumptions  to  these  nonlinear 
differential  equations  transforms  them  from  a  complex  collection  of 
mathematical  terms  into  a  workable  set  of  equations.   The  usefulness  of 
the  latter  rests  on  their  simplicity,  permitting  their  analysis  and 
optimization  via  classical  control  methods.   Tn  addition,  after 
linearization  these  streamlined  equations  yield  plasma  transfer 
functions  whose  coefficients  have  explicit  physical  associations. 
Temporal  response  characteristics  can  then  be  directly  attributed  to 
the  various  plasma  parameters.   The  equations  are  now  presented  and  the 
simplifying  assumptions  described.   The  development  of  the  plasma 
transfer  functions  then  follows. 

The  equations  of  particle  and  energy  conservation  for  a  fusioning 
plasma  are  given  by 


2 

—  =  S  +  X  (kT  -  kT)  -  — <DT>  -  —  2.1 

dt        s   o         2        t 
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and 


Hlf!kT  =    [S   +  Kg(kTQ   -    kT)][Tg   +   KT    (kT0   -   kT) 


dt 


2 
+  Q  !L-<VT>   -  bRn2(kT)l/2  -  3nkT  2>2 


lE 


where 

n  =  fuel  particle  density  [cm   ] 
kT  =  plasma  temperature  [keV] 

kT  =  desired  operating  temperature  [keV] 

3   —1 
<DT>  =  reactivity  for  the  deuterium-tritium  reaction  [cm  sec   ] 

_o     _  1 

S  =  fuel  injection  rate  [cm  sec   ] 

K  =  fuel  injection  rate  feedback  gain 

T  =  fuel  injection  energy  [keV] 

KT  =  fuel  injection  energy  feedback  gain 
s 

T  =  particle  confinement  time  [sec] 

Ty   =  energy  confinement  time  [sec] 

Q  =  alpha  particle  energy  =  3520  keV,  and 

bR  =  bremsstrahlung  emission  constant  =  3.36  x  10    cm  keV1/,2sec 

The  terras  on  the  right-hand-side  of  equation  2.1  represent  (in 
order)  particle  injection  from  an  external  source,  feedback  injection 
from  an  external  source,  fuel  ion  consumption  via  fusion  reactions,  and 
diffusion  leakage.   In  equation  2.2,  the  right-hand  terras  represent 
power  input  to  the  plasma  from  external  particle  injection  (including 
the  feedback  contributions),  power  gain  from  fusion  reactions, 
bremsstrahlung  radiation  power  loss,  and  lastly,  nower  loss  through 
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dif fusion. 

The  diffusion  of  particles  from  a  magnetically  confined  plasma 
remains  the  largest  source  of  uncertainty  in  current  fusion  research. 
The  basic  mechanism  for  particle  loss  in  a  linear  magnetic  confinement 
system  is  classical  diffusion.   Classical  diffusion  losses  are  the 
result  of  binary  Coulomb  collisions.   Collisions  cause  particles  to 
step  their  way  across  magnetic  field  lines  with  a  characteristic  step 
size  given  by  the  particle  orbits'  gyroradii.   In  a  toroidal  device, 
the  characteristic  step  size  is  modified  by  the  spatial  variation  of 
the  magnetic  configuration.   The  resulting  diffusion  is  called 
neoclassical.   Experimental  results,  however,  have  demonstrated  that 
particle  transport  in  toroidal  machines  is  dominated  by 
microinstabilities.   These  experimentally-determined  transport  laws 
have  been  labeled  pseudoclassical  diffusion.   Anomalous  behavior  still 
predominates  and  for  this  reason,  confinement  in  this  analysis  is 
modeled  with  general  expressions  of  the  form 

te  =  xcn*(kT)ra  2.3 

and 

t„  =  R^E  2'4 

where  t  ,  I,   m,  and  R  are  confinement  constants  characterizing  the 
desired  confinement  scheme.   The  schemes  investigated  herein  include 
constant  (l   =  m  =  0)  and  Bohm  (£=0ra=-l)  confinement. 

These  four  equations  comprise  the  plasma  model  used  in  the  initial 
analysis  and  design  procedures.  Implicit  in  their  formulation  is  a  set 
of  assumptions.   These  are  first  listed  and  subsequently  discussed. 
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1.  The  D  +  T  +   He4(3.52  MeV)  +  n( 14.06  MeV)  reaction  is  the 
dominant  fusion  reaction. 

2.  The  ion  density  is  maintained  at  50%  -  50%  D-T. 

3.  The  plasma  is  homogeneous. 

4.  All  particles  are  at  the  same  Maxwell-Boltzman  temperature  kT. 

5.  Particles  above  kT  instantaneously  slow  down  to  kT. 

6.  Beam  fusion  events  are  neglected. 

7.  Synchrotron,  recombination  and  line  radiation  are  negligible. 

The  validity  of  the  first  assumption  is  quickly  established 
through  a  comparison  of  the  D-T  reactivity  with  those  of  all  other 
fusion  reactions.   It  is  seen  to  be  two  orders  of  magnitude  larger,  up 
to  a  temperature  of  70  keV,  and  no  less  than  one  and  one-half  orders  of 
magnitude  larger  from  70  to  100  keV,  than  the  reactivities  of  competing 
reactions.   Having  thus  established  the  dominance  of  the  D-T  reaction, 
the  fusion  reaction  rate,  given  by  nDn„<DT> ,  is  maximized,  under  the 
constraint  of  constant  total  number  of  particles, 

for  nD  =  n„  =  n/2 ,  _n  denoting  in  this  case  the  total  fuel  ion  density. 
Since  each  D-T  reaction  removes  one  deuterium  and  one  tritium,  the 
ratio  of  the  two  will  be  maintained.   One  additional  consideration,  the 
fractional  buildup  concentration  of  product  alpha  particles,  given  by 
n<DT>/4,  can  be  shown  to  be  less  than  10%  for  total  ion  densities  less 
than  10^/cra  at  20  keV.   The  first  two  assumptions  given  above  are 
then  seen  to  be  valid  and  self-consistent. 

In  order  to  obtain  the  spatial  distribution  of  density  and 
temperature  in  a  plasma,  it  is  necessary  to  solve  the  macroscopic 
equations  obtained  by  taking  moments  of  the  Boltzmann  transport 
equation  with  a  Fokker-Plank  collision  term  and  then  averaging  over 
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magnetic  surfaces.   Workers  in  the  field  who  have  performed  such 
analyses  have  obtained  profiles  which  are  basically  parabolic  in  nature 
[6J,  with  the  exact  shape  factors  depending  strongly  on  the  transport 
mechanisms  and  boundary  conditions  assumed.   Alternately,  profiles 
obtained  by  solving  the  set  of  equations  consisting  of  the  electron  and 
power  balance  equations  with  radial  heat  conduction,  Faraday's  law, 
Ampere's  law,  Ohm's  law,  and  charge  neutrality  also  show  extreme 
sensitivity  to  the  scaling  laws  used  [7].   The  inclusion  of  spatial 
effects,  with  the  attendant  complexities,  would  quickly  dispose  of  one 
of  our  initial  objectives  for  this  model,  that  of  explicit  physical 
interpretation  of  the  ensuing  mathematical  terms.   This  analysis  will 
therefore  assume  plasma  homogeneity  and  contend  with  global, 
volume-averaged  parameters. 

In  the  presence  of  a  multispecies  plasma,  a  proper  treatment  of 
the  slowing  down  behavior  of  a  fast  particle  as  it  undergoes  Coulomb 
scattering  involves  the  solution  of  a  Fokker-Plank-type  equation  [81, 
requiring  lengthy  numerical  computations.   Similar  results  can  be 
obtained  using  closed  expressions  for  the  time  behavior  of  the  mean 
energy  of  these  fast  particles  as  determined  by  Butler  and  Buckingham 
[9]  and  Houlberg  [10].   Using  the  latter's  expressions,  O'Farrill  and 
Campbell  [11]  and  Hoxie  [12]  computed  slowing  down  times  for  product 
alphas  and  injected  deuterium  ions.   Typical  times  required  to  slov; 
down  to  20  keV  in  a  10  keV  background  plasma  of  density  10   /cm  are 
0.4  seconds  for  a  fusion  alpha  and  0.05  seconds  for  a  75  keV  injected 
deuteron.   These  are  quite  small  when  compared  with  plasma  thermal 
motion  periods  of  the  order  of  10's  of  seconds. 
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As  these  fast  particles  slow  down,  they  will  initially  impart 
their  energy  preferentially  to  the  electrons.   As  an  example,  a  3.02 
MeV  proton  in  a  10  keV  plasma  will  initially  lose  energy  to  the 
electrons  at  a  rate  170  times  greater  than  to  the  ions.   Not  until  the 
particle's  velocity  has  reached  approximately  one-fourteenth  the 
electron  thermal  velocity  will  the  rate  of  energy  loss  be  equivalent  to 
both  species  [9].   The  electrons  and  ions  will  therefore  have  different 
kinetic  temperatures,  resulting  in  an  interchange  of  energy  with  an 
equipartition  time  given  by  Spitzer  [13].   For  a  10  keV  deuterium 
plasma  with  a  20  keV  electron  component  this  time  constant  is  on  the 
order  of  3  x  10   seconds,  again  quite  small  when  compared  to  the 
thermal  stability  time  constants. 

Fast  particles  will  not  just  serve  to  impart  energy  to  the 
background  plasma.   Fast  fuel  ions,  whether  produced  in  secondary 
fusion  reactions  or  beam-injected,  will  also  induce  fusion  reactions. 
Beam  reactions  can  provide  a  major  contribution  to  the  fusion  power 
density  due  to  the  resonance  peak  in  the  D-T  cross  section  at  about  100 
keV.   (In  fact,  this  effect  will  be  employed  to  advantage  in  TFTR,  the 
Tokamak  Fusion  Test  Reactor  currently  under  construction.)  The  fraction 
of  beam  particles,  w,  that  will  be  investigated  here  is  small  enough  (w 
<  0.07)  so  that  equilibrium  conditions  will  not  be  appreciably 
different  except  at  the  very  low  temperatures  (below  6  keV)  [14].   In 
addition,  for  injection  energies  below  60-75  keV,  the  beam  fusion 
contribution  to  the  energy  balance  is  also  negligible  [15].   These  two 
factors  combine  to  warrant  the  exclusion  of  the  beam  interaction  from 
this  analysis. 
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The  isotopes  of  hydrogen  are  essentially  100%  ionized  above  50  eV. 
Therefore  the  onlv  sources  of  recombination  and  line  radiation  would 
arise  from  impurity  ions  present  in  the  plasma.   The  greatest  source  of 
impurity  ions  is  the  first  wall  as  a  result  of  its  bombardment  by 
neutrons,  energetic  charged  particles  and  radiation.   Typical 
impurities  may  be  carbon,  silicon,  oxygen,  beryllium,  and  aluminum, 
with  their  expected  concentrations  varying  widely,  depending  on  the 
impurity  removal  system  considered.   If  the  impurity  concentration  can 
be  kept  below  0.1%,  say  with  a  divertor  of  modest  efficiency,  then  the 
power  radiated  away  by,  for  example,  silicon,  will  be  less  than 
one-fifth  of  the  hydrogenic  brerasstrahlung.   Therefore,  recombination 
and  line  radiation  will  be  neglected. 

The  accurate  estimation  of  synchrotron  radiation  is  not  very 
straightforward.   Synchrotron  radiation  is  emitted  in  a  series  of 
harmonics,  with  the  intensity  of  the  radiation  decreasing  with 
increasing  harmonics.   But  the  plasma  is  opaque  to  radiation  of 
frequencies  lower  than  the  plasma  frequency;  therefore  the  bulk  of  the 
synchrotron  radiation  will  be  reabsorbed  by  the  plasma.   In  addition, 
the  entire  synchrotron  radiation  spectrum  is  of  such  frequency  that  it 
will  be  reflected  by  metallic  surfaces,  resulting  in  the  reentry  into 
the  plasma  of  the  radiation  that  was  not  initially  absorbed. 
Notwithstanding,  the  ratio  of  synchrotron  to  bremsstrahlung  radiation 
as  given  by  simple  expressions,  for  95%  reflecting  walls  in  a  machine 
of  aspect  ratio  of  5  and  toroidal  magnetic  field  strength  of  4  Tesla  at 
a  temperature  of  10  keV,  is  approximately  one-tenth.   The  small 
magnitude  and  the  uncertainty  associated  with  synchrotron  radiation 
therefore  warrants  its  exclusion  from  further  consideration  in  the 
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present  model. 

Having  established  the  range  of  validity  of  the  model,  the  next 
step  in  the  analysis  is  to  obtain  the  plasma  transfer  functions  for 
this  two-input  system.   This  can  be  done  by  expanding  the  nonlinear 
equations  in  a  first-order  Taylor  series  about  an  equilibrium  point. 
The  resulting  time-invariant  linear  differential  equations  describing 
the  time  response  of  the  perturbation  variables  yield  the  desired 
transfer  functions.   The  derivation  of  the  system  transfer  functions 
follows . 


2.2.   Open-  and  Closed-Loop  Plasma  Transfer  Functions 

In  the  steady-state,  the  time  derivatives  in  equations  2.1  and  2.2 
are  zero,  resulting  in  the  following  equations: 


n  2 
S0  -  ^-<HT>0  -  -2-  =  0  2.5 

o 


n  "-  3n  kT 

ST   +  0  _2_<nT>  -  bBn/(kT  )l/2  -  _^ — 2.  =  0        2.6 
o  s     xx  4     o    B  o    o' 


T,  =  Tcno*(kTQ)ra  2.7 


n„     E„ 
o      o 


The  subscript  zero  denotes  each  variable's  value  at  steady-state. 


2.8 
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Perforraing  a  perturbation  expansion  of  equations  2.1  -  2.4  about  the 
equilibrium  conditions  dictated  by  equations  2.5  -  2.8  results  in  the 
desired  linearized  equations.   Therefore,  into  equations  2.1  -  2.4, 
substitute  the  following: 

n  =  nQ  +  n' 
kT  =  kTQ  +  kT' 
S  =  SQ  +  S' 

Ts  =  TsQ  +  V 

<DT>  =  <DT>Q  +  ^^]0kT'  =  <DT>Q  +  <DT>gkT' 


^  =  V1  +  ^^  +  ^"]°kT'  3  TE°_1  +  W'    +  TE3kTkT' 


where  the  following  definitions  and  expansions  apply: 


CkT0  +  kT')l/2-   CkT0)l/2(l    +  £f)1/2   =    (kTo)1/2  +   2(kTT)l/2 


t  =         E      1      =   -oT   n   ~i~l(kT   ~)~m  =   -         l 

TE8n   -       3n     Jo  *Tcno  llcV  ^T7 


3x   _1 

TF,,         =  E        ]n    =    -»t    n    ~£(kT     )-m"1     =    -  m 

LdKl  ^b-T     JO  C    O        v        O'  uy    T, 
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1  3TE 


3t„  l  ,  3lr 


n3n     3n  Jo   R   3n  Jo     RnQTE   ~   nQTn 


■1 

1    0LE 


rn3kT  =  -wT\o  =  r  ~wT\°  RkT0xE  kT0x 


The  resulting  equations  are 


lUl  =  S'  +  f- -  nn<DT>nV  +  f- - —  -  _^-XDT>,  -  KclkT'   2.9 

n_  o    t„ 


T        -    3kT  s 

^dT       [  3^  J  3n   ls 


(D     +  6kT  9  D         , 

01     6        °    )<DT>0  "  §VkTo>1/2  +  (*  "   iXl^Mjn' 

o  E„ 


.  r A  +  6kToln  <DT>      Vo     m(R-  1)  -R 
+  [( 12 Jn°<DT>3  "  6(kT0)"2  +  R^ 


Ts0  "  3kTo       Sn 

-  Ks(-^- )  -  KT   °  ]kT'  2.10 

S     3nn         Xs3nAJ 


For  coavenience  rewrite  equations  2.9  and  2.10  as 


—  =   S'    +  an'    +  bkT'  2.11 

dt 


rllcT' 

S'    +  2oT   '    +  en'    +  dkT'  2.12 


dt 


'1°     T  s2xs 


with 
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a  = -  n0<DT>0  2.12a 


2 

mn  n   *• 

b  =  - —  -  _<DT>,   -   Kc  2.12b 

RkT0tE  2  3  b 


c   =    [-«     6        °)<DT>Q  -  |bB(kT0)l/2  +   U  ~    l)[JL-±)Vr0     2.12c 


.   _    fQg  +  6kTo,n    _,  _Vo  m(R-    1)    -   R 

d  "  t — 12 — Jno<DT>3  ~  6(kT  yen  + sr 

0  tto 


-  glKg  -  g2KT  2.12d 


2.12e 


2.12f 


Taking  the  Laplace  transform  of  equations  2.11  and  2.12  gives 

sN(s)  =  S(s)  +  aN(s)  +  bT(s) 

sT(s)  =  g^Cs)  +  g2Ts(s)  +  cN(s)  +  dT(s), 

where  now  N,  T,  S,  and  T  are  functions  of  the  Laplace  transform 
variable  s;  they  represent  the  Laplace  transforms  of  the  perturbation 


Si 

= 

\ 

"  3kTQ 

3no 

§2 

= 

So 
3nn 

-21- 

variables  n' ,  kT' ,  S' ,  and  T  '  respectively,  and  are  not  to  be  confused 

with  any  of  the  function-of-time  variables  of  equations  2.1  to  2.4. 

The  desired  transfer  functions  can  now  be  determined.   The 

open-loop  transfer  functions  (i.e.,  K  =  KT  =  0  in  coefficients  b  and 

s 

&         * 

d  above  and  denoted  as  b  and  d  )  are  given  by 

Gl(.)  B   «5l  =  §l(S'a  +  C/gl) 

S<s>    s2  +  (-a  -  d*)s  +  (ad*  -  cb*) 


G  (s)  =  **L  -  §2(S  '  a) 

2      Vs)   s2  +  (-a  -  d*)s  +  (ad*  -  cb*) 


2.14 


and  are  presented  in  functional  block  diagram  form  in  Figure  2-1.   As  a 
consequence  of  the  principle  of  superposition  of  linear  systems,  the 
temperature  perturbation  can  be  computed  independently  for  each  of  the 
two  system  inputs  and  the  results  added  to  yield  the  total  plasma 
temperature  deviation. 

The  closed-loop  plasma  system  is  given  in  block  diagram  form  in 
Figure  2-2.   Manipulations  of  the  block  diagram  equations  can  be  used 
to  yield  an  expression  for  T(s)  given  by  (dropping  the  function-of-s 
notation) , 

SGj  +  TSG2 

T  = 


1  +  G,H,  +  G9H 


Once  again,  transfer  functions  can  be  obtained  for  each  of  the  inputs 
individually,  and  the  principle  of  superposition  can  then  be  used  to 
obtain  the  total  plasma  temperature  variation.   For  simple  proportional 
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S(s) 


Ts(s) 


T(s) 


Figure  2-1.   Open-loop  plasma  system  block  diagram 
representation. 
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Figure  2-2.   Closed-loop  plasma  system  block,  diagram 
representation. 
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feedback,  i.e.,  Hi  =  K  and  H2  =  K.p  ,  the  closed-loop  transfer 
functions  are  given  by 


T         G,  g, (s  -  a  +  c/g, ) 

j  =±  =  i =  i i 2.15 

S   1  +  GjHj  +  G2H2   s2  +  (_a  _  d)s  +  (ad  _  cb) 


T  Go  gn(s  -  a) 

j  =  _L  =  1 =  £ 2.16 

Tg    1  +  G1H1  +  G2H2    s2  +  (_a  _  d)s  +  (ad  _  cb) 


The  effect  of  proportional-integral-derivative  feedback  control  is 
evaluated  and  discussed  in  Chapter  5. 

The  working  model  and  transfer  functions  described  in  this  chapter 
provide  the  necessary  basis  for  the  analysis  of  a  fusioning  plasma  as  a 
control  system  and  the  design  of  a  suitable  controller.   The  open-loop 
analysis  is  presented  in  Chapter  3,  followed  by  the  classical 
closed-loop  analysis  in  Chapter  4. 


CHAPTER  3 
CLASSICAL  OPEN-LOOP  ANALYSIS  OF  PLASMA  MODEL 


The  transfer  functions  obtained  in  the  previous  chapter  are  of  the 
form  of  a  classical  second-order  system  with  a  finite  zero.   This  lends 
them  handily  to  analysis  via  classical  methods,  e.g.,  root  locus,  Bode 
and  Nyquist  plots,  etc.,  and  to  characterization  via  classical 
specifications,  e.g.,  peak  overshoot,  settling  time,  bandwidth, 
resonance  peak  and  frequency,  etc.   In  fact,  analytical  expressions  can 
be  found  for  the  classical  specifications  for  the  model  at  hand  and 
these  are  derived  in  Appendix  C.   Due  to  the  premium  placed  on  computer 
time,  and  the  desire  to  explore  the  broadest  parameter  range  possible, 
the  analytical  expressions  were  used  in  the  calculations  performed  in 
the  analyses  of  this  chapter.   In  the  first  section  of  this  chapter, 
the  coefficients  of  the  open-loop  plasma  transfer  functions  are 
associated  with  physical  processes  present  in  a  fusioning  plasma.   The 
system  performance  is  then  presented.   These  sections  are  followed  by  a 
summary  of  the  open-loop  results. 


3.1.   Open-Loop  Plasma  Transfer  Function  Characteristics 


The  open-loop  plasma  transfer  functions  correspond  to  equations 

2.13  and  2.14  with  K  =  KT  =  0.   These  are  now  repeated. 
s     xs 
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G,(s)e4^t  = 3-1 

S(s)   s2  +  (-a  -  d*)s  +  (ad*  -  cb  ) 


T(s)  _  82(s  "  a)^ 

(ad"  -  cb") 


G2(s)  ~  YU)  ~    2  ,  , 7*; — 7~*        '* 

LsKS)        s  +  (-a  -  d  )s  + 


3.2 


The  supporting  equations  2.12a-f  with  K  =  KT  =0  reduce  to  the 

s    ls 

following: 


a  =  A_L_i  -  no<DT>o  3.2a 


b  = —  -  ^-<DT>,,  3.2b 

RkToxE     2     9 


Q  +  6kT  ,  „    , 

:  a  6     )<DT>0  -  |vkV1/2  +  (*  -  1}^lr^-)kTo  3-2< 


.     Q  +  6kT  bnn        ,D    n    „ 

v    12    J  °    3   6(kT  )l/Z       Rt 


3. 2d 


Ts  -  3kTo 

81--^ 3.2e 

3n„ 


S 

«2  =  T2-  3'2f 

'    3n 


These  are  in  a  form  that  can  be  compared  to  that  of  the  classical 
second-order  system  with  a  finite  zero  given  by 
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„   t.)    -         S'(S  +  Zl>  -        Bl(S  +  H)   -  3.3 

s2  +  25ii>  s  +  u)n2        (s  +   a)2  +  wd2 


where 

g.  =  gain  of  i—  transfer  function,   i  =  1,2 

-z.  =  location  of  zero  of  i—  transfer  function 

-o  -   3^a   =  location  of  poles 

£  =  damping  ratio 

<jo  =  undamped  natural  frequency 

a   =  £wn  =  damping  factor  (0  <  £  <D,  and 


9 
Uj  ■  (U-/1  -  5   ■  damped  natural  frequency. 

A  comparison  of  equations  3.1  and  3.2  with  3.3  can  be  used  to 
equate 


§1  = 


Ts  /3  -  kTQ 

_2 3.4 


So 
52  =  3n" 


3.5 


:1  =  "a  +  G/gl  =  2(T   -  3kT  )no<DT>o  "   T   -  3kT 


(£  -  l)(Tg   -  3RkTQ) 
R*E<Ts   "  3kV 


3.6 
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z2  =  -a  =  no<DT>0  -  |-±  3.7 


^2  =  ad*  _  cb*  =   bkBTn°2«DT>0  -  2kT0<DT>3) 
o 


+  {[2(1  -  m)<DT>n  +  (H   -  l)kT  <DT>J[— iL_- 
o  o     d  L   12Rt 


(Q„  +  6RkTQ)n0. 


E  kTo 


U^o,  +  Vo(4m-  I*"        I   -  1         3  8 
6R^EokTo      6RTEo(kT0)^2     RTe  2 


a*  i  Q  +  6kT  b„n 


2u>„     2Util  °    °   l    12    ;  °  *   3   6(kTQ)l/2 


Rt,,        J 


Equations  3.4  to  3.9  provide  the  means  by  which  physical  associations 
can  be  attached  to  the  control  system  terminology. 

3.1.1.   Open-Loop  Gains 

The  expression  for  g^  is  discussed  first.   A  Maxwellian  plasma  of 

3 
total  density  nt  at  temperature  kT  has  on  the  average  — ntkT  of  kinetic 

energy.   Assuming  the  ions  and  electrons  to  be  at  the  same  temperature 

and  the  plasma  to  be  electrically  neutral  and  hydrogenic,  the  average 

kinetic  energy  in  terms  of  the  ion  density  is  3nkT.   A  unidirectional 

beam  of  ions  at  temperature  kT,  however,  has  average  kinetic  energy  of 

only  nkT.   Therefore  a  beam  of  particles  with  injection  energy  T  =  3kT 
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injected  Into  a  plasma  at  temperature  kT  would  have  no  effect  on  the 
plasma  temperature,  regardless  of  its  intensity.   It  follows  then,  that 
the  beam  could  either  cool  or  heat  the  plasma,  depending  on  whether  its 
injection  energy  is  below  or  above  3kT,  kT  being  the  temperature  of  the 
plasma.   This  is  the  statement  made  by  equation  3.4  for  gj ,  that  the 
gain  of  the  feedrate  transfer  function  is  proportional  to  the 
difference  in  energy  between  the  injected  beam  and  the  background 
plasma. 

The  expression  for  g2»  equation  3.5,  also  permits  straightforward 
physical  interpretation.   It  is  seen  to  be  the  ratio  of  injected  ion 
density  to  plasma  ion  density.   Therefore,  the  magnitude  of  the  effect 
of  injection  energy  will  be  proportional  to  this  ratio  of  densities. 

Figure  3-1  plots  the  gains  gi  and  g2  as  a  function  of  plasma 
temperature  for  0  and  75  keV  injection  energy,  normalized  by  the  value 
of  the  corresponding  input  function.   As  was  shown  in  Section  2.2,  the 
closed-loop  transfer  functions  are  proportional  to  the  gains  g^  and  g2« 
Since  the  total  response  of  the  plasma  is  the  sum  of  the  feedrate  and 
injection  energy  transfer  functions,  the  plots  in  Figure  3-1  show  that 
considerable  flexibility  is  available  for  affecting  the  closed-loop 
plasma  performance. 

3.1.2.   Zeroes 

In  this  section,  the  physical  basis  of  each  zero  found  in  the 
plasma  transfer  functions  is  presented.   The  effect  of  a  zero  in  the 
response  characteristics  of  a  second-order  system  can  be  determined 
from  the  derivations  provided  in  Appendix  C.  The  effect  is  noted  where 
appropriate  within  the  discussions  of  the  classical  specifications. 
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With  T  =  0,  the  expression  for  z,  reduces  to 


y0<DT>0 ,    2hB%      t- 1  ,  in 

21  =  "    6kT     +3(kT)l^"lT~ 


The  three  terras  on  the  right-hand-side  (RHS)  of  this  equation  can  be 

associated  with  the  thermonuclear  fusion  power,  bremsstrahlung 

radiation  power,  and  energy  confinement  (loss  rate),  respectively. 

Figure  3-2  is  a  plot  of  these  terms  and  their  sum,  as  a  function  of 

plasma  temperature.   It  can  be  seen  that  following  an  initial 

singularity  corresponding  to  the  ignition  temperature,  the  zero  is  a 

smooth,  slowly-varying  function  of  temperature.   Throughout  the  range 

of  temperature  under  investigation,  the  bremsstrahlung  contribution  to 

the  zero  is  negligible,  while  the  fusion  power  terra  is  approximately 

twice  that  of  the  confinement  term  and  of  opposite  sign.   The  terms 

combine  to  produce  a  left-hand-plane  (LHP)  zero  of  magnitude  on  the 

order  of  0.1  -  0.2  sec   .   With  T  =  75  kev,  equation  3.6  for  z, 

s         '   ^  l 

applies  and  its  terras  are  plotted  in  Figure  3-3.   The  presence  of  the 

T   -  3kT  terra  in  the  denominator  accounts  for  the  appearance  of  the 

singular  point  in  Figure  3-3.   The  zero  approaches  positive  infinity  as 

kT  approaches  T  /3  and  reappears  at  negative  infinity  as  kT  becomes 
o 

slightly  greater  than  T  /3.   The  reason  why  the  associated  transfer 
o 

function  does  not  similarly  "blow  up"  is  that  its  gain  contains  the 

same  T   -  3kT   terra  in  the  numerator  and  therefore  approaches  zero  at 
so      ° 

the  same  rate.   As  mentioned  previously,  the  implications  of  a  zero  on 
the  system  response  will  receive  further  discussion  under  the 
appropriate  system  specifications. 
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With  or  without  injection  energy,  the  expression  for  %2   is  given 
by  equation  3.7.   This  zero  is  comprised  of  a  fuel  burnup  terra  and  a 
particle  confinement  term.   For  zero  injection  energy,  confinement 
effects  dominate  by  an  order  of  magnitude  throughout  nearly  the  entire 
range  of  temperature.   The  result,  shown  in  Figure  3-4,  is  a  fairly 
slowly-varying  RHP  zero  of  magnitude  0.1  -  0.2  sec"  .   With  7  5  keV 
injection  energy,  due  to  lower  ni  requirements,  both  the  burnup  and 
confinement  terms  are  lower  in  magnitude  than  for  zero  injection 
energy.   However,  as  the  temperature  is  lowered,,  the  reduction  in  x^ 
is  greater  than  the  reduction  in  n,  which  coupled  to  the  l/tg 
dependence  of  the  confinement  term  allows  confinement  effects  to 
dominate  the  magnitude  of  the  zero  only  below  30  keV.   This  zero  is  of 
magnitude  0.05  -  0.1  sec"   in  the  RHP  and  is  shown  in  Figure  3-5. 

3.1.3.   Undamped  Natural  Frequency,  Damping  Ratio,  and  the  Poles  of  the 
Open-Loop  Plasma  Transfer  Function 

Examination  of  the  expression  for  the  square  of  the  undamped 
natural  frequency,  equation  3.8,  reveals  four  terms  corresponding  to 
the  mathematical  coupling  of  the  competing  physical  processes.   The 
terms  are  expressed  below  for  easy  reference. 


bB  V 


3.11 


— Xf-Pr(<DT>„  -  2kT  <DT>a) 

6(kT  )1/2     °      °     3 
o 

[2(1   .  n)<DT>o  +   U   .    l)kT0<BT>3)[-°2R  -- 1    °]-  -2gL__^     3.12 


(Q     +   6RkT   )n„        Q   n  <DT> 
x^  o      o  -i      xa  o  ' 
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The  four  terms  can  be  identified  with  the  coupling  of  fusion  processes 
and  bremsstrahlung,  expression  3.11;  fusion  processes  and  confinement, 
expression  3.12;  bremsstrahlung  and  confinement,  expression  3.13;  and 
lastly,  of  particle  and  energy  confinement,  expression  3.14.   These  are 
displayed  as  a  function  of  plasma  temperature  in  Figures  3-6  through 
3-9,  covering  the  cases  of  0  and  75  keV  injection  energy,  with  constant 
and  Bohm  confinement.   The  reason  for  plotting  the  square  of  the 
natural  frequency  becomes  apparent  upon  inspection  of  the  figures,  as 
the  sum  of  the  terms  is  found  to  be  negative  below  the  critical 
temperature.   This  is  sufficient  to  indicate  thermal  instability.   Note 
however,  that  for  stability,  both  the  natural  frequency  and  the  damping 
ratio  must  be  greater  than  zero. 

One  of  the  characteristics  of  the  four  figures  which  comes  to  the 
foreground  immediately  is  the  small  magnitude  of  the  contribution  of 
bremsstrahlung  radiation  to  the  natural  frequency  oscillation  of  the 

thermal  balance.   This  can  be  attributed  to  the  assumption  of  a 

2 
hydrogenic  plasma,  i.e.,  free  of  impurities.   Noting  the  Z  dependence 

of  bremsstrahlung  emission  suggests  the  use  of  the  impurity 

concentration  as  a  potential  control  mechanism;  indeed,  this  has 

recently  been  proposed  for  the  latest  Tokamak  reactor  design  [16]. 
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The  dominating  terms  influencing  the  natural  frequency  are  the 
fusion  processes-particle  and  energy  confinement  cross-product, 
expression  3.12,  and  the  particle  confinement-energy  confinement 
cross-product,  expression  3.14.   The  latter  is  seen  to  be  fairly 
independent  of  the  injection  energies  and  confinement  schemes 
considered.   Its  variation  with  temperature  can  be  attributed  directly 
to  the  required  ni  condition  for  steady  state  operation  and  is  fairly 
constant.   To  a  great  extent  then,  the  former  term  (i.e.,  the  terra 
given  by  3.12)  contains  the  determining  characteristics  responsible  for 
the  variations  of  the  natural  frequency  with  plasma  temperature. 

Rearrangement  and  closer  inspection  of  expression  3.12  for  the 
fusion  processes-confinement  term  reveals  it  to  be  composed  of  a 
particle  burnup-energy  confinement  term  and  an  energy 
production-particle  confinement  term.   These  are  given  respectively  by 

0  \2(\   -   mKDT>  +  (I  -    l)kT  <DT>J  3.15 


2tE 


-[2(1  -  m)<DT>0  +  (I  -    l)kT0<DT>3 


and 


Q  n„ 
1   "      -[-m<DT>0  +   (i  -    l)kT0<DT>g]  3.16 


12RxF  kT 
Go  ° 

The  fractions  outside  the  brackets  are  very  slowly  varying  functions  of 
temperature;  therefore  the  main  contributions  to  the  variations  in 
natural  frequency  can  be  further  localized  to  the  terms  inside  the 
brackets,  i.e.,  the  functional  dependence  of  the  confinement  scheme, 
and  the  reactivity  function  and  its  derivative.   For  both  constant  and 
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Bohm  confinement,  in  the  low  temperature  region,  a  change  in 
temperature  due  to  a  perturbation  induces  a  change  in  the  reactivity  of 
a  magnitude  large  enough  to  perpetuate  the  perturbation  in  its  original 
direction.   As  the  temperature  increases,  the  incremental  increase  in 
the  reactivity  decreases  in  magnitude  such  that  at  the  critical 
temperature  the  accompanying  changes  in  the  confinement  characteristics 
have  compensated  for  the  change  in  reactivity,  initiating  the  natural 
oscillations  of  the  thermal  balance. 

The  stabilizing  influence  of  the  inverse  temperature  dependence  of 
Bohm  confinement  can  be  seen  by  comparing  Figure  3-6  with  3-7,  and 
Figure  3-8  with  3-9.   Note  that  the  figures  all  have  the  same  abscissa 
scale,  but  different  stability  ranges.   The  plasma  can  begin  stable 
natural  oscillations  at  a  temperature  that  is  11  keV  lower  with  Bohm 
confinement  than  with  constant  confinement  with  either  zero  or  75  keV 
injection  energy.   The  effect  of  75  keV  injection  energy  is  a  4  keV 
lowering  of  the  temperature  required  for  stable  natural  oscillations, 
mainly  due  to  the  different  equilibrium  conditions  required  for  steady 
state. 

The  actual  frequency  of  oscillation  will  depend  on  the  presence  of 
damping  modes  in  the  plasma  energy  balance.   Equation  3.9  is  the 
analytical  expression  for  the  damping  ratio  in  terms  of  its 
constituents.   Figures  3-10  through  3-17  depict  the  relative  magnitudes 
of  these  terms  as  they  combine  to  produce  the  system  damping. 
Throughout  the  range  of  temperatures  and  conditions  under 
consideration,  the  strong  damping  effect  of  diffusion  is  delicately 
balanced  by  an  opposing  force  from  the  fusion  reactivity. 
Bremsstrahlung  provides  a  small  but  positive  contribution  to  the 
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overall  damping  of  the  system.   For  constant  confinement  and  no  energy 
injection,  the  system  is  overdamped  except  for  a  narrow  band  in  the  50 
to  60  keV  range.   The  addition  of  injection  energy  energy  lowers  the 
critical  temperature  and  the  required  ni  condition.   This  shifts  the 
damping  balance  in  favor  of  the  fusion  reactivity  and  the  system  is  now 
underdaraped  from  the  critical  temperature  up  to  approximately  40  keV. 
For  Bohm  confinement,  the  situation  is  reversed.   For  zero  injection, 
Bohm  diffusion  is  not  as  strongly  damping  as  constant  diffusion,  and 
the  system  is  underdaraped  throughout  the  temperature  range.   With  75 
keV  injection,  diffusion  damping  increases  and  the  system  becomes 
overdamped  above  20  keV. 

Figures  3-18  and  3-19  show  the  undamped  and  damped  natural 
frequencies  for  0  and  75  keV  injection  and  constant  and  Bohm 
confinement.   For  temperatures  approximately  5  keV  higher  than  the 
critical  temperature,  the  undamped  natural  frequency  is  found  to  be  0.1 
to  0.2  sec   for  constant  confinement  and  0.2  to  0.3  sec"   for  Bohm 
confinement.   The  damped  natural  frequency  is  found  in  the  underdamped 
temperature  bands  previously  identified  and  is,  of  course,  always 
smaller  than  the  undamped  frequency.   These  figures  summarize  the 
discussion  in  the  preceding  paragraphs  and  serve  as  stepping  stones  to 
the  presentation  of  the  system  poles. 

The  poles  of  the  open  loop  transfer  function  are  shown  in  root 
contour  form  in  Figures  3-20  and  3-21  with  temperature  as  the  varying 
parameter  along  the  curves.   Both  figures  have  the  same  scale  and 
range.   The  real  part  of  the  poles  is  plotted  as  a  function  of 
temperature  in  Figures  3-22  and  3-23.   Complex  poles  are  found  for  the 
temperature  ranges  corresponding  to  the  presence  of  the  damped  natural 
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frequencies  of  Figures  3-18  and  3-19.   These  in  turn  correspond  to  the 
damping  ratio  being  in  the  range  of  0  to  1.   With  the  exception  of  the 
zero  injection  energy,  constant  confinement  case,  the  system  poles  are 
complex  in  the  vicinity  of  the  critical  temperature.   The  complex  part 
of  the  pole  is  given  by  the  damped  natural  frequency.   The  real  part  of 
the  pole,  that  portion  which  governs  the  stability  of  the  system,  is 
given  by  the  damping  factor,  a   =  ?%>  which  in  terras  of  system 
parameters,  can  be  expressed  as 


Q  +  6kT  bRn„ 


1-Jt+m  +  R-mR 


3.17 


This  factor  must  be  positive,  i.e.,  the  real  part  of  the  poles  must  be 
in  the  left-hand  s-plane,  for  the  system  to  be  stable. 

Equation  3.17  is  equation  3.9  multiplied  by  m^.      Therefore  the 
relative  magnitudes  of  the  terms  in  equation  3.17  can  be  extracted  from 
Figures  3-10  through  3-17.   By  inspecting  the  figures,  it  can  be  seen 
that  the  culprit  in  causing  the  thermal  instability  at  low  temperatures 
is  the  magnitude  of  the  slope  of  the  fusion  reactivity  in  this 
temperature  region.   Stability  is  achieved  by  the  natural  increase  in 
magnitude  of  the  reactivity  and  the  decrease  in  magnitude  of  its  slope 
with  increasing  temperature.   Bremsstrahlung  and  diffusion  losses 
contribute  to  the  stability  of  the  system  throughout  the  temperature 
range.   The  contribution  of  diffusion  to  stability  diminishes  with 
increasing  temperature  as  a  result  of  the  longer  confinement  times 
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required  for  equilibrium. 

The  parameters  of  this  section  completely  define  the 
characteristics  of  the  open-loop  plasma  behavior.   The  resulting  system 
response  is  presented  in  the  next  section. 


3.2.   Open-Loop  System  Performance 

The  design  of  a  closed-loop  control  system  begins  with  a  given  set 
of  performance  specifications  in  the  time  and/or  frequency  domain. 
This  demands  that  these  specifications  be  compatible  and  realistic.   In 
order  to  obtain  a  "ball  park"  feel  for  what  the  system  response  should 
be  for  a  fusioning  plasma  operating  in  a  closed-loop, 
feedback-stabilized  thermal  equilibrium,  the  open-loop  system 
specifications  are  presented  in  this  section  and  examined  in  light  of 
the  parameters  discussed  in  the  previous  section. 

The  time  domain  specifications  for  the  plasma  thermal  balance  are 
presented  in  Figures  3-24  through  3-27.   These  figures  show  relative 
temperature  change,  peak  overshoot,  M   time-to-peak,  t  ,  and  settling 
time,  t  .   Results  are  presented  for  0  and  75  keV,  with  constant  and 
Bohm  confinement,  for  step  inputs  in  feedrate  and  injection  energy. 
The  frequency  domain  specifications  of  bandwidth,  resonance  peak  and 
resonance  frequency  are  presented  in  Figures  3-28  through  3-30. 

The  transfer  functions  of  the  plasma  thermal  balance  involve  input 
functions  that  may  differ  by  many  orders  of  magnitude  over  a  broad 
range  of  equilibrium  conditions.   To  provide  a  normalized  basis  for 
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comparing  the  system  responses  from  such  diverse  conditions,  the  Final 
Value  Theorem  was  used  to  provide  weighting  factors.   The  weighting 
factors  were  obtained  from  the  final  value  of  the  transfer  functions 
assuming  step  inputs  of  a  magnitude  equal  to  the  steady  state  value  of 
the  input  functions.   The  resultant  weighting  factors,  shown  in  Figure 
3-24,  can  then  be  used  in  conjuction  with  the  specifications  of  Figures 
3-25  through  3-30,  which  are  presented  for  unit  step  inputs,  in 
relative  comparisons  of  system  performance. 

In  general,  Figure  3-24  indicates  that  the  transfer  function  with 
respect  to  feedrate  responds  to  a  step  input  with  changes  in  steady 
state  temperature  that  are  3  and  8  times  larger  than  those  of  the 
transfer  function  with  respect  to  injection  energy  for  constant  and 
Bohm  confinement,  respectively.   This  implies  that  the  magnitude  of  the 
plasma  response  is  more  sensitive  to  changes  in  feedrate.   Physically, 
this  can  be  explained  as  follows.   An  increase  in  S  results  in  an 
increase  in  n,  which  in  turn  increases  kT  by  the  square  of  the  increase 
in  n  due  to  the  enhanced  reaction  rate.   An  increase  in  Tg ,  however, 
produces  only  a  proportional  increase  in  kT.   This  behavior  will  be 
more  fully  exploited  in  Chapter  4. 

In  addition  to  the  permanent  temperature  change  caused  by  a  step 
change  in  the  input,  the  plasma  temperature  will  experience  an  initial 
transient  component.   The  magnitude  of  this  initial  transient  is 
characterized  by  its  peak  overshoot  (M  ) .   The  time  behavior  of  this 
component  is  described  by  the  time  required  to  reach  its  peak  (t  ),  and 
the  time  required  for  its  magnitude  to  settle  to  5%  of  its  final  value 
(t  ).   These  three  specifications  are  shown  in  Figures  3-25  through 
3-27. 
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Peak  overshoot  values  above  25  to  30  per  cent  are  seldom  allowable 
in  control  engineering  applications.   Inspection  of  Figure  3-25  quickly 
shows  that  in  most  cases,  plasma  behavior  is  unacceptable  below  the 
temperature  range  of  20  to  30  keV.   The  large  overshoots  were  to  be 
expected,  as  peak  overshoot  is  strongly  dependent  on  the  system 
damping,  and  as  presented  earlier,  damping  is  very  low  below  this 
temperature  range.   An  interesting  point  to  note  is  the  presence  of 
overshoot  in  overdamped  regions  of  Figure  3-25.   In  a  simple  classical 
second  order  system,  overshoot  does  not  occur  for  overdamped 
conditions.   However,  due  to  the  presence  of  a  finite  zero  in  the 
plasma  transfer  functions,  overshoot  can  occur  in  an  overdamped  case  if 
the  condition  is  met  that  the  zero  not  be  located  between  the  poles,  as 
derived  in  Appendix  C.  This  is  not  a  sufficient  condition,  only  a 
necessary  one.   This  condition  is  satisfied  for  the  75  keV  injection 
cases  for  both  constant  and  Bohm  confinement. 

The  overshoot  problem  in  both  the  underdamped  and  overdamped  cases 
is  less  severe  for  constant  confinement  than  for  Bohm  confinement  due 
to  the  smaller  pole-zero  separation  in  the  former  case.   This  can  be 
seen  with  the  aid  of  equations  C.ll  and  C.12  and  Figures  3-5,  3-22,  and 
3-23.   As  with  the  final  steady-state  temperature  changes  (Figure 
3-24),  the  maximum  overshoot  is  also  more  sensitive  to  changes  in 
feedrate  than  in  injection  energy. 

Along  with  the  overshoot  problem,  Figures  3-26  and  3-27  show  the 
plasma  to  have  very  poor  speed  of  response.   In  the  vicinity  of  the 
critical  temperature  (when  the  damping  is  generally  small),  and  at 
other  temperatures  where  the  damping  is  high,  settling  times  can 
rapidly  become  greater  than  100  seconds.   Present  conceptual 
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designs[17]  project  pseudo-steady-state  burn  times  on  the  order  of  500 

seconds  for  Tokamak  reactors.   If  this  were  indeed  the  case, 

approximately  one-fifth  of  the  power  production  phase  of  operation 

would  be  required  to  arrive  at  the  desired  steady-state  operating 

conditions.   Control  is  clearly  suggested. 

In  the  frequency  domain,  bandwidth,  resonance  peak  and  resonance 

frequency  are  shown  in  Figures  3-28  through  3-30.   The  system  bandwidth 

and  resonance  frequency  are  seen  to  increase  very  rapidly  in  the 

vicinity  of  the  critical  temperature.   Somewhat  above  the  critical 

temperature,  the  bandwidth  remains  essentially  constant  with 

temperature  .   The  resonance  frequency  bears  a  functional  relationship 

to  the  damping  ratio  and  decreases  as  the  system  becomes  overdamped. 

Thus,  no  resonance  frequency  is  found  for  G,  with  To  =0  and  constant 

o 

confinement.   Damping  in  this  case  is  greater  than  or  nearly  equal  to  1 
at  all  temperatures.   The  variation  of  the  resonance  frequency  of  the 
transfer  functions  for  the  other  cases  can  similarly  be  related  to  the 
variation  in  the  damping  ratio.   Lastly,  the  resonance  peak  is  seen  to 
be  unacceptable  in  most  temperature  regions  of  interest,  exceeding  the 
1.1  to  1.5  design  range  of  practical  control  systems.   These 
characteristics  correlate  very  well  with  the  speed  of  response  noted 
above  from  time  domain  considerations. 

This  section  has  presented  the  classical  specifications  for  the 
open-loop  plasma  system.   In  most  cases,  these  have  been  shown  to  be 
unacceptable.   This  can  be  attributed  to  the  extreme  sensitivity  of  the 
plasma  in  the  vicinity  of  the  critical  temperature  resulting  from  the 
rapid  change  of  plasma  parameters  in  this  temperature  region.   In  the 
next  chapter,  the  system  is  closed  using  proportional  feedback  and  the 
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effect  on  the  system  characteristics  and  specifications  is  noted. 
Derivative  and  integral  feedback  considerations  are  presented  in 
Chapter  5. 


3.3.   Summary 

In  this  section  the  selected  plasma  model  was  cast  in  the  form  of 
a  classical  second  order  system  with  a  finite  zero,  and  analytical 
expressions  for  the  system  characteristics  and  specifications  were 
derived  in  terms  of  the  plasma  conditions.   Expressions  were  derived 
for  the  system  characteristics  of  gain,  the  finite  zero,  the  damped  and 
undamped  natural  frequencies,  the  damping  ratio,  and  the  damping 
factor.   Expressions  were  derived  for  the  impulse  and  step  response,  as 
well  as  the  system  specifications  of  time-to-peak,  peak  overshoot, 
settling  time,  bandwidth,  resonance  frequency,  resonance  peak,  and 
phase  margin.   These  were  evaluated  for  the  cases  of  zero  and  75  keV 
injection  energy,  assuming  both  constant  and  Bohm  confinement  laws. 
The  injection  energy  case  was  examined  because  this  is  the  currently 
preferred  form  of  plasma  heating.   Bohm  confinement  was  evaluated 
because  the  transport  expected  in  the  reactor  regime  is  expected  to 
exhibit  a  Bohm-type  of  dependence. 

Physical  associations  were  drawn  between  the  processes  present  in 
a  fusioning  plasma  and  the  system  characteristics  and  transfer  function 
coefficients.   The  gain  of  the  plasma  system  with  respect  to  the  fuel 
injection  rate  was  found  to  be  proportional  to  the  difference  between 
the  beam  energy  and  the  temperature  of  the  background  plasma.   With 
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respect  to  the  injection  energy  input,  the  gain  of  the  plasma  system 
was  determined  to  be  proportional  to  the  ratio  of  the  fuel  ion  density 
in  the  beam  to  the  fuel  ion  density  on  the  background  plasma. 

The  dominant  processes  influencing  the  natural  frequency  of  the 
plasma  where  found  to  be  the  functional  dependence  of  the  confinement 
scheme,  and  the  magnitude  of  the  reactivity  function  and  its 
derivative.   At  low  temperatures  such  as  between  the  ignition  and 
critical  temperatures,  a  perturbation  which  induces  a  change  in  the 
energy  balance  perpetuates  itself  in  the  original  direction  and  induces 
unstable  behavior.   Above  the  critical  temperature,  the  change  in  the 
fusion  energy  production  rate  caused  by  a  perturbation  is  accompanied 
by  an  opposing  change  in  the  diffusion  rate.   These  opposing  effects 
result  in  the  natural  frequency  of  oscillation  of  the  plasma. 

The  actual  frequency  of  oscillation  depends  on  the  presence  of 
damping  modes  in  the  energy  balance.   Throughout  the  range  of 
temperatures  and  conditions  studied,  diffusion  exerts  a  strong  damping 
influence  on  the  thermal  balance.   The  fusion  reactivity  on  the  other 
hand  has  a  strong  disruptive  effect  on  the  thermal  balance,  up  to  high 
temperatures  in  the  range  of  50  to  60  keV.   In  this  temperature  range, 
the  slope  of  the  reactivity  function  begins  to  decrease  in  magnitude, 
ultimately  reversing  in  sign.   Above  the  temperature  at  which  it 
reverses  in  sign,  the  reactivity  also  adds  to  the  system  damping. 

An  additional  zero  in  a  second  order  system  generally  decreases 
the  time  to  peak  and  increases  the  overshoot.   The  zero  in  the  present 
transfer  function  arises  from  the  need  for  two  independent  differential 
equations  to  aptly  describe  the  plasma  thermal  balance.   The  location 
of  the  zero  is  determined  from  the  relative  magnitudes  of  terms 
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associated  with  the  rates  of  energy  production  from  fusion  and  energy 
loss  via  diffusion.   The  contribution  of  bremsstrahlung  radiation 
emission  to  the  location  of  the  transfer  function  zero  and  the  other 
system  characteristics  of  natural  frequency,  damping  ratio,  and  the 
damping  factor  was  determined  to  be  minor.   For  the  specific  case  of 
complex  poles,  the  damping  factor,  i.e.,  the  real  part  of  the  poles,  is 
the  product  of  the  damping  ratio  and  the  natural  frequency,  and  the 
above  discussion  of  these  system  characteristics  applies.   When  both  of 
these  are  greater  than  zero,  the  system  is  stable. 

The  present  analysis  provided  an  assessment  of  the  relative 
magnitude  of  the  contributions  of  each  of  the  processes  present  in  the 
plasma  to  the  overall  thermal  response.   The  open-loop  system  analysis 
provided  an  initial  evaluation  of  the  system  characteristics  and 
specifications.   Thus,  the  natural  frequency  of  the  thermal  balance  was 
determined  to  be  in  the  range  of  0.05  to  0.3  sec   for  temperatures 
greater  than  a  few  keV  above  the  critical  temperature.   This 
corresponds  to  a  time  scale  for  the  thermal  response  on  the  order  of  3 
to  20  seconds.   This  time  scale  verifies  the  assumption  made  in  Chapter 
2  of  the  exclusion  from  the  analysis  model  of  effects  which  occur  on 
more  rapid  time  scales,  such  as  the  slowing  down  of  alpha  particles  and 
beam-injected  fuel  ions. 

The  damping  ratio  of  the  open-loop  thermal  balance  was  found  to  be 
highly  sensitive  to  the  plasma  parameters.   For  the  case  of  no 
injection  energy  and  a  constant  confinement  law,  the  damping  ratio  is 
much  greater  than  1  in  the  vicinity  of  the  critical  temperature, 
decreasing  to  approximately  1  at  the  higher  temperatures.   With  no 
injection  energy  and  the  Bohm  confinement  law,  the  system  is 
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underdamped  in  the  vicinity  of  the  critical  temperature,  increases  with 
increasing  temperature,  and  levels  off  at  ~  0.5  to  0.6.   With  injection 
energy  and  either  constant  or  Bohm  confinement,  the  system  is 
underdamped  in  the  vicinity  of  the  critical  temperature,  increases  with 
temperature,  and  then  levels  off  in  the  overdamped  range  of  1.0  to 
1.25.   These  considerations  play  an  important  role  in  the  selection  of 
the  plasma  and  reactor  operating  conditions  and  the  design  of  the 
control  system. 

The  manifestation  of  the  system  characteristics  on  the  performance 
of  the  open-loop  thermal  balance  was  evaluated  via  the  system 
specifications  of  peak  overshoot,  time-to-peak,  settling  time(5%), 
bandwidth,  resonance  peak,  and  resonance  frequency.   This  work  is  the 
first  evaluation  of  the  response  characteristics  of  the  plasma  thermal 
balance.   Such  an  evaluation  is  a  required  first  step  in  determining 
the  type  of  control  best  suited  for  this  application.   The  results  of 
this  phase  of  this  study  are  summarized  in  Table  3-1  and  described 
below. 

•   Zero  injection  energy,  constant  confinement.   This  system  is 
overdamped  in  the  entire  temperature  range  of  applicability  of  the 
model.   It  therefore  does  not  exhibit  overshoot.   This  condition  is 
accompanied  by  very  poor  speed  of  response,  with  settling  times  greater 
than  100  seconds  at  temperatures  below  30  keV.   In  agreement  with  the 
time  domain  specifications,  the  system  bandwidth  is  0.01  sec"   slightly 
above  the  critical  temperature  and  remains  very  low,  less  than  0.2 
sec-  at  temperatures  above  30  keV,  with  no  resonance  peak. 
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•  Zero  injection  energy,  Bohm  confinement.   The  strong  effect  of  the 
confinement  law  on  the  sytem  response  is  evident  from  a  comparison  of 
the  specifications  of  this  case  with  those  of  the  constant  confinement 
case.   With  Bohm  confinement,  the  system  exhibits  unacceptably  high 
overshoot,  greater  than  25%  below  40  keV,  and  a  resonance  peak  greater 
than  2  throughout  the  applicable  temperature  range.   Its  speed  of 
response  is  improved  over  the  constant  confinement  case  in  that  the 
time-to-peak  is  15  to  25  seconds,  and  the  settling  time  is  on  the  order 
of  30  to  60  seconds  above  25  keV.   The  bandwidth  is  also  increased  to 
~0.4  -  0.7  sec-  at  temperatures  above  15  keV.   A  resonance  peak  of 
magnitude  greater  than  2  is  found  in  the  0.2  to  0.3  sec-   frequency 
range  below  30  keV. 

•  7  5  keV  injection  energy,  constant  confinement.   The  introduction  of 
injection  energy  into  the  plasma  balance  serves  to  initiate  plasma 
oscillations  below  35  keV.   Peak  overshoot  is  unacceptably  high  up  to 
25  -  30  kev,  and  decreases  rapidly  between  25  and  35  kev.   There  is  no 
significant  overshoot  above  35  kev.   In  general,  the  speed  of  response 
is  improved  over  the  no  injection  case,  though  it  is  still  not  as  good 
as  with  Bohm  confinement  and  zero  injection  energy  at  temperatures 
below  25  keV.   Time-to-peak  and  the  settling  time  are  both  greater  than 
1  minute  at  temperatures  below  20  -  25  keV,  decreasing  to  30  seconds  at 
30  keV.   The  bandwidth  is  not  affected  greatly  by  the  75  keV  injection 
energy.   Between  20  and  25  keV,  a  resonance  peak  of  magnitude  greater 
than  2  appears  in  the  frequency  range  of  0.01  to  0.1  sec 

•  7  5  keV  injection  en^ergy^_Bolp_conMnement .   This  combination  of 
input  and  plasma  parameters  demonstrates  the  most  favorable  open-loop 
system  response.   Its  stable  operating  temperature  range  is  the  widest, 
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starting  as  low  as  10  keV.   Acceptable  overshoot  is  found  starting  at 
15  keV.   In  addition,  reasonable  settling  times  of  less  than  30  seconds 
are  also  attained  at  temperatures  greater  than  15  keV,  and  the  system 
bandwidth  is  in  the  1  sec-   range.   The  magnitude  of  the  resonance  peak 
is  reduced  significantly  with  75  keV  injection  energy.   The  resonance 
frequency  is  in  the  range  of  0.1  -  0.3  sec-   above  10  keV. 

This  brief  summary  demonstrates  the  broad  range  of  system  response 
that  results  from  open-loop  operation.   Though  the  system  response  is 
improved  via  the  lse  of  injection  energy  and  as  a  result  of  the  Bohm 
confinement  law,  the  response  can  be  further  improved  through  the  use 
of  closed-loop  feedback  control. 


CHAPTER  4 
CLASSICAL  CLOSED-LOOP  ANALYSIS  OF  PLASMA  MODEL 


Analysis  of  the  thermal  balance  in  the  previous  chapter  showed 
that  open-loop  operation  of  a  fusioning  plasma  requires  undesirably 
high  temperatures  for  stability.   It  was  also  demonstrated  that  even  at 
the  elevated  temperatures,  the  open-loop  system  possesses  poor 
performance  characteristics,  including  excessive  overshoot  and  slow 
speed  of  response.   The  open-loop  plasma  system  can  be  transformed  into 
a  closed-loop  system  through  the  addition  of  feedback,  i.e.,  by 
modifying  the  input  variables  in  accordance  with  the  behavior  of  the 
system  variables.   The  closed-loop  plasma  system  may  then  be  able  to 
sustain  stable  performance  at  lower  temperatures  which  would  be 
open-loop  unstable.   The  closed-loop  system  may  also  be  able  to  improve 
the  system  response  at  temperatures  which  would  be  open-loop  stable, 
but  would  possess  poor  response  characteristics. 

The  effect  of  closing  the  loop  on  the  stability  and  system 
response  of  the  plasma  thermal  balance  is  presented  in  this  chapter. 
In  the  first  section,  the  loop  is  closed  using  proportional  feedback. 
The  closed-loop  transfer  functions  are  developed,  and  the  system 
characteristics  are  derived  and  evaluated.   The  system  performance  is 
then  presented  and  the  improvements  over  open-loop  operation  are  noted. 
The  chapter  concludes  with  a  summary  of  the  closed-loop  results. 
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4.1.  Closed-Loop  Plasma  Transfer  Function  Characteristics 

The  closed-loop  transfer  function  equations  are  given  by  equations 

2.13  and  2. 14  with  K  and  KT  not  equal  to  zero.   These  are  repeated 
s       s 

below  for  convenience. 


T(s) 


gj(s  -  a  +  c/gj) 


G  (*)    =    l^s>    =  L ± 4.1 

Gl(-S')  "  qf<^    ?  .  ,     ..   .  ,    , 


S(s) 


r   +  (-a  -  d)s  +  (ad  -  cb) 


G  (s)  =       =  - 4«2 

Ts(s)   s2  +  (-a  -  d)s  +  (ad  -  cb) 


The  supporting  equations  2.12a-f  remain  as  before. 


a  =  - -  n0<DT>Q  4.2a 


b  =  —  -  — <DT>-  -  K  4.2b 

RkTQTE     2     3    s 


Q  +  6kT  9  r  -  l 


e  =  { j<DT>0  -  fbB(kT0)l/2  +  U  -    D(^— ^)kTn   4.2c 


.    rQa+6kT0l n  <DT>       bBno    +  m(R-  1)  -R 
d  =  I 12 )no<DT>3  "  6(kT  )l/2  +  R^ 


;iKs  -  g2K  4. 2d 


Ts   -  3kTQ 

8l--^ 4'2e 

1       3n 
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_^_  4.2f 

3n„ 


Starting  with  the  above  equations,  expressions  for  the  system 
characteristics  can  be  derived  for  the  closed-loop  transfer  function  in 
the  same  manner  as  previously  used  for  the  open-loop  case.   After  some 
manipulation,  the  expressions  are  given  by 


T_  /3  -  lcTf 
so 


So 
g2  =  3n~ 


(Qa+2Ts0}  2bBn0(kT0)l/2 

21  =  _a  +  C/gl  =  2(TS  -  3kTQ)n°<DT>°  "   Ts   -  3kTQ 


(£  -  1)(T_   -  3RkT) 
so 


4.3 


4.4 


Rtp(T   -  3kTj 


4.5 


lE^s 


z2  =  -a  =  n0<DT>0  -  f-±  4.6 


2     2 

U)„   =  U), 


open  loop 


Qa  +  2Tso 
+  [  ( T-^<DT>c 


(£  -  l)(3RkTQ  -  Ts  ) 

-h«VU2  + 5E^ "K 


S0<DT>0    S  (£  -  1) 

+  _2 °-_2 KT  4.7 

1    3        3noV  J  Ts 


-80- 

§1      §2 
5  _  ^open  loop  +  2^5  +  J^TS  4*8 


Comparing  the  closed-loop  equations  4.3  through  4.8  with  the 
open-loop  equations  3.4  through  3.9,  it  can  be  seen  that  only  the 
natural  frequency  and  damping  ratio,  and  therefore  the  poles  of  the 
system  have  changed  by  the  addition  of  terms  proportional  to  the 
feedback  coefficients  K  and  KT  .   These  additional  terms  can  be  used 
to  effect  limited  control  over  the  system  performance.   The  effect  of 
the  new  terms  on  the  plasma  thermal  balance  characteristics  is 
discussed  below.   The  closed-loop  system  performance  is  presented  in 
Section  4.2. 

The  characteristic  of  the  plasma  thermal  balance  which  is  of 
primary  importance  is  its  stability.   As  previously  discussed,  the 
stability  of  the  present  system  is  determined  by  the  position  of  the 
closed-loop  poles.   Using  equations  4.1  and  4.2,  and  the  supporting 
equations  4.2a  through  4. 2d,  the  closed-loop  poles  can  be  determined  as 
a  function  of  the  feedback  fractions.   Alternately,  using  the  same 
equations,  the  feedback  fractions  required  to  guarantee  stable 
performance  can  be  determined  for  any  given  plasma  temperature.   This 
latter  approach  was  adopted  for  presentation  and  the  results  of  the 
analysis  are  presented  in  Figures  4-1  through  4-4.   These  figures 
graphically  portray  the  temperature  range  which  can  be  stabilized  with 
the  use  of  proportional  feedrate  and/or  injection  energy  feedback. 
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Figure  4-1  shows  that  with  zero  injection  energy,  the  minimum 
stable  temperature  can  be  lowered  to  18  keV  from  the  minimum  open-loop 
temperature  of  25  keV  using  only  2  to  3%  feedrate  feedback  fraction, 
assuming  constant  confinement.   With  the  Bohm  confinement  law  and  no 
injection  energy,  Figure  4-2  shows  that  the  minimum  stable  temperature 
can  be  lowered  to  11  keV  from  the  open-loop  value  of  14  keV  using  a  -5% 
feedrate  feedback  fraction.   Figures  4-1  and  4-2  further  show  that 
there  is  a  broad  range  of  feedback  fraction  that  can  be  used,  even  at 
temperatures  that  are  open-loop  stable,  to  modify  and  improve  the 
system  performance.   The  selection  of  the  operating  feedback  fractions 
is  discussed  below  and  in  subsequent  sections. 

With  75  keV  injection  energy,  both  the  feedrate  and  injection 
energy  feedback  fractions  can  be  selected  so  as  to  establish  stability. 
The  stabilized  regimes  for  constant  and  Bohm  confinement  are  plotted  in 
Figures  4-3  and  4-4,  parametrically  with  the  injection  energy  feedback 
fraction.   The  figures  show  that  for  either  confinement  scheme,  with 
proper  selection  of  the  feedback  fractions,  stable  operation  is 
possible  to  as  low  as  5  keV  (the  ideal  ignition  temperature  at  which 
the  fusion  yield  begins  to  exceed  the  minimum  ideal  radiative  losses). 
This  temperature  can  be  compared  to  the  open-loop  critical  temperatures 
of  21  keV  and  10  keV  for  constant  and  Bohm  confinement,  respectively. 
Again,  the  figures  show  that  a  broad  selection  of  feedback  fractions  is 
available  with  which  to  adjust  the  system  performance.   The  factor  of 
10  change  in  the  ordinate  scale  is  to  be  noted  in  comparing  Figures  4-1 
and  4-2  with  4-3  and  4-4. 
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An  interesting  characteristic  of  the  stability  range  observed  in 
Figure  4-1  through  4-4  is  that  at  temperatures  above  the  critical 
temperature  with  no  injection  energy  and  above  temperatures  greater 
than  approximately  25  keV  with  75  keV  injection  energy,  there  is  an 
upper  limit  to  the  feedrate  feedback  fraction  above  which  unstable 
behavior  is  again  encountered.   This  limit  can  be  explained  with  the 
use  of  equation  4.8  for  the  system  damping,  and  the  plots  of  the  gains 
gi  and  go  presented  in  Figure  3-1  of  the  previous  section.   For 
stability,  both  u>  and  E,   must  be  greater  than  zero.   However,  for  zero 
injection  energy,  Figure  3-1  shows  g^  to  be  less  than  zero  for  all 
temperatures.   With  75  keV  injection  energy,  Figure  3-1  shows  that  the 
gain  gi  becomes  negative  above  25  kev  and  will  therefore  result  in  a 
decrease  in  the  system  damping  with  increasing  feedrate  feedback 
fraction.   For  a  sufficiently  high  fraction,  i.e.,  the  limit  prescribed 
in  Figures  4-1  through  4-4,  the  damping  ratio  falls  below  zero  and  thus 
destabilizes  the  system.   Additional  discussion  of  the  effect  of 
feedback  on  the  damping  ratio  is  presented  below. 

Given  the  range  of  feedback  that  can  stabilize  the  thermal 
balance,  the  selection  of  the  feedback  fraction  is  guided  by  the  choice 
of  system  performance  specifications.   As  developed  in  Appendix  C, 
these  specifications  are  determined  by  the  system  characteristics  of 
gain,  the  zero,  natural  frequency,  and  the  damping  ratio.   Of  these, 
only  the  natural  frequency  and  the  damping  ratio  are  affected  by 
closed-loop  feedback  as  given  by  equations  4.7  and  4.8.   Thus,  the 
natural  frequency  and  damping  ratio  (and  the  system  specifications)  are 
functions  of  both  the  temperature  of  the  plasma  and  the  feedback 
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representative  temperatures  in  the  following  operating  scenarios: 

1.  Zero  injection  energy,  constant  confinement  -   20,  25,  and  30  keV. 

2.  Zero  injection  energy,  Bohm  confinement  -   15,  20,  and  25  keV. 

3.  75  keV  injection  energy,  constant  confinement  -   10,  20,  and  30  keV. 

4.  75  keV  injection  energy,  Bohm  confinement  -   10,  20,  and  30  keV. 

The  square  of  the  natural  frequency  of  the  closed-loop  thermal 
balance  is  given  by  equation  4.7.   Examination  of  the  new  terms  in 
equation  4.7  reveals  that  the  three  physical  processes  of  fusion  gain, 
bremsstrahlung,  and  confinement  losses,  all  contribute  to  the  effect  of 
feedrate  feedback  on  the  undamped  natural  frequency,  while  only  fusion 
terms  and  confinement  losses  contribute  to  the  effect  of  injection 
energy  feedback.   The  three  terms  contributing  to  the  effect  of 
feedrate  feedback  are  plotted  in  Figure  4-5,  while  the  two  terms 
contributing  to  the  effect  of  injection  energy  feedback  are  plotted  in 
Figure  4-6.   Figure  4-5  applies  for  zero  and  75  keV  injection  energy, 
for  both  constant  and  Bohm  confinement.   Figure  4-6  applies  for  both 
confinement  types,  but  only  for  the  case  of  75  keV  injection  energy. 
The  five  terms  in  the  two  figures  are  normalized  by  the  value  of  the 
input  control  variable.   The  net  effect  of  feedback  on  the  square  of 
the  natural  frequency  is  the  product  of  the  multiplier  shown  in  the 
figures  and  the  feedback  coefficient.   A.s  in  the  open-loop  case, 
bremsstrahlung  is  a  weak  contributor  to  the  natural  frequency,  while 
the  burnup  and  confinement  terms  are  the  primary  contributors. 
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Inspection  of  Figures  4-5  and  4-6  shows  that  the  net  effect  of 
feedback  on  the  natural  frequency  is  weaker  at  the  lower  temperatures 
of  interest  (<10  keV)  but  increases  to  appreciable  levels  at  higher 
temperatures  (>20  keV) .   The  figures  also  show  that  feedrate  feedback 
has  a  stronger  effect  on  the  natural  frequency  than  injection  energy 
feedback,  by  factors  of  2  to  5,  in  the  temperature  range  of  5  to  30 
keV.   These  considerations  are  quantitatively  demonstrated  in  the  plots 
of  the  closed-loop  natural  frequency  versus  feedback  fraction  shown  in 
Figures  4-7  through  4-13. 

For  the  cases  of  no  injection  energy,  Figure  4-7  shows  that  the 
natural  frequency  is  in  the  range  of  0.02  to  0.2  sec   .   In  each  case, 
the  curve  drawn  in  the  figure  spans  the  stable  feedrate  feedback  range. 
The  general  trend  is  seen  to  be  a  rapid  increase  of  the  natural 
frequency  in  the  vicinity  of  the  minimum  feedrate  feedback  fractions 
necessary  for  stability,  followed  by  a  slowing  down  to  a  more  moderate 
rate  of  increase  at  the  higher  feedrate  feedback  fractions.   At 
temperatures  for  which  stability  can  only  be  provided  by  a  very  narrow 
band  of  feedrate  feedback  fractions,  e.g.,  20  keV  and  constant 
confinement,  the  natural  frequency  shows  high  sensitivity  to  the 
feedrate  feedback  fraction  and  does  not  reach  a  plateau  region.  As  the 
band  for  stabilizing  feedrate  feedback  increases,  e.g.,  20  to  25  keV 
and  Bohm  confinement,  the  sensitivity  decreases  at  the  higher  feedrate 
feedback  fractions. 

With  75  keV  injection  energy,  the  natural  frequency  is  further 
increased  to  the  0.1  to  1.0  sec"  range  (Figures  4-8  through  4-13). 
This  is  due  to  the  combination  of  both  increased  permissible  feedrate 
feedback  resulting  from  the  lower  system  gain  and  the  additional 
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contribution  of  the  injection  energy  feedback  fraction.   The  trend  of 
the  variation  of  the  natural  frequency  with  feedrate  feedback  fraction 
is  otherwise  very  similar  to  that  discussed  above  for  zero  injection 
energy.   Inspection  of  Figures  4-8  to  4-13  shows  that  the  effect  of 
injection  energy  on  the  natural  frequency  is  small.   This  can  be  seen 
from  equation  4.7  and  the  plots  of  Figures  4-5  and  4-6.   The  figures 
show  that  terms  multiplying  the  feedrate  feedback  fractions  are  5  to  10 
times  greater  than  the  terms  multiplying  the  injection  energy  feedback. 
Equation  4.7  shows  this  advantage  to  be  due  mainly  to  the  burnup  term 
multiplying  Kg. 

The  damping  ratio  of  the  closed-loop  thermal  balance  is  given  by 
equation  4.8.   The  new  terms  in  the  expression  for  damping  essentially 
consist  of  the  system  gains  divided  by  twice  the  natural  frequency. 
Since  the  gains  are  not  changed  by  the  addition  of  feedback,  the  causes 
of  the  change  in  damping  are  related  to  those  causes  effecting  the 
changes  in  the  natural  frequency,  as  discussed  above.   The  effect  of 
feedback  on  the  system  damping  is  displayed  in  Figures  4-14  and  4-15. 
These  are  plots  of  the  two  new  terms  in  equation  4.8,  normalized  by  the 
value  of  the  input  variable. 

The  two  new  terms  in  equation  4.8  contain  the  natural  frequency  as 
a  factor  in  the  denominator.   Since  the  application  of  feedback  alters 
the  natural  frequency,  Figures  4-14  and  4-15  are  plotted  using  the 
open-loop  natural  frequency.   In  order  to  isolate  the  contribution  of 
feedback  to  the  system  damping,  the  data  in  these  figures  must  be 
combined  with  the  open-loop  natural  frequency  curves  of  Figures  3-6 
through  3-9.   Thus,  the  singular  points  in  Figures  4-14  and  4-15  can  be 
associated  with  the  crossover  temperatures  of  Figures  3-6  through  3-9 
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at  which  the  natural  frequency  is  identically  zero.   At  temperatures 
away  from  the  singular  point,  since  the  natural  frequency  is  generally 
slowly  varying  with  temperature,  the  damping  multiplier  follows  the 
slope  of  the  gain  function  (shown  in  Figure  3-1). 

The  net  effect  of  feedback,  on  the  damping  ratio  is  summarized  in 
Figures  4-16  through  4-22.   With  zero  injection  energy,  Figure  4-16 
shows  that  at  all  temperatures,  the  system  is  overdamped  at  low 
feedrate  feedback  fraction.   The  damping  ratio  decreases  with 
increasing  feedrate  feedback  fraction  until  unstable  behavior  results 
at  sufficiently  high  feedrate  feedback.   With  75  kev  injection  energy, 
Figures  4-17  through  4-22  show  that  damping  always  increases  with 
increasing  injection  energy  feedback.   Using  equation  4.8,  this  can  be 
attributed  to  the  gain  gn   being  greater  than  zero  for  all  temperatures, 
and  the  small  effect  of  injection  energy  feedback  on  the  natural 
frequency.   Physically,  with  increasing  injection  energy  feedback 
fraction,  the  system  energy  balance  is  more  strongly  dominated  by  the 
input  control  function  and  thus  less  sensitive  to  plasma  parameter 
perturbations.   At  temperatures  below  which  the  system  gain  g^  becomes 
negative,  the  system  is  initially  underdamped,  and  the  damping  ratio 
monotonically  increases  with  feedrate  feedback  for  constant  confinement 
(see  e.g.,  Figure  4-17).   At  intermediate  temperatures  and  low 
injection  energy  feedback  fraction,  the  system  damping  with  constant 
confinement  is  similarly  behaved  (see  e.g.,  Figure  4-18).   As  the 
injection  energy  feedback  is  increased,  stability  is  achieved  at  lower 

feedrate  fractions  and  the  KT  term  in  equation  4.8  initially 

s 

dominates.   In  these  constant  confinement  cases,  the  system  is 
initially  overdamped,  passes  through  a  minimum  in  the  system  damping, 
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and  subsequently  increases  with  increasing  feedrate  feedback  fraction. 
The  behavior  of  the  system  with  constant  confinement  at  intermediate 
temperatures  is  also  observed  with  Bohm  confinement  at  low  and 
intermediate  temperatures  because  stability  is  again  achieved  at  lower 
temperatures  and  feedrate  feedback  fractions.   At  high  temperatures 
(above  25  keV),  the  system  is  overdamped  at  low  feedrate  feedback 
fractions  with  both  constant  and  Bohm  confinement,  and  the  damping 
ratio  monotonically  decreases  with  increasing  feedrate  feedback  until 
the  onset  of  system  instability  (see  Figures  4-19  and  4-22). 

The  present  formulation  of  the  problem  of  thermal  stability 
permits  the  investigation  of  the  full  range  of  applicability  of 
feedback.   Past  works  have  been  able  to  assure  thermally  stable 
operation  only  by  observing  the  position  of  the  closed-loop  poles  for 
very  specific  plasma  conditions.   The  general  applicability  of  the 
present  method  in  determining  the  pole  positions  and  the  necessary 
feedback  fractions  is  displayed  in  Figures  4-23  through  4-29.   These 
figures  present  the  variation  in  the  pole  positions  with  feedback  for 
the  cases  presented  previously  in  this  section.   The  methodology  can  be 
used  to  establish  the  stability  range  of  any  other  injection  energy 
and/or  confinement  scheme.   In  past  works,  no  statements  have  been  nor 
could  have  been  made  as  to  the  preferability  of  one  set  of  feedback 
conditions  over  another.   By  correlating  the  feedback  fractions  with 
the  position  of  the  closed  loop  poles  and  the  system  characteristics  of 
natural  frequency  and  the  damping  ratio,  control  can  conceptually  be 
exercised  over  the  temporal  response  of  the  system  in  a  systematic  and 
predictable  fashion. 


-110- 


o 

CZ 

c 

n 

1 1 

o 

E-< 

J-l 

(  1 

o 

en 

CO 

en 

y-i 

L_ 

O      • 

pj 

10     AJ 

•o    OJ 

o 

(1 

a)  e 

1 

CD 

0)     (1) 

<±-i    c 

f    1 

•H 

LJ 

j->    d 

L_ 

5    o 

D 

LJ 
E-h 

C    B 
O   J= 
•H     O 

4J  m 

o 

(1 

03 

ct: 

(    ) 

(0    CO 

> 

LJ 

0)    c 

L_ 

r-H      CO 
O     ■!-> 

CD 

d.  en 

C 

D.    O 
O     O 

ro 


t'O-  £"0- 

ssnod  3Hi  jo  iyyd  ny3d 


-Ill- 


-^ 

CZ 
i — i 

c 
o     • 

Hh 

•w  > 

l_j 

o  ^ 

az 

CO 

Vj    O 

L_ 

^!     II 

o 

00 

V 

CO    H 

~c£ 

0)       ~ 

CD 
Q 

0)     J-J 

u-i     C 

0) 

r,  l 

£    E 

LJ 

4-1     0) 

L_ 

M-l 

CO 

ll  1 

o   o 

~o 

az 

i-> 

(0     4J 

en 

•w   c 

CO    J-> 

r     n 

>    to 

LJ 

C 
0J    o 

L_ 

r-t       O 
O 

■*• 

> 

"o 

a  a) 
o  ^ 

I'O 


1*0-  E"0- 

S310d  3HI  JO  I^Ud  1U3d 


S'0- 


-112- 


ro 


ro-  e*o- 

S310d  3H1  JO  lyyd  1U3d 


S"0- 


-113- 


I'O 


l"0-  £"0- 

S310d  3HI   JO  IdUd  1U3d 


-114- 


to 


E-H 

•H 

n 

J_> 

LO 

CE 

o 

CO      • 

n 

rr 

lJ    > 

L_ 

^ 

cj   O 

L_J 

XI 

en 

T3     II 
0) 
OJ   H 

CD 

Q 

M-i   J«i 

LJ 

X!       - 

Li.l 

U    i-> 

L_ 

•H     C 

S 

D 

E-h 

c    <u 
o   c 

rr 

J-)    H-l 

C£ 

to    c 

•H     O 

f   i 

u   o 

tfl 

r.i 

*1 

L_ 

CD     O 

-H   « 
O 
CX     « 

> 

a.  ai 

O   ^i 

o 

Lo 


0*0       E'O-      "r"0-       9'0- 

S310d  3HI  JO  IdUd  1U3d 


B"0- 


0-I- 


-115- 


> 


o-o  e'o-  *"o-  9'0- 

S310d  3HI  JO  IdUd  1U3d 


8"0- 


Q'l- 


-116- 


5  0) 


>  B 

<v   o 


> 

o  M 


O'O 


E'O-      fr"0-      9*0- 

S310d  3HI  JO  IdUd  HU3d 


-117- 
The  next  section  presents  the  system  performance  for  operation 
within  the  stahilized  regime.   The  information  presented  provides  the 
basis  for  choosing  the  feedback  coefficients  so  as  to  achieve  the 
desired  system  response. 


4.2.  Closed-Loop  System  Performance 

The  thermal  behavior  of  a  fusioning  plasma  has  been  shown  to  be 
unstable  below  the  critical  temperature.   Section  3.1  presented  a 
method  of  determining  critical  temperatures  for  a  wide  variety  of 
operating  parameters.   In  the  majority  of  cases,  even  in  the  stable 
regime,  the  characteristics  and  classical  control  specifications  of  the 
open-loop  system  were  shown  to  be  unacceptable  in  Section  3.2.   The 
previous  section  has  shown  how  feedrate  and  injection  energy  feedback 
can  be  used  to  stabilize  the  plasma  burn  below  the  critical 
temperature.   This  was  accomplished  by  controlling  the  position  of  the 
closed-loop  poles.   Stable  system  performance  was  found  for  a  broad 
range  of  the  feedback  fraction.   This  section  presents  the  variation  in 
system  performance  with  feedback  fraction.   The  results  can  be  used  to 
choose  the  feedback  coefficients  so  as  to  produce  the  desired  plasma 
temporal  response  in  the  stabilized  region  and  to  improve  the  system 
response  above  the  critical  temperature.   The  time  domain 
specifications  are  presented  first  and  are  shown  in  Figures  4-30  to 
4-61.   These  are  followed  by  the  frequency  domain  specifications  in 
Figures  4-62  to  4-85. 
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To  systematically  investigate  the  system  specifications  over  the 
range  through  which  feedback  will  stabilize  the  operating  point, 
reference  is  made  to  Figures  4-1  through  4-4.   The  presentation  in  this 
section  proceeds  by  selecting  several  operating  temperatures  of 
interest  for  each  case  of  zero  and  75  keV  injection,  with  constant  and 
Bohm  confinement,  and  then  examining  the  system  response  following  step 
changes  in  the  input  over  the  range  of  feedback  fractions  through  which 
the  system  is  stable.   Thus  a  parametric  analysis  is  presented  for  each 
specification  as  a  function  of  operating  temperature,  injection  energy, 
confinement  scheme,  and  feedrate  and  injection  energy  feedback 
fraction.   The  representation  of  such  a  broad  spectrum  of  results  is 
facilitated  by  the  method  of  analysis  chosen  for  this  study. 

As  in  the  open-loop  analysis,  the  final  value  theorem  was  used  to 
compare  the  step  responses  of  transfer  functions  involving  a  broad 
range  of  equilibrium  conditions,  inputs,  and  feedback  coefficients. 
Results  are  shown  in  Figures  4-30  to  4-37  and  are  presented  in  terms  of 
the  absolute  temperature  change  that  would  be  brought  about  by  a  step 
change  of  a  magnitude  equal  to  the  full  steady  state  value  of  the 
input.   The  figures  show  that  the  trend  in  absolute  temperature  change 
is  similar  for  all  the  conditions  investigated.   That  is,  temperature 
changes  induced  by  changes  in  feedrate  are  larger  than  changes  induced 
by  corresponding  changes  in  injection  energy.   Also,  as  the  feedback 
fraction  is  increased,  the  temperature  change  brought  about  by  a  step 
input  decreases.   This  latter  trend  is  observed  for  both  feedrate  and 
injection  energy  feedback.   In  effect,  this  states  that  the  feedback  is 
serving  its  purpose  and  is  attempting  to  minimize  changes  in  the 
operating  temperature.   Since  a  perturbation  analysis  was  used  to 
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obtain  these  results,  the  figures  are  actually  applicable  for  step 
changes  that  result  in  less  than  a  10%  change  in  the  steady  state 
values.   As  in  the  open-loop  analysis,  the  closed-loop  results  are  in 
good  agreement  with  time  development  runs. 

The  first  time  domain  specification  presented,  in  Figures  4-38 
through  4-45,  is  peak  overshoot.   Feedback  has  a  pronounced  effect  on 
peak  overshoot.   Whereas  peak  overshoot  is  unacceptable  at  most 
temperatures  of  interest  in  open-loop  operation,  closed-loop 
proportional  feedback  now  yields  overshoots  which  can  be  limited  to 
less  than  20%  by  proper  selection  of  the  feedback  coefficients. 

With  zero  injection  energy  and  Bohm  confinement,  feedback 
fractions  of  -4%  at  15  keV,  -3%  at  20  keV,  and  -2%  at  25  keV  are 
sufficient  to  provide  acceptable  overshoot  (see  Figure  4-38).   With 
constant  confinement,  the  problem  is  obtaining  acceptable  conditions 
for  overshoot,  time-to-peak,  and  settling  time  simultaneously.   In  this 
case,  the  overshoot  is  very  sensitive  to  the  feedback  fraction. 
Therefore,  control  which  may  be  desirable  to  improve  the  speed  of 
response  may  result  in  unacceptable  overshoot. 

With  75  kev  injection  energy,  both  feedrate  and  injection  energy 
can  be  used  to  control  the  overshoot.   As  shown  in  Figures  4-39  through 
4-45,  increasing  the  injection  energy  feedback  fraction  can  always  be 
used  to  lower  the  overshoot.   This  follows  from  the  fact  that  the 
system  damping  increases  with  the  injection  energy  feedback,  whereas 
overshoot  decreases  with  increasing  damping  ratio.   As  discussed 
previously,  the  system  damping  increases  with  feedrate  feedback  at  low 
temperatures  (<20  keV)  for  constant  confinement.   At  these 
temperatures,  overshoot  decreases  with  increasing  feedrate  feedback 
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f raction  (see  Figure  4-39).   As  the  temperature  increases  to 
approximately  20  keV,  overshoot  decreases  with  increasing  feedrate 
feedback  at  low  injection  energy  feedback  fractions,  and  increases  at 
high  fractions  (see  Figure  4-40).   At  still  higher  temperatures  (above 
25  keV) ,  the  overshoot  increases  with  increasing  feedrate  feedback  (see 
Figure  4-41).   The  general  trends  with  Bohm  confinement  are  similar 
except  that  they  they  appear  at  approximately  5  keV  lower  temperatures 
than  with  constant  confinement. 

Inspection  of  the  overshoot  present  under  certain  operating 
conditions,  in  particular  those  represented  in  Figures  4-38,  4-41, 
4-44,  and  4-45,  shows  that  the  permissible  feedback  range  for 
acceptable  overshoot  can  be  quite  narrow  and  exhibit  high  sensitivity. 
This  implies  that  although  operation  with  only  proportional  feedback  is 
possible,  only  limited  flexibility  is  available  in  specifying  the 
feedback  necessary  to  minimize  overshoot.   The  reason  for  the  highly 
sensitive  behavior  can  be  seen  as  follows.   Referring  to  equation  C.ll 
in  Appendix  C,  it  can  be  seen  that  peak  overshoot  is  dependent  on  the 
location  of  the  zero,  the  pole-zero  separation,  and  exponentially  on 
the  system  damping.   It  was  shown  in  the  last  section  that  the  location 
of  the  zero  is  not  affected  by  feedback.   It  was  also  shown  in  the  last 
section  that,  in  general,  the  poles  are  slowly-vaying  functions  of  the 
feedback  fractions  (see,  e.g.,  Figures  4-23  to  4-29).   Thus  the  main 
contributor  to  the  variation  in  peak  overshoot  is  its  exponential 
dependence  on  the  system  damping.   The  strong  variation  in  peak 
overshoot  thus  follows  directly  from  the  strong  variation  of  the  system 
damping  with  feedback  fraction  as  shown  previously  in  Figures  4-16  to 
4-22. 
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The  application  of  proportional  feedback  results  in  marked 
improvements  in  the  speed  of  response.   As  evidenced  by  Figures  4-46  to 
4-53,  the  time-to-peak  can  be  reduced  to  the  seconds  and  tens  of 
seconds  range,  from  the  open-loop  values  of  tens  and  hundreds  of 
seconds  which  are  shown  or  can  be  extrapolated  from  Figure  3-26. 
Except  at  low  temperatures,  this  can  be  achieved  with  values  of 
feedback  fractions  below  10  to  20%. 

The  variation  of  the  time-to-peak  with  feedback  fraction  is  rather 
straightforward  and  similar  for  all  the  cases  investigated.   That  is, 
t  decreases  monotonically  with  both  feedrate  and  injection  energy 
feedback  at  all  temperatures,  with  zero  or  75  keV  injection  energy,  for 
constant  or  Bohm  confinement.   With  no  injection  energy  and  constant 
confinement,  the  time-to-peak  is  reduced  from  hundreds  of  seconds  to 
the  tens  of  seconds  range  (see  Figure  4-46),  although  as  mentioned 
previously,  the  reduction  comes  at  the  expense  of  increased  overshoot. 
With  no  injection  energy  and  Bohm  confinement,  the  generally  limited 
decrease  in  t  obtained  with  the  use  of  positive  feedrate  feedback 
fractions  would  not  be  sufficient  to  warrant  the  increased  overshoot. 
Alternately,  if  increased  t  is  permissible,  decreased  overshoot  is 
possible  with  negative  feedrate  feedback  fractions  (see  the  discussion 
on  p. 61).   With  75  keV  injection  energy,  very  modest  values  of  the 
feedback  fraction  are  sufficient  to  reduce  the  time-to-peak  below  10 
seconds  except  at  low  temperatures  (see  Figures  4-47  to  4-53).   Note 
that  the  time-to-peak  with  75  keV  injection  energy  is  dictated  by  the 
response  of  G2,  since  its  overshoot  is  many  times  greater  than  that  of 
Gi. 
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As  in  the  case  of  the  time-to-peak,  the  behavior  of  the  system 
settling  time  is  similarly  improved  for  operation  with  closed-loop 
feedrate  and  injection  energy  feedback.   Thus,  settling  times  obtained 
during  open-loop  operation  of  tens  and  hundreds  of  seconds  are  reduced 
to  seconds  and  tens  of  seconds  as  shown  in  Figures  4-54  through  4-61. 
The  strong  effect  of  feedback  on  the  settling  time  can  be  traced  with 
the  help  of  equations  4.8  and  C.13.   From  equation  C.13,  it  can  be  seen 
that  the  settling  time  is  inversely  proportional  to  the  damping  factor, 
o.   Equation  4.8  can  be  solved  in  terms  of  the  damping  factor  by 
multiplying  by  w  to  yield 


0  -  »n«  =  "n^open  loop  +  ^s  +  ^Tg  4'9 


The  effect  of  the  feedback  fraction  on  a  for  various  injection  energies 
and  plasma  temperatures  can  then  be  explained  with  equation  4.9  and  the 
aid  of  the  plots  of  the  system  gains  given  in  Figure  3-1. 

With  zero  injection  energy,  gi  is  negative  for  all  temperatures. 
Thus  with  constant  confinement,  since  K  must  be  greater  than  zero  for 
stability  above  25  keV  (as  shown  in  Figure  4-1),  the  settling  time 
cannot  be  decreased  below  some  minimum  which  is  in  the  hundreds  of 
seconds.   At  temperatures  above  25  keV,  as  the  range  of  stabilizing 
feedrate  feedback  fraction  increases,  the  settling  time  can  be  reduced 
to  the  tens  of  seconds  range.   With  Bohm  confinement,  the  settling  time 
can  be  lowered  to  20  -  40  seconds  with  a  negative  feedback  fraction, 
while  simultaneously  decreasing  the  peak  overshoot,  but  increasing  the 
time- to- peak. 
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With  75  keV  injection  energy,  gi  becomes  negative  above  25  kev. 
Below  this  temperature,  increased  feedback  fraction  can  always  be  used 
to  decrease  the  settling  time  (see  e.g.,  Figures  4-55,  56,  59,  and  60). 
This  is  a  beneficial  trend,  since  as  discussed  previously,  overshoot  is 
similarly  decreased  with  increased  feedback  fraction,  for  both  constant 
and  Bohm  confinement.   Above  25  keV,  the  negative  sign  of  gi  reverses 
the  above  trend,  and  at  sufficiently  high  feedrate  feedback  fraction, 
the  settling  time  increases  with  increased  feedrate  feedback  (see 
Figures  4-57  and  4-61).   At  all  temperatures,  increased  injection 
energy  feedback  can  always  be  used  to  simultaneously  reduce  the 
settling  time  and  peak  overshoot.   This  can  again  be  attributed  to  the 
gain  g2  being  greater  than  zero  at  all  temperatures.   The  last  terra  in 
equation  4.9  thus  decreases  with  Km  ,  increasing  the  damping  factor  and 
decreasing  the  settling  time. 

In  the  frequency  domain,  the  specifications  of  bandwidth, 
resonance  peak  and  resonance  frequency  are  shown  in  Figures  4-62  to 
4-85.   The  system  bandwidth  can  be  increased  significantly,  from  one  to 
two  orders  of  magnitude  above  the  open-loop  values,  to  values  from  1  to 
10  sec-  and  higher  in  some  cases,  with  the  use  of  proportional 
feedback  (see  Figures  4-62  to  4-69).   This  is  consistent  with  the 
improved  speed  of  response  discussed  previously,  but  also  indicates 
that  the  plasma  will  exhibit  increased  sensitivity  to  noise  and 
fluctuations  in  the  input  parameters.   As  shown  in  Figures  4-78  to 
4-85,  the  resonance  frequency  is  similarly  affected  by  feedback  and 
increases  by  one  to  two  orders  of  magnitude  to  the  0.1  to  1.0  sec 
range.   The  plots  of  resonance  peak  versus  feedback  fraction,  Figures 
4-70  through  4-77,  are  similar  in  nature  to  the  curves  presented 
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earlier  for  peak  overshoot.   This  might  have  been  expected  given  the 
close  relationship  between  overshoot  and  the  resonance  peak  for  a 
classical  second  order  system.   Although  improvements  over  the 
open-loop  values  of  the  resonance  peak  are  potentially  available  in 
certain  feedback  ranges,  the  large  magnitude  of  the  resonance  peak  is 
still  of  major  concern  in  the  event  that  the  system  stimulates 
frequencies  in  the  vicinity  of  the  resonance  frequency.   At  low 
injection  energy  feedback  fractions,  the  resonance  peak  can  be 
maintained  below  1.5  only  within  a  very  narrow  range  of  feedrate 
feedback.   The  range,  however,  can  always  be  increased  by  increasing 
the  injection  energy  feedback  fraction. 

The  next  section  summarizes  the  results  and  observations  of  this 
chapter,  and  offers  conclusions  on  the  effect  of  proportional  feedback. 
Chapter  5  discusses  the  additional  effect  of  additionally  applying 
derivative  and  integral  control. 


4.3.   Summary 

To  investigate  the  improvement  in  stability  and  system  performance 
obtained  by  using  feedback  control,  expressions  for  the  system 
characteristics  were  derived  for  the  closed-loop  transfer  function  in 
the  same  manner  as  for  the  open-loop  case.   The  loop  was  closed  using 
proportional  feedback  because  this  basic  type  of  control  is  very 
effective,  has  wide  applicability,  is  appropriate  for  an  initial 
investigation  such  as  the  present,  and  is  often  used  in  more  advanced 
control  applications  in  combination  with  additional  types  of  control. 
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The  relative  contributions  of  the  fusion  processes  to  the  effect  of 
feedback  on  the  system  characteristics  were  evaluated.   As  in  the 
open-loop  case,  bremsstrahlung  is  a  weak  contributor,  and  burnup  and 
confinement  terms  are  the  primary  contributors.   The  natural  frequency 
and  the  damping  ratio,  and  thus  the  poles  of  the  system  are  affected  by 
terras  proportional  to  the  feedback  coefficients. 

The  poles  of  the  closed-loop  transfer  function  determine  the 
stability  of  the  system.   By  proper  selection  of  the  feedback 
fractions,  and  through  the  combined  effects  of  feedback  on  the  natural 
frequency  and  the  damping  ratio,  the  real  part  of  the  system  poles  can 
be  forced  into  the  negative  real  plane  and  thus  closed-loop  stability 
can  be  achieved.   The  novel  formulation  of  the  problem  as  presented  in 
this  study  permitted  the  feedback  fractions  required  to  achieve 
stabilization  to  be  graphically  depicted  in  a  simple  form,  dependent 
only  on  the  input  parameters,  the  confinement  law,  and  temperature. 
Thus,  given  the  input  parameters  as  determined  from  power  level  or 
plasma  physics  requirements,  and  the  confinement  law,  the  stabilized 
temperature  range  and  the  required  feedback  can  be  quickly  obtained. 

With  zero  injection  energy  and  constant  confinement,  proportional 
feedrate  feedback  control  lowers  the  critical  temperature  to  18  keV 
from  the  previous  value  of  25  keV  without  feedback  control.   With  Bohm 
confinement,  the  critical  temperature  is  lowered  to  11  keV  from  the 
previous  14  keV.   With  75  keV  injection  energy,  the  effect  of 
proportional  feedrate  feedback  control  in  lowering  the  critical 
temperature  is  dependent  on  the  injection  energy  feedback  fraction.   By 
appropriate  choice  of  these  two  parameters,  the  critical  temperature 
can  be  effectively  lowered  to  the  ignition  temperature,  5  keV.   The 


-182- 
selection  of  the  feedback  fractions  is  then  motivated  not  just  from  a 
requirement  to  operate  in  a  stable  mode,  but  also  by  a  desire  to  obtain 
good  system  response. 

Good  system  response  is  obtained  by  proper  selection  of  the 
natural  frequency  and  system  damping.   For  feedrate  feedback,  burnup 
serves  to  increase  the  natural  frequency,  while  diffusion  decreases  it. 
For  injection  energy  feedback,  both  processes  can  be  used  to  increase 
the  natural  frequency.   The  combined  application  of  both  feedrate  and 
injection  energy  feedback  can  be  used  to  increase  the  natural  frequency 
by  factors  of  3  to  5  over  the  open-loop  values  at  most  temperatures  of 
interest. 

The  most  dramatic  display  of  the  effect  of  feedback  is  the  broad 
range  of  damping  ratio  which  is  attainable  with  the  combination  of  both 
feedrate  and  injection  energy  feedback.   The  entire  range  from  zero 
damping  to  overdamping  is  possible  at  nearly  all  temperatures  of 
interest  by  proper  selection  of  the  feedback  coefficient.   In  some 
ranges  of  feedback  fraction,  however,  the  system  damping  shows  great 
sensitivity  to  the  feedback  coefficients.   This  implies  that  small 
variations  in  the  feedback  fractions  will  have  a  pronounced  effect  on 
the  damping  ratio. 

The  broad  range  of  control  available  with  proportional  feedrate 
and  injection  energy  feedback  can  be  used  to  attain  a  broad  range  of 
system  performance.   The  performance  specifications  of  the  plasma 
thermal  balance  with  closed-loop  proportional  feedback  are  summarized 
in  Tables  4-1  and  4-2.   These  are  described  below  in  the  same  format 
used  for  the  open-loop  specifications. 
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•  Zero  injection  energy,  constant  conf ineraent .   Uhereas  the  open-loop 
system  exhibits  no  overshoot  up  to  temperatures  in  the  40  to  50  keV 
range,  the  closed-loop  overshoot  can  be  controlled  to  values  less  than 
0.25  with  feedrate  feedback  fractions  less  than  0.02  at  temperatures  as 
low  as  20  keV.   At  the  low  overshoot  values  and  a  temperature  of  20 
keV,  the  tlrae-to-peak  remains  rather  high  at  greater  than  100  seconds, 
though  the  settling  time  is  improved  to  times  below  1  minute.   At  these 
conditions  the  system  bandwidth  remains  below  0.3  to  0.4  sec   .   At 
high  feedrate  feedback  fractions,  a  resonance  peak  which  can  attain 
values  greater  than  10  is  now  also  present  at  a  frequency  in  the  range 
of  0.01  to  0.2  sec"  .   Its  magnitude  decreases  to  less  than  2  at  low 
feedrate  feedback  fractions. 

•  Zero  injection  energy,  Bohm  confinement.   The  high  overshoot 
present  in  the  open-loop  system  can  be  reduced  to  an  acceptable  level 
by  proper  selection  of  the  feedrate  feedback  fraction.   The  associated 
effect  on  the  speed  of  response  is  that  the  time-to-peak  can  increase 
into  the  one-minute  range  at  15  keV  with  large  negative  feedrate 
feedback  fractions.   The  settling  time  to  within  5%  of  the  final  value 
improves,  decreasing  at  these  conditions  to  under  40  seconds.   The 
bandwidth  can  be  varied  by  two  orders  of  magnitude  with  feedrate 
feedback,  and  a  resonance  peak  of  significant  magnitude  (greater  than 
5)  is  located  in  the  0.01  to  0.1  sec-   freauency  range  for  high 
feedback  fractions. 

•  75  keV  injection  energy,  constant  confinement.   This  set  of  operating 
condition  still  presents  difficulties  in  achieving  acceptable  overshoot 
at  temoeratures  below  20  keV.   However,  improvement  is  possible  with 
sufficiently  high  injection  energy  feedback  fractions  and  low  feedrate 
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feedback  fractions  at  temperatures  above  20  keV.   For  these  latter 
conditions,  a  time-to-peak  of  5  to  10  seconds  is  attained,  with 
settling  times  less  than  ~30  seconds.   The  bandwidth  is  increased  by 
nearly  an  order  of  magnitude,  up  to  ~2  sec   ,  and  the  resonance  peak 
can  be  significantly  reduced  (into  the  0.2  to  0.3  range,  at  frequencies 
in  the  range  of  0.2  to  0.6  sec-  )  by  appropriate  selection  of  the 
feedback  coefficients. 

•  75  keV  injection  energy,  Bohm  confinement.   As  in  the  case  of  the 
open-loop  system,  this  operating  condition  exhibits  the  most  desirable 
performance  specifications.   Even  at  temperatures  as  low  as  10  keV, 
neak  overshoot  can  be  maintained  below  25%,  the  time-to-peak  below  20 
seconds,  and  the  settling  time  less  than  30  seconds.   Bandwidth  is 
increased  to  the  range  of  2  to  10  sec-  ,  and  the  resonance  peak  remains 
approximately  the  same  as  in  the  open-loop  case. 

The  results  presented  in  this  chapter  show  that  a  significant 
improvement  in  system  performance  can  be  attained  using  only 
proportional  feedback.   The  improvement  is  adequate  for  many  conditions 
of  interest,  for  example,  10  to  20  keV  with  75  keV  injection  and  Bohm 
confinement.   However  some  conditions  still  encounter  difficulties  at 
low  temperatures,  such  as  the  case  of  constant  confinement  with  or 
without  injection  energy.   The  next  chapter  continues  to  explore 
methods  of  improving  the  system  performance  by  investigating  the 
proportional- integral-derivative  controller. 


CHAPTER  5 
PROPORTIONAL-INTEGRAL-DERIVATIVE  CONTROL 


The  classical  control  theory  analysis  of  the  preceding  chapters 
has  shown  how  the  unacceptable  temporal  behavior  of  the  thermal  balance 
of  a  fusioning  plasma  can  be  improved.   The  improvement  was 
accomplished  using  only  proportional  feedback  with  the  fuel  feedrate 
and  injection  energy  as  the  controlling  mechanisms  and  the  plasma 
temperature  as  the  control  variable.   The  improvement  in  system 
response  with  proportional  feedback  was  shown  to  be  sufficient  for  many 
operating  conditions  of  interest  but  inadequate  for  others. 

The  selection  of  the  proper  control  system  for  a  given  application 
is  usually  determined  by  a  compromise  between  the  quality  of  the 
control  desired  and  the  admissible  cost  of  the  control  system.   For  a 
capital  intensive  energy  source  such  as  fusion  is  expected  to  be,  the 
value  of  highly  reliable  and  precise  plasma  operation  will  probably 
allow  for  a  large  initial  capital  expenditure  for  the  control  system  so 
as  to  ensure  high  quality  control.   It  is  thus  desirable  to  investigate 
other  control  schemes  that  may  further  improve  system  response 
characteristics,  extend  the  stable  plasma  temperature  regime,  allow 
greater  operational  flexibility,  and/or  provide  more  precise  control  of 
the  thermal  balance.   This  chapter  presents  extensions  of  the  work  of 
the  previous  chapters  to  indicate  how  these  improvements  might  be 
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accomplished  with  the  use  of  the  proportional-integral-derivative  (PID) 
controller. 

The  major  advantage  of  proportional  control  is  its  powerful 
stabilizing  action,  simplicity,  and  wide  range  of  adjustment.   As 
discussed  in  Appendix  B,  its  major  disadvantage  is  that  proportional 
control  is  subject  to  an  inherent  steady-state  error  under  changing 
load  conditions.   This  error  can  be  eliminated  with  integral  control, 
in  which  the  magnitude  of  the  control  action  is  proportional  to  the 
time  integral  of  the  error.   This  type  of  control  can  thus  provide  a 
permanent  change  in  the  input  variables  so  as  to  eventually  eliminate 
any  deviations  of  the  controlled  variable  from  its  desired  or  setpoint 
value.   The  response  with  proportional  control  can  also  be  improved  by 
supplementing  it  with  derivative  control,  in  which  the  magnitude  of  the 
control  action  is  proportional  to  the  derivative  of  the  error  function 
and  opposite  in  sign.   Thus,  since  derivative  control  is  such  so  as  to 
oppose  all  change  in  the  controlled  variable,  this  control  provides  an 
additional  stabilizing  influence  on  the  system.   Derivative  control  is 
also  characterized  by  rapid  response,  since  it  functions  by 
anticipating  system  response.   The  combined  effect  of  these  three 
control  modes  are  obtained  from  the  proportional-integral-derivative 
(PID)  controller. 

Proportional  control  was  modeled  in  the  formulation  of  the  plasma 
particle  and  energy  conservation  equations  of  Chapter  2  by  terms  of  the 
form 


CSp  =  Ks<kTo  "  kT>  5.1 

CT   =  KT  (kT  -  kT)  5.2 
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in  equations  2.1  and  2.2,  where  Cg  and  C™   represent  the  proportional 
control  contributions  from  feedrate  and  injection  energy  feedback, 
respectively,  Kg  and  K^.  are  the  feedback  coefficients  used  previously, 
kTQ  is  the  desired  or  setpoint  operating  temperature,  and  kT  is  the 
actual  operating  temperature.   Similarly,  integral  and  derivative 
control  can  be  implemented  by  introducing  the  following  terms  into 
equations  2.1  and  2.2, 

CSi  '  Ksi  Jt  <kTo  "  kT>dt  5-3 

CT  .  =  K-,  ,  /   (kT  -  kT)dt  5.4 

s      s   t 


d  kT 
~dT 


CSd  =  Ksd  -rr-  5.5 


CV  "  KTsd  ~iT  5'6 


whe  re , 

Cg^  =  integral  control  contribution  using  fuel  feedrate  with 
feedback  coefficient  K  . 

Cq,  .  =  integral  control  contribution  using  injection  energy  with 

feedback  coefficient  K~,  . 
s 

Cgj  =  derivative  control  contribution  using  fuel  feedrate  with 

feedback  coefficient  K  j 

Cj   ^  =  derivative  control  contribution  using  injection  energy  with 

feedback  coefficient  K™  , 
Tsd 

and  the  integrals  are  taken  over  all  time  t. 
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The  derivation  of  the  new  system  transfer  functions  can  now 
proceed  as  in  Section  2.2.   Thus,  a  perturbation  expansion  can  be  used 
to  linearize  equations  2.1  and  2.2  including  the  new  terms.   The 
linearized  equations  can  then  be  Laplace  transformed,  and  the  resultant 
expressions  solved  for  the  system  transfer  functions.   Alternately,  new 
expressions  for  the  feedback,  transfer  functions  H,  and  H2  can  be 
derived  by  taking  the  Laplace  transforms  of  equations  5.1  to  5.6  to 
give 


ll  "  ^s  T  "T-  T  ^sd= 


H,  =  K„  +  —  +  L,s  5.7 


KT_i 
H2  =  Rj  +  — —  +  KT  ds  5.8 

s     s       s 


Following  the  sign  conventions  established  in  Figure  2-2,  these 
expressions  can  be  substituted  into  equations  2.15  and  2.16,  together 
with  the  expressions  for  Gi  and  Go  given  by  equations  2.13  and  2.14  to 
yield  the  closed-loop  transfer  functions  directly. 

Following  either  approach,  the  system  transfer  functions  are  given 
by 


T      gis(s  -  a  +  c/g,) 
Jj  si- 5.9 

S   <s  <*  +  §lKsd  +§2KT  d>  + 


s2[-a  -  d*  +  gxKs  +  g2KT  +  (c  -  ag1)Kgd  -  g2aKT   ]  + 

s  sd 

s[ad*  -  cb*  +  (c  -  agl)Ks  -  g2aKT  +  glKgi  +  g2KT  1]  + 
[(c  -  ag1)Kgi  -  g2aKT  i]} 
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and 


T  g2s(s  -  a) 


J2  =  _  =  — 5.10 

Ts   {s3(l  +  glKgd  +g2KT  d)  + 


s2[-a  -  d*  +  glKs  +  g2KT  +  (c  -  ag^Kgj  -  g2aKT  ]  + 
s[ad  -  cb  +  (c  -  ag^Kg  -  g2aKT  +  g^si  +  g2KT  ±]  -i 
[(c  -  agl)Ksl  -  g2aKT  ±]} 

These  expressions  can  be  used  to  determine  the  system 
characteristics  and  response  performance  in  terms  of  the  six  individual 
feedback  gains  by  following  procedures  similar  to  those  presented  in 
the  previous  chapters.   For  the  present  work,  it  suffices  to  provide 
numerical  examples  of  plasma  performance  that  display  the  general 
trends  of  the  PID  controller;  the  presented  results  would  thereby  serve 
as  guidance  for  future  efforts. 

The  main  features  of  the  PID  controller  that  are  demonstrated 
below  are  elimination  of  offset  error  by  integral  action  and  improved 
system  response  and  stability  from  derivative  action.   These  effects 
are  illustrated  at  the  optimum  plasma  operating  temperature,  15  keV, 
for  the  case  of  75  keV  neutral  beam  injection  and  constant  confinement. 
To  minimize  the  number  of  variable  input  parameters  and  still 
demonstrate  the  power  and  versatility  of  PID  control,  only  the  effect 
of  injection  energy  feedback  will  be  employed,  i.e.,  it  will  be  assumed 
that  Kg  =  Kgi  =  Kgd  =  0.0.   The  effect  of  PID  control  is  then  developed 
by  successively  adding  integral  and  derivative  control  to  the 
proportional  component  of  injection  energy  feedback. 
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To  observe  the  temporal  plasma  behavior,  the  nonlinear 

differential  equations  2.1  and  2.2,  including  the  control  functions  of 

equations  5.1  to  5.6,  were  numerically  solved  in  time.   At  15  keV,  the 

result  of  a  +10%  temperature  perturbation  is  shown  in  Figure  5-1  for 

three  values  of  the  injection  energy  proportional  feedback  coefficient 

only.   The  equilibrium  is  seen  to  be  unstable  at  the  desired 

temperature  of  15  keV  with  K™  =0.0  and  0.1  (the  plasma  finds  a  new 

s 

equilibrium  at  28  keV  or  extinguishes),  and  stable  with  KT  =  0.2. 

s 

Note  that  the  system  response  with  KT  =0.2  compares  favorably  with 

s 

the  predicted  response  obtained  by  studying  the  stability  diagram  of 
Figure  4-3,  and  interpolating  Figures  4-39  and  4-40,  4-47  and  4-48,  and 
4-55  and  4-56  for  the  expected  overshoot,  time-to-peak,  and  settling 
time. 

A  series  of  dynamic  simulations  of  the  plasma  equilibrium  with 
parametric  variation  of  the  feedback  coefficients  ^  .  and  KT  ,  can 
thus  be  used  to  demonstrate  the  main  features  of  the  PID  controller. 
The  results  are  presented  in  Figures  5-2  to  5-7  and  discussed  below. 

The  effect  of  varying  the  integral  control  feedback  coefficient, 

ICj.  . ,  on  the  plasma  temperature  temporal  response  is  shown  in  Figure 

5-2.   For  a  10%  step  change  in  the  operating  temperature,  from  15  to 

16.5  keV,  proportional  feedback  of  magnitude  K™   =0.2  alone  is  unable 

sP 

to  achieve  the  desired  final  temperature.   By  adding  the  effect  of 

integral  control,  the  16.5  keV  final  temperature  can  be  achieved.   At 

low  values  of  KT  .  the  behavior  is  overdamped  and  slow  to  respond,  and 

at  high  feedback  fractions  it  is  underdamped  and  exhibits  overshoot  and 

oscillatory  behavior.   The  value  of  KT  .  =  0.2  would  be  a  reasonable 

s 

compromise  between  acceptable  overshoot  and  adequate  speed  of  response. 
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With  K™  „  =  0.2  and  K™  .  =  0.2,  the  fractional  contribution  of  the 
%P  Tsi 

two  control  modes  to  the  total  effect  of  injection  energy  feedback  are 

displayed  in  Figure  5-3.   The  case  of  KT  •  =  0.0  is  shown  for 

s 

comparison.   There  are  two  important  features  shown  in  the  figure.   The 
first  is  that  at  steady  state,  the  only  active  control  element  is  that 
arising  from  the  integral  action.   Thus,  integral  action  is  totally 
responsible  for  achieving  the  desired  final  temperature.   The  second 
feature  is  that  since  the  area  underneath  the  curves  serves  as  a 
measure  of  the  energy  requirements  of  the  control,  it  can  be  concluded 
that  proportional-integral  control  is  more  energy-efficient  than 
proportional  control  alone. 

The  effect  of  derivative  control  on  the  response  of  the  plasma 
temperature  is  shown  in  Figure  5-4,  again  for  a  10%  step  change  in 

temperature  at  15  keV  with  KT  =  0.2.   The  figure  shows  that  derivative 

s 

control  attains  a  faster  speed  of  response  at  the  expense  of  increased 

overshoot.   Subsequent  to  its  first  peak,  however,  the  response  does 

not  oscillate,  and  attains  its  final  temperature  at  times  up  to 

one-half  of  the  no-derivative-control  case.   It  can  also  be  noted  from 

Figure  5-4  that  derivative  action  does  not  aid  in  eliminating  offset 

error. 

Using  the  values  of  KT   =  KT  .  =  KT  .  =  0.2,  Figure  5-5  shows  the 
s^     s      s 

fractional  contributions  to  the  total  feedback  of  the  individual 
proportional  and  derivative  control  modes.   This  figure  shows  that 
derivative  control  makes  its  contributions  at  the  early  times,  when  the 
temperature  is  changing  most  rapidly.   Its  effect  diminishes  as  the 
temperature  levels  off.   The  figure  also  shows  that  the  increased  speed 
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of  response  is  obtained  at  a  cost  of  an  increase  in  the  injection 

energy  requirement  from  the  neutral  beams. 

Finally,  the  effect  of  the  PID  controller  is  demonstrated  in 

Figures  5-6  and  5-7.   Thus,  it  can  be  seen  that  with  KT  «■  Km  ,    =  K™  . 

Ts    V    Tsd 

=  0.2,  the  plasma  attains  the  desired  final  temperature  of  16.5  keV 
with  a  noticeably  improved  speed  of  response,  and  with  a  negligible 
increase  in  the  overshoot.   The  price  is  an  increased  power  requirement 
during  early  times,  but  a  lower  overall  energy  requirement  is  achieved. 

The  PID  controller  gains  chosen  for  this  example  may  not 
necessarily  be  the  optimum  selections.   Simultaneous  variation  of  the 
feedback  coefficients  either  experimentally  or  analytically  using 
equations  5.9  and  5.10  could  potentially  improve  the  effectiveness  of 
the  controller.   Still,  Figures  5-6  and  5-7  serve  to  indicate  the 
improvements  achievable  with  PID  control. 

In  addition,  displays  similar  to  Figures  5-6  and  5-7  could  be  used 
in  conjunction  with  optimal  control  methods  in  establishing  the 
relative  weighting  criteria  that  would  be  assigned  to  deviations  from 
desired  or  setpoint  conditions  and  input  control  requirements.   These 
would  be  used  in  establishing  a  quantitative  definition  of  optimality. 
The  general  approach  of  using  modern  control  methods  in  the  control  of 
the  plasma  thermal  balance  is  presented  in  Appendix  B. 


CHAPTER  6 
CONCLUSIONS 


A  review  of  the  literature  on  the  problem  of  fusion  thermal 
stability  has  shown  that  many  control  mechanisms  have  been  proposed, 
implemented  and  analyzed  which  permit  operation  of  the  plasma  power 
balance  in  otherwise  unstable  temperature  regimes.   Though  these  control 
mechanisms  can  achieve  the  goal  of  stabilizing  the  thermal  balance,  the 
issues  of  overall  system  performance  have  not  been  addressed.   Certainly 
stability  is  a  necessary  condition  for  operation,  but  perturbations  and 
disturbances  associated  with  real-world  systems  additionally  require  that 
the  system  possess  desirable  response  characteristics,  such  as  acceptable 
overshoot  and  adequate  speed  of  response. 

This  work  represents  an  analysis  of  the  thermal  balance  of  a 
fusioning  plasma  from  a  control  system  perspective.   By  applying  the 
techniques  of  classical  control  theory,  in  addition  to  determining  the 
requirements  for  stability,  the  response  characteristics  of  the  thermal 
balance  have  been  evaluated.   Thus,  where  previous  studies  were  directed 
at  determining  the  critical  temperature  for  stability,  the  present  methods 
were  additionally  directed  at  obtaining  the  system  response  in  the  stable 
regime. 

The  response  of  the  thermal  balance  was  evaluated  using  standard 
inputs  to  permit  comparison  of  the  system  behavior  to  accepted  control 
system  guidelines.   The  results  show  that  the  stable  regimes  and 
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response  characteristics  of  open-loop  equilibria  are  highly  dependent 
on  the  governing  transport  law  and  the  injection  energy  of  the  incoming 
fuel.   With  constant  confinement,  operation  without  control  would  be 
limited  to  undesirably  high  temperatures.   Even  at  the  elevated 
temperatures,  the  open-loop  response  would  be  either  exceedingly 
sluggish  with  no  injection  energy,  or  would  exhibit  high  overshoots 
with  only  75  keV  injection  energy.   With  Bohm  confinement,  stable 
operation  is  possible  in  a  broader  temperature  range  with  or  without 
injection  energy,  but  overshoot  is  unacceptable  at  all  temperatures 
without  injection  energy,  and  also  at  temperatures  below  nearly  20  keV 
with  75  keV  injection. 

Although  it  may  be  possible  to  obtain  acceptable  system  response 
through  the  selection  of  the  proper  injection  energy,  this  mode  of 
operation  would  place  additional  and  unnecessary  constraints  on  plasma 
physics  considerations  and  reactor  operations.   Restrictions  on  the 
injection  energy  are  unnecessary  because  the  system  response  can  be 
improved  via  a  dedicated  control  system  employing  closed-loop  feedback. 
This  work  has  shown  that  closed-loop  proportional  feedback  using  the 
feedrate  and  injection  energy  can  be  used  to  extend  the  stable 
operating  temperature  range.   The  results  also  show  that  significant 
improvements  in  system  response  are  possible  using  only  proportional 
feedback. 

The  stabilized  temperature  regime  is  a  function  of  the  confinement 
lav/,  injection  energy,  and  feedback  coefficients.   A  range  of  feedrate 
and  injection  energy  feedback  fractions  is  available  with  which  to 
stabilize  the  fusion  power  balance.   Using  only  proportional  feedback 
and  no  injection  energy,  the  critical  temperature  decreases  to  13  keV 
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and  11  keV  with  constant  and  Bohm  confinement  respectively,  in  contrast 
to  the  corresponding  open-loop  values  of  25  and  14  keV.   With  injection 
energy,  stable  operation  is  available  to  as  low  as  5  keV,  the  ideal 
ignition  temperature,  by  proper  selection  of  the  feedback  coefficients. 

The  flexibility  in  selecting  the  feedback  coefficients  provides  a 
broad  range  of  control  over  the  system  performance.   With  constant 
confinement,  the  sluggish  open-loop  performance  with  no  injection 
energy  can  be  moderately  improved,  and  the  high  overshoots  with 
injection  energy  can  be  significantly  reduced.   However,  these 
improvements  are  possible  only  at  temperatures  above  20  keV  and  require 
high  feedback  fractions.  With  Bohm  confinement,  the  high  overshoots 
present  in  open-loop  operation  with  zero  injection  can  be  reduced  to 
acceptable  levels  at  temperatures  as  low  as  20  keV,  with  a  simultaneous 
decrease  in  the  settling  time  to  under  30  seconds.   With  75  keV 
injection,  acceptable  overshoot  can  be  obtained  at  temperatures  as  low 
as  10  keV,  with  the  time-to-peak  below  20  seconds  and  settling  times 
less  than  30  seconds.   Even  with  Bohm  confinement,  it  is  still 
difficult  to  simultaneously  satisfy  overshoot  and  speed  of  response 
requirements  at  low  temperatures  and  low  feedback  fractions. 

The  present  techniques  are  of  sufficient  accuracy  to  provide  an 
adequate  representation  of  the  actual  system  response  without  requiring 
dynamic  simulation  of  the  thermal  balance.   Thus,  the  techniques  can  be 
used  in  investigations  of  the  control  potential  of  other  mechanisms 
such  as  the  radiation-enhanced  control  of  Reference  17,  and  the 
promising  ripple-transport  control  mechanism  of  Reference  30.   The 
control  methodology  presented  here  can  be  then  used  to  ensure 
acceptable  temporal  behavior  of  the  plasma  output  power. 
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Though  the  system  response  is  significantly  improved  with  only 
proportional  feedback,  additional  improvement  is  available  with  the 
control  mechanisms  studied  in  this  work,  by  employing 
proportional-integral-derivative  (PID)  control.   This  work,  provides 
some  preliminary  indications  of  the  potential  benefits  of  PID  control, 
and  the  indications  are  that  more  efficient  control  is  possible  while 
simultaneously  improving  the  system  response.   Either  an  analytical  or 
parametric  analysis  would  be  suitable  for  further  examining  plasma 
behavior  with  PID  control. 

It  is  expected  that  the  complex  interrelationships  and  large 
number  of  variables  that  are  potentially  available  to  effect  control  of 
the  thermal  balance  of  a  fusioning  plasma  will  require  the  use  of 
modern  control  theory  and  optimization  techniques.   The  plasma 
equations  used  in  this  work  are  well  suited  for  adapting  to  state 
variable  form,  and  the  state  variable  representation  can  expand  to 
accommodate  additional  control  mechanisms  and  system  variables. 

If  fusion  achieves  its  ultimate  potential  as  a  viable  energy 
source,  the  control  of  the  power  balance  is  a  technical  area  which  will 
have  been  satisfactorily  addressed.   This  work  has  shown  how  only  two 
control  variables,  the  feedrate  and  injection  energy  proportional 
feedback  fractions  (and  variations  including  integral  and  derivative 
control)  can  significantly  change  plasma  performance.   It  is  hoped  that 
this  work  serves  as  a  framework  on  which  other  control  schemes  can  be 
analyzed  and  evaluated,  using  increasingly  sophisticated  techniques, 
which  lead  to  achieving  satisfactory  control  for  fusion  thermal 
stability. 


APPENDIX  A 

INTRODUCTION  TO  FUSION 

AND 

REVIEW  OF  THE  LITERATURE  ON  THERMAL  STABILITY 


A.l .   Introduction 

The  fusion  reaction  receiving  the  majority  of  attention  at  the 
present  time  is  the  reaction  of  deuterium  and  tritium  to  yield  a 
neutron,  helium,  and  energy  as  follows: 


2     3"+      1 

1°  +  1T  *  2He  +  0n  +  17*6  MeV 


The  energy  production  rate  per  unit  volume,  or  power  density,  Pf ,  of  a 
plasma  sustaining  this  fusion  reaction  can  be  expressed  as 

Pf  =  nQnT<av>E  A.l 

where  n^  and  n™  are  the  deuterium  and  tritium  fuel  ion  densities,  <ov> 
is  the  reaction  rate,  i.e.,  the  probability  of  fuel  ions  with  a  given 
velocity  distribution  of  undergoing  fusion  per  unit  time,  and  E  is  the 
energy  given  off  per  reaction,  17.6  MeV.   The  temperature  of  the  plasma 
enters  the  determination  of  fusion  power  density  via  the  reaction  rate 
and  the  relationship  between  the  plasma  temperature  and  the  particle 
velocity  distribution.   Typical  values  of  the  reaction  rate  for  a 
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Maxwellian  velocity  distribution  are  10    to  10"   cm  /sec  for 
temperatures  in  the  range  of  10  to  100  keV.   Thus,  from  equation  A.l, 
in  order  to  obtain  plasma  power  densities  comparable  to  present  day 
alternative  energy  sources  such  as  fission  power  or  coal,  e.g.,  in  the 
10  to  100  watts/cm  range,  fuel  ion  densities  on  the  order  of  10   to 
10  /cm  are  required. 

The  leading  approaches  to  achieving  such  densities  and 
temperatures  are  the  magnetic  confinement  schemes.   In  general,  these 
schemes  preferentially  direct  ions  along  magnetic  field  lines  and/or 
suppress  cross  flow  so  that  by  sufficiently  compressing  the  field 
lines,  or  sufficiently  suppressing  cross  flows,  high  ion  densities  can 
be  achieved.   A  figure  of  merit  often  used  to  evaluate  the  magnetic 
confinement  schemes  is  the  ratio  of  the  particle  kinetic  pressure  to 
the  magnetic  energy  density  (magnetic  pressure).   This  ratio  is  denoted 
by  3  and  is  expressed  as 


5"  A.2 

ii 

2y„ 


where  n  is  the  total  particle  density,  kT  is  the  kinetic  temperature,  B 
is  the  magnetic  field  strength,  and  u  is  the  magnetic  permeability 
constant  in  vacuum.   The  constant  g  is  a  characteristic  of  the 
confinement  method  effected  by  the  particular  machine  configuration  and 
hardware.   For  example,  the  Tokamak  type  of  machine  has  a  typical  8  of 
5%,  and  the  mirror  type  has  a  typical  g  of  50%. 
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Using  equations  A.l  and  A. 2,  and  assuming  that  the  plasma  has 

equal  densities  of  deuterium  and  tritium  ions  in  a  fully  ionized  state, 

nelectrons   n  .  .       ,    ,.  .     ,  c 

i.e.,  nn  =  nT  =  =  — ,  the  densities  can  be  eliminated  from 

D    r       2        4 

the  equations,  resulting  in  an  expression  for  the  power  density  in 
terms  of  machine  characteristics  and  temperature  only.   This  expression 
is  given  by 


p  =  E6V^ov>_  A>3 

U02   (kT)2 


For  a  given  fuel  cycle,  the  energy  released  per  fusion  event,  E,  is  a 

known  constant,  and  for  a  given  machine,  g  is  approximately  a  constant. 

Therefore,  for  a  given  maximum,  technologically-limited  magnetic  field 

strength,  the  power  density  is  optimized  by  optimizing  the  quantity 

.   The  optimum  occurs  in  the  vicinity  of  12  to  16  keV  and,  as 

(kT)2 

described  in  the  next  section  of  this  Appendix,  this  temperature  range 

has  been  found  to  be  thermally  unstable. 

An  appreciation  for  the  reasons  for  operating  at  the  temperatures 

determined  from  peak  considerations  can  be  obtained  as  follows. 

(kT)2 
The  most  comprehensive  description  of  a  commercial  fusion  reactor 

design  to  date  is  contained  in  the  STARFIRE  report[16].   The  STARFIRE 

reactor  operates  at  an  average  fuel  ion  temperature  of  24  keV,  has  a 

fusion  power  density  of  6.75  W/cm  ,  and  its  nuclear  fusion  island  cost 

is  estimated  at  $1.2  billion.   This  latter  figure  includes  only  those 

systems  necessary  to  confine  and  heat  the  plasma,  and  convert  the 

neutron  energy  into  the  currently  useful  form  of  heat.   The  net  cost  of 

electricity  was  estimated  to  be  35  mills/kWh.   If  the  temperature  were 
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increased  to  40  keV,  for  example,  and  assuming  that  magnetohydrodynamic 
(MHD)  plasma  stability  requirements  can  be  satisfied  without  additional 
penalty,  the  power  output  of  the  plant  would  decrease  by  45%,  resulting 
in  a  simultaneous  increase  in  the  cost  of  electricity  to  53  mills/kWh. 
This  increase  includes  credit  for  a  decrease  in  the  power  production 
equipment  cost.   Alternately,  the  net  power  output  can  be  maintained  by 
increasing  the  reactor  size  -  at  an  additional  cost  of  $326  million  for 
the  reactor  equipment  only.   By  comparison,  the  entire  instrumentation 
and  control  system  for  the  STARFIRE  complex  was  estimated  to  cost  $29 
million,  or  approximately  2%  of  the  original  nuclear  equipment  cost. 
Control  to  permit  operation  at  the  desired  optimum  condition  is  thus 
economically  justified. 

An  alternative  to  operating  at  lower  power  density  is  to  operate 
at  higher  magnetic  field.   In  the  STARFIRE  reactor,  the  superconducting 
toroidal  field  coils  are  designed  to  produce  an  on-axis  magnetic  field 
of  5.8  Tesla,  and  a  maximum  field  of  11.1  Tesla.   If  the  operating 
point  of  the  thermal  balance  were  increased  to  40  keV,  the  power 
density  can  be  maintained  at  the  previous  6.75  W/cm  by  increasing  the 
magnetic  field  strength.   If  the  reactor  maintains  the  same  value  of  g, 
according  to  equation  A. 3  a  new  field  strength  of  12.9  Tesla  is 
required  to  retain  the  previous  power  density.   This  is  beyond  the 
capability  of  anticipated  fabrication  and  metallurgical  advances  in 
superconducting  magnet  technology  during  the  next  two  to  three  decades. 
Thus,  the  rationale  has  been  established  for  desiring  to  operate  fusion 
reactors  at  temperatures  approaching  the  12  to  16  keV  range,  and  for 
exploring  methods  of  effecting  effective  control. 
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A.2.   Literature  Review 

As  presented  in  Chapter  1,  the  first  consideration  of  thermal 
stability  and  its  control  was  provided  by  Mills  [2].   In  Reference  2, 
Mills  formulated  the  following  particle  and  energy  equations, 


dn. 
IT  ' 


n±  ^}k(Ti  +  V  =  ni2p(l  _  PK<*v>DTcE*  -  |sk(Ti  +  Te) 


where  p  is  the  ratio  nm/nf ,  E  is  the  energy  of  the  alpha  particle,  c 
is  the  fraction  of  the  charged  particle  energy  retained  in  the  plasma, 
and  S  is  the  particle  feedrate.   Setting  temporal  derivatives  to  zero, 
Mills  obtained  the  conditions  for  steady  state.   To  study  the  stability 
of  the  equilibrium,  he  solved  his  steady  state  conditions  for  an 
expression  mathematically  equal  to  a  constant  and  proceeded  to  find  its 
logarithmic  variation.   Examination  of  the  behavior  of  the  resultant 
equation  indicated  that  the  equilibrium  was  stable  against  temperature 
fluctuations  when  the  plasma  temperature  was  above  28  keV.   He 
suggested  several  methods,  including  feedrate,  fuel  mixture,  radiation 
loss,  and  confinement  time,  by  which  control  might  be  exercised  for 
operation  below  28  keV.   He  demonstrated  the  method  which  he  deemed 
best,  confinement  time  control  inversely  proportional  to  the 
temperature  (Bohm-type  confinement) ,  to  achieve  stability  to  as  low  as 
7  keV. 
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This  work  was  performed  using  a  very  simple  model.   Only  the  D-T 
reaction  was  included  and  the  ion  species  were  assumed  to  be  either 
tritium  or  "everything  else".   Only  cold  fuel  was  injected,  fuel  burnup 
was  neglected,  and  radiative  losses  were  artificiallv  lumped  and 
evaluated  as  a  constant  fraction  of  the  fusion  energy  output.   The 
fusion  energy  was  assumed  to  thermalize  instantaneously,  and  the 
differences  in  the  electron  and  ion  temperatures  were  determined  in  a 
simple  fashion.   Still,  the  analysis  was  very  useful  in  obtaining 
results.   In  his  second  paper  [3],  Mills  extended  and  improved  this 
model  by  including  three  ion  species  -  deuterium,  tritium,  and  the 
product  alpha,  two  types  of  radiative  losses  -  brerasstrahlung  and 
synchrotron,  and  explicit  energy  transfer  between  all  particles. 
Using  this  second  model  Mills  repeated  his  Bohm  confinement 
control  and  found  that  the  stabilizing  influence  of  this  type  of 
control  depended  on  the  confinement  constant  assumed.   He  concluded 
that  mixture  control  was  a  better  control  method  and  presented 
time-dependent  behavior  demonstrating  stable  operation  with  a 
relay-type  (on-off)  application  of  mixture  control.   Nonetheless,  some 
confinement  time  control  was  necessary  to  prevent  unreasonable 
temperature  overshoot.   This  was  the  first  work,  to  be  performed  on 
fusion  thermal  stability. 

Ohta  et  al.  [4]  approached  the  problem  using  a  similar  model  but 
with  a  different  control  methodology.   The  equations  this  time  were 


in=  -  JL+  s 

dt     t 
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dnT  =  n2f  (T)   -  f£  +  STe 


dt 


TE 


Qr<av> 

f(T)   =  — : -   1.12xlO~°T1/2 

12 


where  n  (cm"  )  and  T  (keV)  are  the  density  and  temperature  of  the 
plasma.   Ions  and  electrons  were  assumed  to  be  at  the  same  temperature. 
The  particle  and  energy  confinement  times  are  distinguished  as  tn  and 

_o     _  1 

Tg.   The  fuel  injection  rate  is  S  (cm  sec  L)   and  the  injection  energy 
is  3T  (keV).   The  fuel  cycle  was  D-T,  and  the  ion  species  were  assumed 
to  be  50-50  deuterium-tritium.   Bremstrahlung  was  explicitly  included. 
The  confinement  time  was  assumed  to  be  proportional  to  n  T111,  with  &  and 
m  as  constants  defining  the  diffusion  type.   Different  constants  of 
proportionality  were  used  for  particle  and  energy  confinement. 
Synchrotron  radiation  was  neglected.   High  energy  injection  of  the 
incoming  fuel  was  included. 

The  nonlinear  particle  and  energy  equations  were  linearized  by  a 
perturbation  expansion  about  an  equilibrium  condition,  i.e.,  n  =  n  + 
n' ,  T  =  T  +  T' ,  etc.  The  results  were  linear  differential  equations 
in  the  perturbed  quantities.  The  following  equations  were  explicitly 
solved 


1  dn'  =  , J_    1  ,n'    _  J_  ;T 
nodt      Tno   Tnl  no   TT1  To 


1  dn'  +  1  dT'    , J_  _  1  ^n'    r_l_         1  si' 
nQdt    TQdt  "  lTl   ^JnQ    l"^"^JT0 
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—  -  nf 

■nl 


-213- 

-1 


0  3n  Jo 


1       ^n"1 
1        n 


TT1    °  3T  J° 


J_  =  _Lf2nf(T)  -  nT- 


-1 


Tl   ToL  3n 


T0   n„L    3T        3T  J° 


This  yielded  the  conditions  necessary  to  insure  that  perturbations  decay 
in  tirae,  thus  indicating  the  stability  of  the  equilibrium.   The  conditions 
are  Dresented  in  Table  A-l,  along  with  the  resultant  critical  temperatures, 
For  operation  in  the  unstable  regime,  Ohta  et_  a_l.  [41  suggested  feedback 
stabilization.   Confinement  time,  bremsstrahlung,  feedrate,  and  injection 
energy  were  suggested,  and  the  latter  two  were  deemed  preferable. 
Feedback  proportional  to  the  density  perturbation  was  ineffective  and 
feedback  proportional  to  the  temperature  perturbation  was  investigated 
in  detail.   The  feedback  term,  incorporating  a  delay  time,  was  added 
to  the  above  perturbation  equations  and  a  generalized  set  of  conditions 
necessary  for  feedback  stabilization  was  obtained.   The  conditions  are 
given  by 


Ay2  +  By  +  C  =  0 
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TABLE  A-l 
STABILITY  CRITERIA  AND  CRITICAL  TEMPERATURE  (Ref.  4) 


Confinement 
Scheme 


Stability  Criteria 


(a) 


Critical  Temperature 
(keV) 


Charged 

Particle 

Heating(b) 


150  keV 
Neutral  Beam 
Injection^0) 


t   =  constant 


^--r=^-    >  0 


28 


21 


T  .rl 


(Bohm) 


E 


(Classical) 


Fn    3T, 


>  -2  and  -  ■—- 
o  ^2 


F^lT^T  and  Ui   ^  "rJ 


14 


42 


33 


<a>F(T)  -  "-1- 


(b), 
(c) 


f(T) 

0,  R  =  1 

Ts  =  50  keV,  R  »  10 
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with 


A  =  1  +  (TS/TQ  -  DctAt 


1     1     1                  aTt^ 
B  =  -L---L  +  -i_-  (T  /T  -  l)a  +  


T        T  S'    0  T 

T2   TEo  TEo 


c  .  _L<-L  - -L)  -  f?i 

Tno  TEo   T2    TEo 


where 


52  -  1  - 


The  system  is  stable  if  A>0,  B>0,  C>0,  or  A<0,  B<0,  C<0.   Note  that  the 
conditions  for  stability  migrate  as  the  equilibrium  conditions  change. 

The  dynamic  behavior  of  a  fusioning  plasma  was  then  studied  by 
Horton  and  Kammash  [18]  using  a  model  similar  to  Mills'  second  model. 
The  strength  of  this  model  stemmed  from  its  refined  treatment  of  the 
confinement  laws  (three  Tokamak  reactor-type  scaling  laws  were 
included),  and  charged  particle  heating  (energy  transfer)  mechanisms. 
Fuel  burnup  was  included.   Stable  operating  regimes  for  a  range  of 
reactor  parameters  were  obtained  by  observing  the  time  dependent 
behavior  of  the  system.   The  three  confinement  schemes  investigated 
demonstrated  differing  characteristic  behavior-  negligible  overshoot 
with  Bohm  diffusion,  damped  oscillation  with  anomalous  diffusion,  and 
quasi-steady  (no  true  steady  state)  with  neoclassical  diffusion. 
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Powell  and  Hahn  [19]  added  a  novel  feature  to  the  previous  models 
-  the  interaction  of  the  plasma  with  the  confining  magnetic  field. 
Their  contention  was  that  results  obtained  previously  were  too 
pessimistic.   They  contended  that  as  the  temperature  of  a  plasma 
increases  due  to  a  runaway  (unstable)  condition,  the  plasma  would  tend 
to  expand,  thereby  doing  work  on  the  magnetic  field.   A  reverse 
situation  would  exist  if  the  temperature  were  decreasing.   This 
exchange  of  energy  in  the  form  of  work  might  serve  as  a  stabilizing 
influence  on  the  energy  balance,  leading  to  lower  critical 
temperatures. 

The  plasma  model  assumed  both  pure  deuterium  or  D-T  fuel 
injection,  allowance  for  the  four  dominant  fusion  reactions,  equal 
electron  and  ion  temperatures,  quasistatic  expansion  of  the  plasma,  and 
a  sharp  plasma-magnetic  field  boundary.   Power  input  to  the  magnetic 
coils  was  included  but  its  effect  was  investigated  in  a  subsequent 
study.   Volume  was  introduced  in  the  equations,  as  was  the  magnetic 
field  strength  (implicitly,  through  the  sharp  boundary  assumption). 
Results  were  presented  for  the  deuterium  injection  case  only.   A  stable 
energy  balance  was  obtained  for  temperatures  above  100  keV  and  the 
effect  of  including  the  magnetic  field  was  estimated  as  a  20  keV 
reduction  in  the  critical  temperature. 

In  the  subsequent  paper  [20],  Powell  significantly  extended  the 
results  and  conclusions  of  the  previous  model,  including  the  effect  of 
including  the  input  power  to  the  magnetic  field  coils.   He  first 
established  the  association  of  the  input  control  triplet  (B,  S~,  t), 
with  a  unique  equilibrium  point  in  n-T-V  space.   The  unique  association 
so  ascertained,  the  conclusion  was  drawn  that  field  strength,  feedrate, 
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and  confinement  time  constituded  a  satisfactory  set  of  control 
variables.   He  introduced  the  following  control  functions: 


«!L   =  C(Bf2  -  Bh 
dt      r 


ds 


2 
dt 


=  c(s2f  -  s2) 


£-  C(Tf  -  T) 
dt      r 


where  C  is  a  constant,  and  Be,  S2f,  and  t^  are  the  desired  values  of 
the  magnetic  field  strength,  deuterium  feedrate,  and  confinement  time. 
These  control  functions  were  employed  to  transfer  the  plasma  from  one 
equilibrium  state  to  another  without  instabilities  or  other 
misbehavior.   Overshoot  was  minor.   Other  conclusions  of  this  study 
were  that  scores  of  d3/namic  simulations  would  be  required  to  give  a 
complete  picture  of  the  system  response,  that  the  control  functions  may 
have  difficulties  at  low  temperatures,  and  that  optimization  of  the 
control  is  an  important  area  deserving  study. 

The  first  explicit  eigenvalue  analysis  of  the  thermal  balance  was 
presented  by  Stacey  [21].   The  plasma  model  included  the  D-T  fuel  cycle 
only,  two  ionic  species,  instantaneous  particle  slowing  down  and  energy 
equilibration,  explicit  energy  transfer  between  particles,  and  separate 
ion  and  electron  temperatures.   Confinement  time  and  radiation  losses 
were  modeled  in  great  detail.   The  equations  were  linearized  about  an 
equilibrium  point,  Laplace  transformed,  and  the  eigenvalue  dependence 
on  temperature  was  observed.   The  largest  eigenvalue  of  the  system  was 
shown  to  be  positive  at  low  temperatures,  becoming  negative  and 
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indicating  stability  at  temperatures  consistent  with  previous  results. 

The  results  were  confirmed  by  dynamic  simulation.   For  control  of  the 

plasma  at  an  unstable  point,  Stacey  investigated  both  feedrate  and  fuel 

mixture.   He  used  three  relay-type  algorithms  responding  to  the 

electron  temperature.   Delay  time  in  control  application  was  included. 

The  results  depicted  a  reactor  operating  in  an  oscilltory  mode,  with 

the  inability  to  restrict  the  oscillations  to  a  narrow  band  unless  the 

delay  was  reduced  to  zero. 

The  same  type  of  eigenvalue  analysis  performed  by  Stacey  was 

performed  by  Fujisawa  [22],  who  included  a  separate  energy  equation  for 

the  fast  alpha  group.   Energetic  injection  of  the  fuel  was  included  in 

the  model,  but  confinement  expressions  were  simpler  and  synchrotron 

radiation  was  ignored.   Fujisawa  linearized  his  equations  and  cast  them 

into  the  form  —  =  Ax  +  Bu,  where  A  and  B  are  the  system  matrices,  x  is 
dt 

the  state  vector  representing  particle  densities  and  energies,  and  u  is 
the  input  vector  (used  as  a  scalar  in  this  study).   The  system 
stability  was  examined  by  observing  the  eigenvalues  of  A.   For 
operation  in  the  unstable  region,  Fujisawa  investigated  the  effect  of 
external  feedback  control.   The  control  mechanism  investigated  was  a 
pure  capacitive  delay  of  the  fuel  injection  flow  rate.   The  effect  of 
this  control  on  the  dominant  eigenvalue  resulted  in  stabilized  behavior 
in  a  limited  temperature  range.   Additional  calculations  indicated 
improved  temporal  behavior  of  the  stabilized  plasma  by  the  use  of  a 
phase  advance  compensation  network. 

Bassioni  and  Husseiny  [23]  reverted  to  the  earlier  model  of  Ohta 
et  al.  [4]  and  examined  the  controllability  of  the  linearized  system. 
Fuel  burnup  was  included  as  the  only  revision  of  the  earlier  model. 


-219- 

Casting  the  linearized  system  into  the  form  of  Fujisawa,  the  earlier 
stability  criteria  of  Ohta  et^  al_.  [4]  shown  in  Table  A-l  were  rederived 
with  good  agreement.   In  addition,  the  matrix  G  =  [B  |  AB]  was  found  to 
be  non-singular.   This  indicated  the  controllability  of  the  system  with 
the  input  vector  u.   This  vector  consisted  of  particle  injection, 
external  heating,  and  bremsstrahlung  radiation  loss  control.   Invoking 
a  constant  feedback  law  including  time  delay,  the  characteristic 
equation  of  the  system  was  then  obtained.   Using  the  Pontryagin 
criterion,  the  domain  of  linear  stability  was  found.   More  general 
results  were  derived  using  the  Mikhailov  criterion. 

The  results  obtained  in  Reference  23  are  very  general.   Due  to 
their  generality,  it  is  difficult  to  translate  their  conditions  into 
critical  temperatures.   On  the  other  hand,  Usher  and  Campbell  [24,25], 
using  essentially  the  model  of  Ohta  £t  al.  [4],  and  including  burnup , 
reaction  product  accumulation,  and  beam-plasma  interactions,  clearly 
delineated  t^e  critical  temperatures  for  various  fuel  cycles  and 
confinement  models.   Their  findings  are  presented  in  Table  A-2.   The 
differences  from  Table  A-l  are  due  to  improvements  in  the  plasma 
modeling.   For  operation  below  the  critical  temperature,  delayed 
feedback  control  was  employed.   The  results  permitted  the 
characterization  of  the  confinement  schemes  as  to  ease  of  control, 
required  magnitude  of  the  feedback,  and  the  allowable  delay  times. 
Maya  [14]  then  quantitatively  demonstrated  that  the  use  of  the  type  of 
feedback  of  Ohta  et  al.  [4]  was  limited  as  to  how  far  below  the 
critical  temperature  stabilized  operation  was  possible.   O'Farrill  and 
Campbell  [11]  then  presented  critical  temperatures  for  two  component 
plasmas.   They  used  an  extensive  model  for  the  beam-plasma  interaction 
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TABLE  A-2 

CRITICAL  TEMPERATURES  (keV)  FOR 
DIFFERENT  FUEL  CYCLES  (Ref.  24) 


Fuel 
Cycle 

BOHM 

Constant 

Classical 

Maximum 
Temperatures 

D-T 

14 

24 

38 

200 

D-He3 

68 

98 

140 

660 

D-D 

40 

82 

260 

>1000 
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and  extended  previous  results  by  including  non-Maxwellian  fusion 
reactions.   Table  A- 3  summarizes  the  extended  ranges  of  operation 
resulting  from  the  beam-plasma  interaction. 

Nonlinear  techniques  were  applied  to  the  analysis  of  thermal 
stability  by  Ferrell  and  Kastenberg  [27],   Using  the  model  of  Ohta  et_ 
al.  [4],  the  Mean  Value  Theorem  was  used  to  cast  the  system  into  a  form 
readily  analyzed  by  comparison  theorems.   Upper  and  lower  bounds  of 
stability,  as  well  as  the  domain  of  allowable  perturbations  were  thus 
obtained. 

A  multigroup  formulation  of  the  slowing  down  of  the  alpha 
particles  was  included  in  the  work  of  Ohnishi  and  Wakabayashi  [28]. 
The  energy  transfer  of  the  product  alphas  and  energetic  injected  ions 
to  the  thermal  background  plasma  was  treated  in  great  detail.   The 
linearized  state  equation  was  obtained  as  before  and  the  eigenvalues  of 
the  system  were  obtained.   The  effect  of  the  instantaneous  slowing  down 
approximation  could  thus  be  evaluated.   The  results  showed  that  in  the 
unstable  regions,  the  multigroup  approximation  resulted  in  nearly  a  30% 
reduction  in  the  magnitude  of  the  largest  eigenvalue  over  that  obtained 
assuming  instantaneous  slowing  down.   These  results  demonstrated  the 
stabilizing  influence  of  the  finite  time  required  by  the  fusion  alphas 
to  transfer  their  energy  to  the  bulk  background  plasma. 

The  use  of  toroidal  field  ripple  as  a  control  mechanism  of  the 
thermal  balance  was  proposed  by  Petrie  and  Rawls  [29].   Magnetic  field 
ripple  is  the  non-axisymmetric  component  of  the  toroidal  magnetic  field 
created  by  the  finite  number  of  coils  used  to  generate  the  confining 

magnetic  field.   The  resultant  ripple-enhanced  thermal  conductivity, 

7/? 
with  its  T  -  dependence,  was  proposed  as  an  effective  method  of 
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TABLE  A-3 

EFFECT  OF  THE  ENERGETIC  DEUTERON  FRACTION 

ON  THE  STABLE  OPERATING  BULK  TEMPERATURES 

OF  THE  TCT  FOR  THREE  CONFINEMENT  MODELS  (Ref.  11) 


Stability  Ranged 
Diffusion  Model        w  (keV) 

T_  a  Constant 


a   »  0,  m  =  0 


Bohm 


0.0 

27  to  30 

0.1 

1  to  5  and  25  to  30 

0.2 

1  to  7  and  23  to  30 

0.3 

1  to  11  and  18  to  30 

0.4 

to  1.0 

1  to  30  all  stable 

TEaT-l 

0.0 

1: 

>  to  30 

I   =  0,  m  =  -1 

0.1 

1 

to 

30  all  stable 

Classical 

TE  a   n-lTl/2 

0.0 

All  unstable 

I   =  -1,  m  =  1/2 

0.1 

1 

to 

3 

0.2 

1 

to 

4 

0.3 

1 

to 

6 

0.4 

1 

to 

7  and  29  to  30 

0.5 

1 

to 

9  and  26  to  30 

0.6 

1 

to 

15  and  19  to  30 

0.7 

to  1.0 

1 

to 

30  all  stable 

Deuteron  injection  energy  =  150  keV 
nB 


nT 


- 1  «•■  »  "  urhjy  R  "  °-°16 


(a) i   kev  to  30  keV  was  the  range  of  temperatures 
studied. 
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inhibiting  the  runaway  growth  of  the  plasma  temperature  by  increasing 
energy  loss  faster  than  the  energy  gain  from  the  increase  in 
reactivity.   A  set  of  conditions  determined  from  physics  and 
engineering  considerations  were  derived  in  Reference  29  specifying  the 
range  of  ripple  required  to  provide  satisfactory  burn  control.   The 
conditions  were  specified  in  order  to  ensure  that  the  ripple-enhanced 
transport  does  not  appreciably  increase  the  difficulty  of  reaching 
ignition,  that  it  is  the  primary  loss  mechanism  during  burn,  and  that 
the  increased  loss  of  fast  alphas  does  not  degrade  the  reactor 
performance.   The  work  was  performed  using  a  0-D  model,  and 
one-dimensional  transport  analysis  results  were  shown  to  be  in  good 
agreement. 

The  use  of  ripple-enhanced  transport  to  effect  thermal  stability 
was  also  investigated  by  Bromberg  and  Cohn  [30].   As  in  Reference  29, 
since  stability  is  achieved  by  enhancing  an  energy  loss  mechanism,  a 
higher  nx  condition  is  necessary  to  achieve  ignition.   In  Reference  30, 
thermal  stability  by  radial  motion  of  the  plasma  as  determined  by  MHD 
constraints  was  also  investigated.   Physically,  the  radial  motion 
accompanying  a  runaway  temperature  condition  reduces  the  plasma  density 
and  thus  the  alpha  heating  power,  minimizing  the  thermal  excursion.   A 
critical  temperature  as  low  as  25  keV  was  obtained  without  resorting  to 
any  anomalous  ion  loss  mechanism. 

The  use  of  low  energy  neutral  fuel  injection  at  the  plasma  edge 
("gas  puffing")  as  a  thermal  control  mechanism  was  investigated  in 
Reference  31.   This  type  of  control  enjoys  the  advantage  of  requiring 
negligible  auxiliary  power  in  that  some  form  of  gas  injection  for 
refueling  of  the  reactor  during  operation  is  postulated  for  many 
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reactor  designs.   It  effects  control  of  the  thermal  instability  by 
controlling  the  plasma  ion  density  via  diffusion  of  the  injected 
neutrals  into  the  central  reactor  region.   Though  the  oscillations  in 
reactor  power  ranged  as  high  as  20%  in  the  numerical  examples  presented 
in  the  articles,  the  authors  presented  arguments  stating  that  the 
excursions  could,  in  principle,  be  reduced  to  zero  amplitude. 
The  literature  review  presented  herein  describes  all  the 
mechanisms  which  have  been  proposed  for  the  control  of  the  thermal 
instability.   Emphasis  has  been  placed  on  0-D  models  since  even  most  of 
the  analyses  reported  in  the  literature  which  began  with  1-D  models 
averaged  the  plasma  parameters  over  the  appropriate  dimension  in  order 
to  be  able  to  obtain  manageable  mathematics  and  quantifiable  results. 
Illustrative  examples  have  been  presented  where  appropriate.   However, 
this  review  is  not  exhaustive  in  numerical  examples  since  the 
applications  to  specific  reactor  conditions  do  not  provide  any 
additional  insights.   The  view  that  has  been  presented  of  the  thermal 
stability  problem  is  one  in  which  many  mechanisms  have  been  proposed 
that  provide  varying  degrees  of  control.   However,  no  formal  method  of 
specifing,  evaluating,  and  selecting  the  preferred  choice  has  been  put 
forth.   This  issue  is  addressed  in  the  main  body  of  this  report. 


APPENDIX  B 
CONTROL  THEORY  REVIEW 


B.l.   Introduction 


The  most  fundamental  function  of  a  control  system  is  to  adjust  the 
operating  state  of  a  process,  as  determined  by  the  measured  values  of  a 
set  of  variables,  to  match  a  set  of  desired  values  for  the  same  set  of 
variables.   This  can  be  accomplished  by  using  either  open-  or 
closed-loop  control.   An  open-loop  system  is  one  in  which  the 
controller  action  is  not  dependent  on  the  actual  response  of  the 
controlled  system.   A  closed-loop  control  system  is  one  in  which  the 
deviation  of  the  actual  operating  condition  from  the  desired  operating 
condition  determines  the  course  of  action  taken  by  the  controller  to 
return  the  system  to  the  desired  state. 

The  thermal  balance  of  a  fusioning  plasma  can  be  viewed  as  an 
open-loop  control  system  as  follows.   The  values  of  the  input 
parameters  required  to  maintain  a  prescribed  level  of  energy  production 
in  the  plasma  can  be  obtained  by  solving  the  appropriate  governing 
equations  of  mass  and  energy  balance.   In  theory,  the  steady  state 
value  of  the  system  variables  is  dictated  by  the  choice  of  input 
parameters.   In  practice,  however,  external  disturbances  are  quite 
likely  to  cause  deviations  from  the  nominal  operating  conditions.   In 
this  open-loop  system,  the  input  parameters  would  not  be  automatically 
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ad justed  to  offset  the  deviations,  since  the  input  controller  is  not 
affected  by  the  system  response.   Chapter  3  and  the  previous  Appendix 
have  shown  how  the  open-loop  plasma  thermal  balance  is  unstable  to  such 
disturbances  in  the  temperature  range  of  12  -  16  keV,  although 
stability  could  subsequently  be  achieved  at  higher  temperatures.   Even 
at  the  higher,  stable  temperatures  the  thermal  balance  was  shown  to 
have  poor  system  response. 

This  open-loop  system  can  be  transformed  into  a  closed-loop  system 
through  the  addition  of  feedback,  i.e.,  by  modifying  the  input 
(control)  variables  in  accordance  with  the  behavior  of  the  system 
variables.   The  closed-loop  plasma  system  may  then  be  able  to  sustain 
stable  performance  at  temperatures  which  were  previously  unstable.   The 
closed-loop  system  may  also  be  able  to  improve  the  system  response  at 
temperatures  which  were  previously  stable,  but  had  poor  response 
characteristics . 

Other  desirable  features  of  a  control  system  are  that  it  can 
increase  the  ability  of  the  overall  system  to  tolerate  spurious  signals 
and  noise,  i.e.,  it  can  decrease  the  system  sensitivity  to 
perturbations,  and  that  it  can  return  perturbations  of  the  system 
variables  to  their  nominal  values  in  a  reasonable  time,  i.e.,  that  it 
result  in  acceptable  overall  system  response.   A  review  of 
representative  system  specifications  of  these  and  other  desirable 
features,  in  both  the  time  and  frequency  domain,  which  are  used  to 
evaluate  system  performance  are  presented  in  this  Appendix.   Many  of 
these  desirable  characteristics  can  be  achieved  with  closed-loop 
control  systems. 
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The  analysis  and  design  of  control  systems  can  be  divided  into  two 
general  categories,  classical  and  modern  control.   Whether  classical  or 
modern,  the  analysis  phase  consists  of  certain  basic  steps  which 
include  gaining  a  thorough  understanding  of  the  physical  processes 
involved,  describing  the  system  with  a  sequence  of  functional  block, 
diagrams  and/or  differential  equations,  the  application  of  the  analysis 
technique  (whether  classical  or  modern),  and  the  comparison  of  the 
results  of  the  analysis  with  intuition,  past  experience,  and  the 
physical  understanding  of  the  system. 

In  the  classical  control  system  design  method,  if  the  results  of 
the  analysis  reveal  inadequacies  in  system  performance,  equipment  or 
components  can  be  added  to  the  system  in  order  to  improve  the  responce 
and  obtain  the  desired  specifications.   The  process  is  iterative  in 
that  the  analysis  phase  is  reentered  and  the  response  of  the  plant  is 
reevaluated,  including  the  effect  of  the  added  compensation  elements. 
The  latter  elements  are  then  readjusted  in  light  of  the  updated 
analysis  results.   The  iteration  is  terminated  and  the  control  system 
design  considered  complete  when  acceptable  behavior  is  obtained  as 
measured  by  the  set  of  performance  specifications. 

In  the  modern  control  system  design  method,  a  performance  index  is 
first  defined  for  the  system.   Then  the  necessary  structure  to  obtain 
or  minimize  the  specified  performance  index  is  determined  by  following 
a  logical  mathematical  formulation  which  completely  and  uniquely 
satisfies  the  constraints.   These  methods  are  described  in  this 
Appendix,  beginning  with  the  classical  methods. 
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B.2.   Classical  Methods 

Consider  the  linear  open-loop  system  of  Figure  B-l,  and  the 
closed-loop  system  of  Figure  B-2.   The  functions  r(t)  and  c(t) 
represent  the  system  input  (control)  and  output  functions.  The  system 
feedforward  transfer  function,  defined  as  the  Laplace  transform  of  the 
system  impulse  response  for  zero  initial  conditions,  is  given  by  G(s). 
H(s)  represents  the  feedback  transfer  function.   The  open-loop  transfer 
function  is  defined  as  G(s)H(s),  while  the  closed-loop  transfer 
function  is  given  by  the  expression 

C(s)  m  G(s) 

R(s)    1  +  G(s)H(s) 


For  second  order  systems,  the  denominator  of  the  transfer  function 

C(s) 

,  called  the  characteristic  polynomial,  when  set  equal  to  zero, 

R(s) 

yields  the  system  characteristics  of  natural  frequency,  damping  ratio, 
and  the  closed-loop  poles  of  the  system.   These  characteristics  dictate 
the  transient  response  of  the  system. 

The  evaluation  of  system  performance  in  the  time  domain  involves 
observing  the  response  of  the  system  to  standard  inputs  such  as  an 
impulse-  or  step-function.   The  specifications  used  to  judge  the 
analysis  include: 

1.  Maximum  overshoot,  M  -  the  magnitude  of  the  system  response 
following  a  step  input,  at  the  peak  of  the  first  overshoot. 

2.  Time-to-peak,  t  -  the  time  required  to  reach  the  first  peak 
(maximum  overshoot) 
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Figure  B-l.   Open-loop  system. 
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3.  Time  to  first  zero,  t  -  the  time  required  to  reach  the  final 
value  for  the  first  time. 

4.  Settling  time,  t  -  time  for  the  output  response  to  remain 
within  a  specified  percentage,  usually  5%,  of  the  final  value. 

In  the  classical  second  order  system,  the  above  performance 
specifications  can  be  explicitly  expressed  in  terms  of  the  system 
characteristics.   (The  specifications  for  the  plasma  system  are  derived 
in  Appendix  C.)  When  the  real  parts  of  all  the  poles  are  less  than 
zero,  the  transient  response  decays  exponentially,  i.e.,  the  system  is 
stable.   If  any  of  the  poles  have  positive  real  parts,  the  response 
grows  with  time,  indicating  instability.   To  study  the  response  of  the 
classical  second  order  system,  it  is  therefore  necessary  to  study  only 
the  system's  poles. 

The  open  loop  transfer  function  of  Figure  B-l  has  a  constant  gain, 
K,  associated  with  it.   If  this  gain  is  allowed  to  vary,  the  roots  of 
the  characteristic  equation  (1  +  KGH)  will  vary.   The  root  locus  method 
provides  a  graphical  means  of  evaluating  the  roots  of  the 
characteristic  equation  as  the  loop  gain  varies.   Basically,  a  set  of 
construction  rules  [32  -  34]  are  followed  that  permit  the  plotting  of 
the  root  loci.   Mechanical  aids  and  computer  programs  [35]  are 
available  that  provide  accurate  root  locus  plotting.   An  unstable 
linear  system  is  stabilized  by  adjusting  the  gain  such  that  all  the 
roots  have  crossed  the  imaginary  axis  and  entered  the  left-half  plane. 
If  stability  is  not  attainable,  additional  elements  or  modifications 
are  incorporated  into  the  control  system. 
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For  the  analysis  and  synthesis  of  control  systems  in  which  the 
specifications  are  specified  in  the  time  domain,  the  root  locus 
technique  is  an  able  tool.   However,  frequency  response  methods  are 
often  naturally  applicable  due  to  the  presence  of  sinusoidally  varying 
time  signals.   Quite  often,  the  desired  performance  of  a  system  can  be 
translated  into  frequency  domain  specifications  such  as  bandwidth, 
resonance  peak  and  frequency,  and  gain  and  phase  margins.   These  are 
defined  below  for  the  case  of  second  order  systems. 

1.  Bandwidth,  u)g-  the  frequency  where  the  closed  loop  response 
magnitude  is  1//2  of  its  low  frequency  magnitude. 

2.  Resonance  frequency,  u  -  the  frequency  at  which  the  transfer 
function  obtains  its  maximum  magnitude. 

3.  Resonance  peak,  M_-  the  magnitude  of  the  transfer  function  at 
the  resonance  frequency. 

4.  Gain  margin,  GM-  At  the  frequency  where  the  phase  of  the  open 
loop  transfer  function  G(juj)H(ju>)  passes  through  -180°,  GM  is  the 
additional  gain  needed  to  make  the  system  unstable  if  it  is  already 
stable,  or  the  amount  by  which  the  gain  must  be  reduced  to  make  the 
system  stable  if  it  is  already  unstable. 

5.  Phase  Margin,  PM-  At  the  frequency  at  which  the  magnitude  of 
G(ja))H(ju))  is  unity,  PM  is  the  additional  phase  lag  needed  to  make  the 
system  unstable  if  it  is  already  stable,  or  it  is  the  additional  phase 
gain  needed  to  make  the  system  stable  if  it  is  already  unstable. 

Three  of  the  more  common  methods  of  plotting  the  frequency 
response  of  a  system  include  polar  plots,  Bode  plots,  and  gain-phase  or 
log-modulus  plots.   It  is  convenient,  in  plotting  frequency  response, 
to  substitute  jco  for  the  Laplace  transform  variable,  s.  The  polar  plot 


-232- 
is  useful  in  connection  with  the  Nyquist  stability  criterion,  discussed 
below.   Both  the  gain-phase  plot  as  well  as  the  polar  plot  can  be 
obtained  from  the  Bode  plot,  so  only  the  latter  will  be  discussed. 

The  Bode  plot  is  simply  a  graph  of  the  magnitude  and  phase  of  the 
transfer  function  versus  frequency.   Also  called  the  corner  or 
logarithmic  plot,  the  Bode  plot  applies  the  property  of  logarithms  of 
changing  multiplications  into  additions.   Because  of  this  simplicity, 
the  addition  of  poles  and  zeros  to  the  transfer  function  (in  a  design 
step,  for  instance)  is  quickly  incorporated  in  the  graph.   A  gain 
adjustment  merely  raises  or  lowers  the  magnitude  plot  at  every  point  a 
fixed  amount,  while  the  phase  plot  remains  the  same.   The  construction 
rules  for  the  Bode  plot  are  quite  simple  and  straightforward  and  are 
also  available  in  References  32  to  35. 

The  Nyquist  stability  criterion  is  a  graphical  method  of 
predicting  the  stability  of  the  closed  loop  system  using  the  frequency 
response  of  the  open  loop  transfer  function.   Consider  a  closed 
(clockwise)  contour  Q  enclosing  the  whole  right-half  plane  in  s-space. 
Consider  now  a  polar  plot  of  G(s)H(s),  as  s  follows  the  contour  Q.  The 
Nyquist  criterion  states  that  the  net  number  of  (counterclockwise) 
encirclements  of  the  point  (-1,0)  by  the  plot  of  GH  is  equal  to  the 
number  of  poles  of  GH  minus  the  number  of  poles  of  the  closed  loop 
system  inside  the  contour  Q,  i.e.,  in  the  right-half  s-plane.   For 
stability,  it  is  required  that  there  be  no  closed-loop  poles  in  the 
RHP,  therefore,  the  number  of  counterclockwise  encirclements  of  (-1,0) 
should  equal  the  number  of  poles  of  GH. 
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Whereas  numerical  stability  techniques  such  as  the  Routh  criterion 
tell  only  if  a  system  is  stable  or  not,  Bode  and  Nyquist  plots 
additionally  give  a  measure  of  the  relative  stability  of  the  system. 
Relative  stability  is  measured  by  gain  and  phase  margin. 

After  the  specifications  have  been  assessed  for  the  given  system, 
if  the  desired  specifications  have  not  been  met,  gain  adjustments  and 
additional  frequency-dependent  compensation  is  designed  into  the 
system.   It  is  in  this  design  phase,  that  all  the  specifications  must 
be  simultaneously  satisfied,  if  possible.   This  is  accomplished  by  the 
addition  of  physically  realizable  components  in  one  of  two  ways-  in 
cascade  (in  the  feedforward  path),  or  in  the  feedback  loop.   Two 
examples  of  compensation  elements  are  the  phase  lag  and  the  phase  lead 
networks.   Other  types  of  compensation  elements  can  be  found  by  the  use 
of  Truxal's  method,  inverse  root  locus  method,  Chen's  method,  etc. [32] 
These  methods  involve  trial  and  error  procedures,  experience, 
ingenuity,  and  a  great  deal  of  perseverance.   The  satisfactory  system 
design  is  by  no  means  unique. 

The  number  of  techniques  available  with  which  to  effect  control  of 
a  system  is  large.   Table  B-l  presents  a  partial  list  of  generic 
control  types.   Each  control  type  has  certain  advantages  and 
disadvantages,  and  the  list  is  structured  in  the  order  of  increasing 
control  system  sophistication.   These  are  discussed  below,  including 
examples  of  their  applications  to  the  control  of  the  plasma  thermal 
balance. 

The  on-off  controller  is  the  simplest  type  of  control.   The 
control  input  only  has  two  states  in  this  mode  of  operation,  either 
fully  on  or  fully  off,  depending  on  the  value  of  a  predetermined  system 
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Table  B-l 
GENERAL  MECHANISMS  AVAILABLE  FOR  EFFECTING  FEEDBACK  CONTROL 


On-off  (Two-position  control) 

On-off  with  time  delay  or  differential-gap 

Multiposition  control 

Proportional 

Integral 

Derivative 

Proportional  plus  integral 

Proportional  plus  derivative 

Proportional  plus  integral  plus  derivative 
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variable  in  relation  to  its  desired  value.   This  control  is  thus  often 
called  two-position  control.   Its  advantages  include  its  simplicity  and 
practicality.   Its  disadvantages  include  indefinite  cycling,  and  the 
inability  to  ever  achieve  an  exact  correction  or  a  stable,  balanced 
condition.   Variations  of  the  basic  on-off  control  which  attempt  to 
reduce  the  cycling  frequency  and  the  magnitude  of  the  oscillations  are 
the  multiposition  control,  and  on-off  with  time  delay.   This  type  of 
control  was  attempted  for  control  of  the  plasma  thermal  balance  using 
the  fuel  ion  injection  ratio  as  the  control  variable  in  References  3 
and  21,  as  discussed  in  Appendix  A.  Oscillatory  behavior  was  indeed 
observed  in  dynamic  simulations.   In  summary,  cycling  and  overshoot  are 
characteristics  of  the  on-off  control  mode,  and  their  effects  are 
unacceptable  in  the  present  application. 

Proportional  control  offers  a  continuous,  linear  relation  between 
the  adjustment  of  the  input  parameter  and  the  deviation  from  the 
desired  operating  condition,  i.e.,  the  feedback  fraction  is 
proportional  to  the  perceived  error.   This  type  of  control  has  been 
used  by  many  researchers  of  thermal  stability  with  good  results  in 
terms  of  obtaining  a  stable  thermal  balance  (see  Appendix  A).   The 
advantages  of  proportional  control  are  its  powerful  stabilizing  action, 
and  its  reduction  of  the  system  sensitivity  to  external  perturbations. 
Its  disadvantage  is  that  proportional  control  is  balanced  for  only  one 
set  of  plasma  conditions.   Thus,  for  permanent  departures  from  this 
condition,  whether  it  be  planned,  such  as  a  load  change,  or  the  result 
of  component  drift,  proportional  control  will  only  approach  the  desired 
condition,  but  will  not  reach  it.   This  can  be  understood  as  follows. 
To  reach  equilibrium  at  the  new  condition,  a  change  in  the  input  is 
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required.   However,  with  proportional  control,  the  change  in  the  input 
is  proportional  to  the  perceived  deviation  from  the  desired  condition. 
Thus,  if  the  system  did  reach  the  new,  desired  condition,  the  change  in 
the  input  would  decrease  to  zero.   The  net  result  is  that  an  offset 
from  the  desired  condition,  or  an  error,  will  always  be  present  with 
proportional  control  except  for  one  reference  set  of  conditions.   This 
error  can  be  eliminated  by  manually  readjusting  the  input  variables  or 
the  new  set  point.   However,  this  manual  readjustment  would  need  to  be 
repeated  for  each  load  change,  and  this  procedure  would  not  be 
practical  in  a  complex  application.   An  alternative  method  of 
eliminating  offset  is  to  use  integral  control. 

With  integral  control,  the  magnitude  of  the  change  effected  on  the 
input  parameter  is  proportional  to  the  integral  of  the  error.   Thus, 
there  is  no  offset  error  with  integral  control  since  the  control  input 
will  be  continuously  adjusted  until  the  error  remains  at  zero.   The 
disadvantage  of  integral  control  is  that  its  initial  response  may  be 
slow,  since  it  may  take  some  time  for  the  integral  of  the  error  to 
increase  to  a  significant  value.   For  this  reason,  integral  control  is 
often  combined  with  proportional  control  to  improve  its  speed  of 
response. 

To  improve  the  response  of  slow  systems,  derivative  control  is 
often  used.   Derivative  control  is  accomplished  by  making  the 
adjustments  to  the  control  input  proportional  to  the  rate  of  change  of 
the  error.   Thus,  in  a  sense,  derivative  control  predicts  the  error  and 
compensates  accordingly,  often  reducing  the  maximum  error.   Since  the 
compensation  of  derivative  control  is  such  so  as  to  oppose  all  change, 
this  type  of  control  has  a  great  stabilizing  effect  of  the  system. 
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However,  since  this  type  of  control  only  affects  the  input  during  a 
change  in  the  error,  derivative  control  alone  can  not  correct  a  steady 
state  error.   It  is  therefore  usually  coupled  with  either  proportional 
control,  to  give  what  is  known  as  proportional-plus-derivative  control 
action,  or  with  proportional  and  integral  control,  to  give  the 
proportional- integral-derivative  (P-I-D)  controller.   This  latter 
control  scheme  combines  the  advantages  of  the  three  previous  control 
modes  in  that,  following  a  system  disturbance  or  load  change,  in 
addition  to  the  temporary  contribution  to  the  change  in  input  from  the 
proportional  control,  integral  control  introduces  the  required 
permanent  change  in  the  input  to  give  zero  steady  state  error; 
derivative  action  gives  good  initial  speed  of  response,  decreases 
overshoot  and  cycling,  and  enhances  system  stability.   Even  with  the 
current  advancements  in  modern  control  theory,  the  P-I-D  controller 
remains  as  a  standard  on  which  comparisons  of  control  system 
performances  are  based. 

The  original  development  of  classical  control  theory  has  given  way 
since  the  1960's  to  theoretical  investigations  into  better  approaches 
to  control.   With  the  aid  of  computers,  microprocessors,  and  software 
developments,  hardware  limitations  have  been  removed,  and  the  advanced 
control  techniques  are  finding  widespread,  practical  applications. 
Modern  control  theory  analysis  and  design  methods  are  covered  in  the 
next  section,  which  begins  with  a  discussion  of  the  notion  of  state 
space. 
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B.3.   Modern  Methods 

The  state  of  a  system  is  defined  by  Kalman  [36]  as  follows: 

State.   The  state  of  a  system  is  a  mathematical  structure 
containing  a  set  of  n   variables  x^Ct),  x2(t),  **',  xn(t),  called  the 
state  variables,  such  that  the  initial  values  xAt    )  of  this  set,  and 
the  system  inputs  u.(t)  are  sufficient  to  uniquely  describe  the 
system's  future  response,  for  t  >   t  . 

From  this  definition,  the  n-dimensional  state  vector  is  formed 
whose  components  are  the  state  variables,  the  state  space  is  defined  as 
the  n-dimensional  space  in  which  the  components  of  the  state  vector 
represent  the  coordinate  axes,  and  the  state  trajectory  is  the  path 
produced  in  state  space  by  the  vector  with  the  passage  of  time.  It  is 
noted  that  there  is  no  unique  way  of  selecting  the  system  variables 
with  which  to  represent  the  state  of  the  system. 

Consider  the  linear  dynamical  system  represented  by 

x(t)  -  A(t)x(t)  +  B(t)u(t)  B.l 

v(t)  =  C(t)x(t)  +  D(t)u(t)  B.2 

where  x,    u,    and  y_   are  the  state,  input,  and  output  vectors  of 
dimensions  n,  r,  and  m,  respectively,  and  A,  B,  C,  and  D  are  n  x  n,  n  x 
r,  and  ra  x  r  matrices.   Also  consider  the  matrix  differential  equation 

X  =  A(t)X  B.3 

where  X  and  A  are  both  n  x  n  matrices,  and  the  i-th  column  of  X(t)  is 
the  solution  to 
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x  =  A(t)x  B.4 

subject  to  the  initial  condition  given  by  the  i-th  column  of  X  .   A 
solution  of  B.3,  X(t),  subject  to  the  initial  condition  X(t  )  =  XQ,  is 
called  a  fundamental  matrix  of  B.4,  provided  X  is  nonsingular.   The 
solution  of  B.3  subject  to  X(tQ)  =  I,  where  I  is  the  identity  matrix, 
is  called  the  state  transition  matrix  of  B.4  and  is  denoted  as  T( t , t  ), 
Thus  the  following  holds 


T(t,tQ)  =  A(t)T(t,tQ),  T(tQ,t0)  =  I 


The  state  transition  matrix  has  several  interesting  properties, 
Among  these  are[37]: 

1.  T(t,tQ)  is  unique  for  all  t. 

2.  The  solution  of  B.4,  with  x(tQ)  =  xQ  is  x(t)  =  T(t,t0)xQ. 

3.  For  all  t,  tQ,  and  tp  T(t  .t^K  tj  ,tQ)  =  T(t,tQ). 

4.  T(t,tQ)  is  nonsingular  for  all  t,  and  T_1(t,t0)  =  T(t  ,t). 

5.  If  X(t)  is  any  fundamental  matrix  of  B.4,  then  T(t,t  )  = 

x(t)x"1(to). 

Employing  the  state  transition  matrix,  the  state  transition 
function  of  the  system  represented  by  B.l  and  B.2  is  given  by[37] 


x(t)  =  TCt.t^XQ  +  /   T(t,t')B(t')u(t')  dt' 


The  output  vector  y(t)  is  then  given  by 


y(t)  =  C(t)T(t,t0)xQ  +  C(t)J   T(t,t')B(t')u(t')  dt'  +  D(t)u(t) 
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and  the  x(t)    given  above  represents  the  unique  state  trajectory  through 
the  point  Xq  •   The  C(  t)T(  t  ,t0)xQ  term  in  y<t)  is  known  as  the 
zero-input  response,  and  the  remainder  of  y(t)  is  the  zero-state 
response. 

Continuing  with  definitions,  the  m  x  r  matrix  whose  k   column 
represents  the  response  of  a  system  in  the  zero-state  (x  =0),  to  a 
unit  impulse  at  time  t'  at  the  k   input,  all  other  inputs  remaining 
identically  zero,  is  called  the  impulse  response  matrix.   The  zero 
state  response  is  then  given  by 


t 
^zs(t)  =  J    H(t,t')u(t')  dt' 

Co 


which  then  results  in 

H(t,t')  =  C(t)T(t,t')B(t')  +  S(t  -  t')D(t'),  t  >   t'  . 

H(t,t')  is  assumed  to  be  zero  for  t  <  t'. 

Since  the  system  studied  in  this  work  yields  linear,  constant 
coefficient  differential  equations,  the  concept  of  time-invariance  is 
presented  here.   A  system  is  time-invariant  if  its  input-output 
relations  do  not  change  with  time,  i.e.,  the  shape  of  the  response  to 
an  input  depends  only  on  the  shape  of  the  input,  and  not  on  the  time 
the  input  is  applied.   A  system  characterized  by  linear  differential 
equations  with  constant  coefficients  is  time-invariant [ 38] .   The 
property  of  time-invariance  allows  considerable  simplification  in  the 
upcoming  discussions  and  will  henceforth  be  assumed  unless  otherwise 
indicated. 
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Time  invariance  removes  the  time  dependence  of  the  matrices  A,  B, 
C,  and  D  in  equations  B.l  and  B.2.   As  a  consequence,  the  state 
transition  matrix  is  neatly  given  by 

T(t,t0)  =  exp  A«(t  -  tQ) 

Substitution  of  this  result  into  the  state  transition  function  x(t)  and 
the  output  y_(t)  provides  explicit  expressions  for  these  entities. 
Also,  the  impulse  response  matrix  now  depends  only  on  t  -  t',  and  is 
customarily  written  as  H(t  -  t'),  or  H(t)  when  t'  =  0. 

There  are  several  methods  of  computing  exp  A«(t  -  t  ),  including 
finding  a  similarity  transformation  to  reduce  A  to  its  Jordan  canonical 
form[39],  as  a  particular  case  of  function  of  a  matrix,  and  using  the 
Laplace  transform[37,40] .   The  result  is,  taking  t'  =  0, 

T(t,0)  =  exp  At  =  R±   exp  (X±t)  B.5 

where  the  R^  are  the  (constant)  residue  matrices  of  T(s),  the  Laplace 
transform  of  T(t,0),  corresponding  to  the  eigenvalues  \j   of  the  matrix 
A.  This  formula  is  valid  for  the  case  of  distinct  eigenvalues.   If  this 
is  not  the  case,  modifications  are  required  and  reference  is  made  to 
[40]. 

The  definition  of  stability  for  a  linear,  time-invariant  system 
may  now  be  easier  to  understand.   The  response  of  the  system  given  by 
B.l  and  B.2  to  any  input  is  the  sum  of  a  term  having  the  form  of  the 
input  plus  a  terra  containing  terms  of  the  form  of  B.5.   If  the  \i    all 
have  real  parts  less  than  zero,  then  if  the  input  is  bounded,  then  the 
response  will  be  bounded,  as  the  R^  exp  (X^t)  terras  will  decay  with 
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time.   In  the  case  of  the  zero-input  response,  which  is  of  interest  in 
the  present  work,  a  system  starting  with  the  condition  x  =  0  will 
return  to  the  x^  state  following  a  perturbation,  if  the  system  is 
stable.   If  an  eigenvalue  is  identically  equal  to  zero  (real  part),  an 
additional  constraint  is  that  it  be  a  simple  zero,  i.e.,  it  have  a 
distinct  eigenvector  associated  with  it.   Also,  to  prevent  states  from 
growing  unbounded,  any  undamped  modes  of  the  system  (eigenvalue  with 
zero  real  part)  must  not  be  coupled  directly  to  the  input. 

As  in  the  case  of  classical  control,  if  the  system  performance 
does  not  meet  the  desired  specifications,  compensation  is  employed  to 
yield  the  desired  transfer  function.   In  modern  control  theory,  the 
concent  involves  the  utilization  of  all  the  system  states  to  provide 
the  necessary  improvement  in  system  response.   In  effect,  the  classical 
control  theory  design  specifications  are  applied  to  the  system,  but 
these  are  satisfied  by  the  use  of  state  variable  feedback.   A  summary 
of  the  design  procedure  is  as  follows  [33] : 

1.  From  the  state  equations,  the  state  variable  block  diagram  is 
drawn. 

2.  The  equivalent  feedback  transfer  function  and  the  overall 
transfer  function  are  found  and  written  in  terms  of  k-,  the  feedback 
gain  for  each  state  variable. 

3.  The  desired  closed-loop  transfer  function  is  synthesized  from 
the  desired  performance  specifications. 

4.  The  closed-loop  transfer  functions  of  steps  2.   and  3.   are 
equated  and  solved  for  the  k^  yielding  the  necessary  feedback 
compensator. 
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Before  proceeding  to  the  optimal  control  methods,  the  concepts  of 
controllability  and  observability  for  a  time-invariant  system  are 
presented. 

The  system  represented  by  equations  B.l  and  B.2  is  said  to  be 
controllable  if  any  (xo,tQ)  can  be  transferred  to  any  (xpt,), 
including  zero  states,  by  some  input  function  applied  from  t   to  ti. 
For  the  time-invariant  case,  the  test  for  controllability  requires  that 
the  n  x  nr  matrix  Q  defined  by  [B  |  AB  |  •••  |  An~^3]  have  rank  n.   A 
system  is  said  to  be  observable  if  no  x  yields  a  null  zero-input 
response.   The  test  for  observability  requires  that  the  nm  x  n  matrix  R 
=  [C,  A'C',  •••,  (An~  )'C]  have  rank  n,  where  the  prime  denotes 
transpose. 


B.4.   Optimal  Control  Theory 

Optimal  control  is  a  particular  branch  of  modern  control  that 
attempts  to  design  a  system  that  is  the  "best"  possible  system.   Many 
optimal  control  problems  do  not  have  easily  obtainable  solutions,  but 
nearly  all  linear  optimal  control  problems  do.   The  results  of  a  linear 
optimal  control  problem  are  often  directly  applicable  to  nonlinear 
systems  operating  on  a  small  signal  basis,  and  the  procedures  followed 
to  obtain  the  linear  optimal  design  can  often  be  carried  over  to  the 
nolinear  problem.   The  linear  optimal  control  design,  in  fact,  can 
tolerate  nonlinearities  without  major  reduction  of  its  properties.   In 
the  remainder  of  this  section  the  linear  optimal  control  problem  will 
be  discussed  in  the  context  of  the  regulator  problem. 
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"Best,"  as  used  in  the  above  paragraph,  is  defined  by  specifying  a 
performance  index  (PI).   The  selection  of  the  PI  is  often  based  on 
mathematical  convenience.   That  is,  the  PI  is  chosen  so  that  the  system 
will  have  a  mathematical  solution.   The  resulting  optimal  control  will 
satisfy  the  performance  index,  but  certain  classsical  specifications 
may  not  be  met.  k   compromise  must  therefore  be  achieved  between 
desirable  system  characteristics  and  the  tolerable  computational 
difficulty  associated  with  the  PI. 

Given  a  starting  time  and  state,  and  a  final  time  and  state,  an 
integral  performance  index  to  yield  an  optimal  control  is  specified  by 


PI  =  /  f  L(x,  u,  t)  dt. 


L(x,   _u»  t)  is  selected  to  be  positive,  making  PI  a  monotonically 
increasing  function  of  t  in  the  interval  t^  <  t  <  tc.   If  the  magnitude 
of  PI  is  a  minimum  over  the  time  interval,  in  transferring  the  state 
from  x.( t^ )  to  x(ti) ,  the  system  is  said  to  be  optimal. 

Adding  constraints  on  the  amount  of  energy  that  can  be  spent  as 
input  (the  magnitude  of  the  control  functional),  and  the  error  or 
deviation  permitted  in  the  system  states  results  in  the  quadratic 
performance  index  (QPI),  having  the  form 


QPI  =  /   (x'Qx  +  u'Zu)  dt. 
o 


In  general,  Q  and  Z  must  be  nonnegative  and  positive  definite  symmetric 
matrices.  In  addition,  the  system  in  question  must  be  controllable  and 
observable. 


-245- 
The  resulting  optimal  control  law  _u  (t)  that  satisfies  the  QPI  has 
a  very  desirable  property.   The  minimization  of  the  QPI  turns  out  to  be 
achieved  by  a  linear  feedback  law  of  the  form  _u(t)  =  K'^(t)  where  K  is 
the  feedback  matrix.   The  derivation  of  the  optimal  control  can  be 
found  in  the  literature  [41].   The  basic  steps  can  be  outlined  as 
follows.   The  minimum  performance  index  is  shown  to  exist  and  be  of  the 
form  _x'  ( t)Pjc(  t)  ,  P  being  a  positive  definite  matrix  (symmetric).   P  is 
shown  to  satisfy  the  matrix  Riccati  equation,  a  nonlinear  differential 
equation  of  the  form 


-P  •  PA  +  A'P  -  PBZ  ^'P  +  Q. 


*        —  1 
And  finally,  the  optimal  control  is  shown  to  be  u_  (t)  =  -Z   B'Px(t). 

The  notion  of  optimality  can  be  satisfactorily  extended  to 

accommodate  classical  design  criteria,  such  as  closed-loop  pole 

placement,  gain  and  phase  margin,  etc.   The  satisfaction  of  the  desired 

system  specifications  is  achieved  by  the  selection  of  the  Q  and  Z 

matrices  in  the  QPI.   Methods  for  choosing  the  elements  of  the  Q  and  Z 

matrices  can  also  be  found  [42], 


APPENDIX  C 
DERIVATION  OF  SPECIFICATIONS  FOR  A  STABLE  CLASSICAL  SECOND  ORDER  SYSTEM 

WITH  A  FINITE  ZERO 


Given  the  transfer  function  of  a  stable  second  order  system  of  the 
form 


2  2 

_,  .      — — (s  +  z)        — — (s  +  z) 
G(s)  .  T(s)  _     z =  z ca 

R(s)   s2  +  25uns  +  co2n   (s  +  a)2  +  u)2d 


where, 

-z  =  location  of  the  zero 

£  =  damping  ratio 

co  =  undamped  natural  frequency 

a   =  £o)  =  damping  factor 

0)^  =  un/l  -  5    damped  natural  frequency 

analytic  expressions  for  the  following  can  be  derived; 

1.  Impulse  response 

2.  Step  Response 

3.  Time-to-Peak,  t 
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4.  Peak  Overshoot,  M 

5.  Settling  time,  tg  (5%) 

6.  Bandwidth,  u)g 

7.  Resonance  frequency,  (M 

8.  Resonance  peak,  >L 

9.  Phase  margin,  PM 

Note  that  the  restriction  of  stability  translates  into  negative 
roots  of  the  characteristic  equation.   This  in  turn  specifies  £  >  0  and 
oin  real  and  positive. 


C.l.   Impulse  Response 

For  the  case  of  a  unit  impulse  input,  R(s)  =  1,  the  expression  for 
G(s)  can  be  partial  fraction  expanded  and  then  inverse  Laplace 
transformed  into  the  time  domain.   The  result,  given  in  most 
compilations  of  Laplace  transform  pairs  is  given  by, 


cl),//(z  -  a)2  +  m  2 


g(t)  =  e  atsin(u)Ht  +<(,),   1  >  5  >  0   C.2 


2 

[(z  -  a^e"0!6  -  (z  -  a2)e-02t],   g  >  1  C.3 


z(tJo  ~  Oi) 


2   °1' 


with 


=  tan 
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-1  wd 


a,  =  a  +  u>,- 


z  -  a 
* 
'1  -  °  T  ud 

* 

02   =  o   -  w^ 

*     /  2 


C.2.   Step  Response 


For  R(s)  =  — ,  the  expression  for  G(s),  following  the  same 
s 

procedure  as  for  the  impulse  response,  transforms  in  the  time  domain  to 


;(t)  =  1  -  —  zp  e  atsin(w,t  +  d>) ,   1  >  E,   >  0         C.4 

zw,  a 


2za^wj  2Z02WJ 


u_  (z  -  a,)        ul  (z  -  a,) 


with 


-1  ud       -iud 
o  =  tan  +  tan  — 

z  -  a  a 


9      o 
zp  =  /(z  -  a)  +  a)j'c   (a  measure  of  the  pole-zero  separation) 


and  Oi,  0o  and  u)j  as  before. 
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C.3.   Time- to- Peak,  t 

To  determine  the  time  to  peak  overshoot  following  a  step  input, 

take  —2 in  equations  C.4  and  C.5  and  equate  to  zero.   For  C.4,  this 

dt 


gives 


~^-=  0     =>     u^e-^VcosUjtp  +  *)   -   ae-^VsinUrftp  +   <),) 


Rearranging  gives 


d 
tanCcjjt..  +  i>)   =  —  C.6 


Substituting  the  expression  for  <j>  given  previously  for  a  step  response 
into  C.6,  and  defining 

¥,  =  tan-1 —  =>  —  =  tan(kTT  -  f, ) 

L  -a  a  l 

and 


-1  "d 
y0   =  tan  x 


an  expression  for  t  can  be  extracted  as 


kir  -  V? 
t   =  ,   1  >  5  >  0  C.7 


It  can  be  shown  that  the  peak  overshoot  occurs  for  k  =  1  or  2. 


-250- 
Proceeding  in  the  same  fashion  for  C.5  results  in 


1       z  -  0i 

t  =  -i—  In  i),  5  >  1  C.8 

P   2w/    z  "  a2 


Since  the  logarithmic  function  exists  only  for  positive  argument,  a 
necessary  restriction  for  t   to  exist  is  that  the  zero  not  be  located 
between  the  poles,  i.e., 

z  >  Oi  and  z  >  an 
or        z  <  Oi  and  z  <  oo. 


C.4.   Peak  Overshoot,  M^ 


Peak  overshoot,  defined  by 


Mp  =  g(tp)  -  1 


can  be  obtained  by  substituting  C.7  into  C.4  or  C.8  into  C.5,  as 
appropriate.   For  1  >  E,   >  0,  substituting  C.7  into  C.4  gives 


CO  r  0  /  \   1 

g(tD)  =  1  +—  zp  e_L— Ct  "  VJsiny,,   k=l         C.9 
P        zu^        d  l 


=  l  -  ^_  W  e-[^2lT  "  Vising,   k=2        CIO 
zoo.        d  A 


Noting  that  for  1  >  E,   >  0,  sinf,  =  ^1  -  g     >  0,  either  expression 
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define  the  peak,  depending  on  whether  the  zero  is  in  the  left  (z  >  0) 
or  right  (z  <  0)  hand  plane.   The  resultant  expression  for  peak, 
overshoot  is  then, 


Mn  =  ^P-e-br^  "  *2>]  C.ll 

P   |  z  |    wd 


with  k  =  1  for  z  >  0,  and  k  =  2  for  z  <  0. 

For  £,   >    1,  no  clarification  is  obtained  by  substituting  C.8  into 
C.5.   The  expression  for  peak  overshoot  will  simply  be  written  as 


Mp  =  -4 K  e  l   P  "  -7 ^T  e~°2  P  C'12 


with  t  as  given  by  equation  C.8. 


C.5.   Settling  Time,  t  (5%) 


Since  the  damping  factor  a,  and  the  oscillation  frequency  o)j,  do 
not  change  due  to  the  presence  of  the  zero,  the  settling  time  for  the 
present  case  remains  unchanged  from  the  case  of  the  simple  second  order 
system  with  no  zero.   The  settling  time  in  either  case  is  given  by 


3 
t„  =  -  C.13 
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C.6.   Bandwidth,  u)g 

The  frequency  function  associated  with  G(s)  is  obtained  by 
substituting  s  =  jui  into  C.l.   Doing  so  results  in 


Un2 

(z  +  jw) 

G(jm)  =     Z   C.14 

un  -  u/  +  25cona)j 


The  bandwidth  frequency  is  defined  as  that  frequency  at  which  the 

•2 

magnitude  of  the  transfer  function  is  -3dB,  or  —  of  its  zero  frequency 

value.   Applying  this  definition  to  the  magnitude  of  G(jcj)  yields 


^-(z2  +  u2)l/2 

|G(jw)l  = =  —  C.15 

[(u>  2  -   w2)2  +  (2au))2]^2    2 


which  results  in  the  quartic  equation 


u>4  +   [4o2-   2un\l  +  JL-)}*2  -   ^   =  0 


Imposing   physical  constraints   on  oig  results    in  there   being  one  and  only 
one   solution   for  Wg  given  by 


iB  =    (-A  +  /A2  +  u>nV/2  C.16 


with 
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2 


A-  2a2-%2(l+-^- 


C.7.   Resonance  Frequency,  u^ 

The  resonance  frequency  can  be  found  by  taking  the  derivative  of 
the  magnitude  of  the  frequency  function,  equation  C.15,  equating  the 
result  to  zero  and  solving  as  follows: 

d|G(ju>)|  m   0  => 
dw 

=>     w[(wn2  -  u>2)2  +   (2au>)2]1/2(z2  +  w2)~1/2  = 

[4a2u)  -   2a)(a)n2  -   u>2)][Un2  -   u>2)2  +   (2aa))2]"1/2(z2  +  u>2)1/2 

This   reduces    to   the  quartic 


j4  +  2z2w2  +  4a2z2  -   2wn2z2  -   ain4  =   0 


which  physical  constraints   limit   to  at  most    the   one   solution  given  by 


u)m  =   {-z2  +   [z4  +   (2u>n2  -   4a2)z2  +  ^ ]  l/2}l/2  c.17 


-254- 
C.8.   Resonance  Peak,  M 

No  gain  in  clarity  is  obtained  by  substituting  C.17  into  C.15. 
Therefore  ML  is  presented  as 


co  2(z2  +  um2)l/1 
\- ,  "     =   m , C.18 


with  u  given  by  equation  C.17. 


C.9.   Phase  Margin,  PM 

For  a  system  that  is  stable,  phase  margin  is  defined  as  the  excess 
phase  gain  available  at  the  gain  crossover  frequency  over  that  required 
for  marginal  stability.   For  a  system  that  is  initially  unstable,  phase 
margin  is  defined  as  the  excess  phase  lag  at  the  gain  crossover 
frequency  over  that  required  for  marginal  stability.   To  find  the  gain 
crossover  frequency,  set  the  magnitude  of  the  frequency  function, 
equation  C.15,  equal  to  one  to  obtain 


ci)  4 
'go"  (^T+2"n2-^2)1/2 


The  phase  of  the  frequency  function  at  to   is 


/G(Ju„r)  =  tan  1-ii-  tan  X- 


S£_        z     .  2      2 


so  that  the  phase  margin  is 
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PM  =  180°  -  |/G(jugc)|,   /G(jugc)  <0  C.19 

=  /G(jcogc)  -  180°,   /G(jugc)  >  0  C.20 
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